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Summary. — A tentative model, a modification of Konopinski-Mahmoud’s 
model concerning the leptons is developed for the classification of the 
processes where leptons take part. By means of this model, we obtain 
a good selection rule, and also the redundant neutral particle has not 
been introduced which always appears in the model based on the charge 
independence hypothesis extended to leptons. We have referred to the 
relation between the N-G model and our scheme and shortly discussed 
the interaction structure with this scheme. And though we obtain a good 
selection rule for lepton processes, we conclude that for a classification 
of leptons, their interaction structure is always to be taken into account. 


1. — Introduction. 


To interpret the contradiction between the copious production of the new 
particles and their metastabilities, several ways (1) have been proposed. 
Their common idea consists in the introduction of new quantum numbers 
by which the appropriate selection rules have been expressed. In this way, 
GELL-MANN and NISHIJIMA (1) have presented a scheme for classifying the 
fundamental particles which belong to the baryon- and meson-families under 
the charge independence hypothesis and characterized these particles by means 


(*) Preliminary reports of the present paper have been published in Prog. Theor. 
Phys., 17 (1957). 

(4) A. PAIS: Physica, 19, 869 (1953); Proc. Nat. Acad. Sci. U.S., 40, 484 (1954); 
M. GeLL-MANN: Phys. Rev., 92, 833 (1953); T. NAKANO and K. NISHIJIMA: Prog. Theor. 
Phys., 10, 581 (1953); K. NISHIJIMA: Prog. Theor. Phys., 13, 285 (1955); B. D'ESPAGNAT 
and J. PRENTKI: Nuovo Cimento, 3, 1045 (1956). 
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of the isotopic spin and strangeness. Thus s far as these particles are concerned 
they could explain successfully the phene:: ‘na concerning strongly interacting 
particles and discriminate whether the in’ ractions are strong or weak. 

On the other hand, from the phenomenc.ogical point of view, three kinds 
of interaction have been found among the baryons, mesons and leptons; 
namely, strong, electromagnetic and weak interactions. Among these three 
kinds of interactions, as to weak interactions, there seems to exist the regu- 
larity relating to all the families; baryon-, meson- and lepton-families, as. 
pointed out by OGAWA (?).. This suggests that the weak interactions have a 
more fundamental character than the strong interactions. Therefore it is 
desirable to formulate a theory of weak interactions governing both the baryons 
and leptons. For this purpose it is desired to introduce some new quantities 
by which the appropriate selection rules can be expressed retaining the GELL- 
Mann-Nishijima scheme. 

There seems to exist two ways of approach for such programs. The first 
is to consider the isotopic spin as having some intrinsic significance for all 
the fundamental particles including leptons (*). This may be justified from 
the fact that in the electromagnetic interaction the charge operator Q=J,+« 
associates always with the 3rd component of the isotopic spin. The second 
is to seek some other new quantum numbers not based on the charge inde- 
pendence hypothesis. This standpoint, which takes the extension of the charge 
independence hypothesis to the lepton family as to be wrong, will be supported 
from the following facts: 1) no strong interactions have been found in the 
transition processes where the leptons take part; 2) the ratios of mutual mass 
difference between the leptons are large and 3) the weak interactions involving 
leptons are weaker than the electromagnetic interactions which violate the 
charge independence. It might be required in this approach to treat the lepton 
processes as dynamical problems rather than kinematical ones, leaving a deli- 
cate problem whether those processes should be forbidden by some selection 
rules or by the detailed structure of the interactions concerned. 

In this note we would like to discuss on a tentative scheme for the classi- 
fication of the fundamental particles including leptons taking into account 
their interaction structures. The first approach is abandoned because it re- 
quires the introduction of the redundant neutral particles, though it is the 
simplest way for the classification and produces good selection rules. We take 
the second approach and assign the 7,-charge to the leptons in an analogous 
manner as in the N-G model in addition to the conservation of the lepton number 
proposed. by KONOPINSKI-MAHMOUD (*). In our approach useful selection 


(2) S. Ogawa: Prog. Theor. Phys., 15, 487 (1956). 

(3) T. OucHI and K. SENBA: Prog. Theor. Phys., 16, 528 (1956); S. Gord: Prog. 
Theor. Phys., 17, 107 (1957); Y. Karayama: Prog. Theor. Phys., 17 (1957). 

(4) E. J. Konopinski and H. M. ManmouD: Phys. Rev., 92, 1045 (1953). 
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rules can be obtained by means of 7,-charge and the difficulty concerning the 
redundant neutral particles can be overcome. 

The conservation of the lepton number has not so clear physical meaning 
as that of the baryon number which relates to the stabilities of matter satis- 
fying the experimental requirement since we take the lepton as Dirac particels. 
However, if it happens that the neutrino is proved to be a Majorana particle, 
our model will be ruled out. Fortunately the analysis carried out by F. REINS 
et al. (5), seems to support our assumption that the neutrino is a Dirac part- 
icle. Hereafter, for the sake of simplicity, only one kind of neutrino is assumed. 

In the following sections, the selection rules by means of 7,-charge will 
be presented so that 


1) |An.|=2 for the processes: 
[baryon] + [meson] + [baryon] + [meson] + [lepton] ; 
2) Am, = 0 for the processes: 
[lepton] — [lepton ]. 


We discuss in detail about the lepton processes and also shortly refer to 
the interaction structure which causes the mutual transition between the 
various families. 


2. — Assignments of a new quantum number for the lepton family. 


About the lepton processes several attempts (5%) have been presented 
-assuming the existence of Fermi interactions among four fermions. But they 
cannot forbid the unwanted processes completely though these processes are 
consistent with the various conservation laws. To avoid this defect YANG and 
TIOMNO (6) tried to introduce the intrinsic « generalized parity » in such a way 
that the generalized parity conservation is violated in the unwanted processes. 
CAIANIELLO (7) reported that all the unwanted processes can be successfully 
eliminated if the coupling constants are permitted to change in a certain way 
corresponding to the changes of co-ordinates. And TIOMNO (*) proposed to 
choose the suitable combination of coupling types by imposing the « mass 
reversal » invariance on the Fermi interactions. 


F. Reins and C. L. Cowan Jun: Nature, 178, 446 (1956). 

C. N. YANG and J. Tiomno: Phys. Rev., 79, 495 (1950). 

E. R. CAIANIELLO: Nuovo Cimento, 8, 534 and 749 (1951); 9, 336 (1952). 
J. TroMmNO: Nuovo Cimento, 1, 226 (1955). 
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In comparison with these attempts Konopinski-Mahmoud’s model (4) is 
quite interesting. They showed that the unwanted processes can be forbidden 
by the conservation of lepton numbers if the positively charged u-meson, 
the negatively charged electron and the neutrino are assumed to be normal part- 
icles. With the accumulation of the experimental facts for new unstable 
particles, the K-M model has been found insufficient to exclude the unwanted 
processes concerning strange particles. 

To complete the K-M model, let us assign the 7,-charge to leptons under 
the hypothesis of the conservation of the lepton number in an analogous 
manner as in the Nishijima model, assuming the relationship. between electric 
charge @ and 7,-charge as follows: 


(1) Q =43(n +m). 


The 7,-charge was assigned to each particles in the following table. 


TaBLE I. — Types of fields. 


Name Lorentz space Q n mM 
U spinor Il 1 1 
e » == ll 1 — 3 
v » 0 1 | —l 


In order to avoid the introduction of the redundant neutral particles appearing 
in the charge-independent model of leptons, the leptons are assumed here to be 
scalar quantities in the isotopic spin space. An anti-particle, i.e., the charge 
conjugate particle, has an opposite 7,-charge to that of the normal particle as 
is the case for the electric charge, and we have 


(2) MEA, mC = — Cm, 


where C is the charge conjugation operator. 
Under these assumptions we shall go on into the detailed investigation of 
lepton processes in the following sections. 


3. — Interaction structure and the classification of fundamental particles. 


As one of the phenomenological approaches to the investigation of the 
structure of interactions, it seems to be most attractive to classify the ele- 
mentary particles into several families 1) baryon-family, 2) lepton-family, 
3) meson-family and 4) photons. 

Among these families the electromagnetic field (photon) interacts with all 
the charged particles with a universal coupling constant «e», satisfying the 
requirement of gauge invariance. Therefore we take this constant «e» to 


cà 
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characterize the electro-magnetic field « A, ». Moreover, for the electromag- 
netic field there exists an invariance principle under gauge transformation, 
but no invariance principle like this has yet been found for other families. 
However there seems to be some peculiar regularity concerning the weak inter- 
action among these families as has been recently stressed by OGAWA (2). The 
coupling constants of these interactions take a unique value of the order of 
magnitude (10-14) in the natural unit regardless of coupling type or the kinds 
of particles concerned. Because this weak interaction occurs not only among 
the baryon- and meson-families, but also among the lepton family, we may 
ask why the particles belonging to all these families have a coupling con- 
stant of an unique order of magnitude, and what is the source of weak 
‘interaction. This question will offer a key to the study of the interaction 
structure among these particles. Further, the fact that the weak interaction 
associates with all the families suggests that the weak interaction has more 
fundamental and elementary character than the strong interaction which appears 
only between baryon- and meson-families. 

It may be a clue to the above question to remind the fact that the lepton 
processes involve at least one neutrino. But with the accumulation of the 
experimental data about new unstable particles, there exists the other weak 
interaction that does not involve neutrino. So, we cannot take the neutrino 
as the unique origin of the weak interaction. On the other hand, on the belief 
that all the weak interactions come from a single base, OGAWA has proposed 
an attractive model in which new unobserved charged Bose fields (simply 
called « B-field ») have been introduced and are assumed to interact with all 
Fermi particles with a unique coupling constant. His fundamental idea con- 
sists rather in the introduction of the B-field which plays a similar role as 
the electromagnetic field does in 
the electromagnetic interaction. 

From the phenomenological 


AN=0, AI-0,AI,-0, 
At,=0 , An=0 + Ate=0 


baryon family 


investigation we know as uni- e (ee . IATET.AN:0) 271 
+ . Q | # = = n =~ 
versal weak interaction the 81,3 SA lARF; 
A è : : ira 1At,1>0 | Se a aa Le [Azle 1 
following three kinds; 1) wea AN:0 = (| AW: 0 
; . ers [Azl= 2] | DIA > [Agl= 2 
fermion-boson interaction, 2) ASA e Rss lee 
universal weak Fermi inter- Usi 3S 


N 
action extended to baryons | ep tontionniy, L A ee lepton , family 
(pv) D ann (e.v) 


and 3) weak boson-boson inter- At,:0 . At:0, An:0 
action. We assume in this —— Strong interaction 
paper for the sake of simplicity ————— Weck interaction 
that the interactions occur Eig:ol. 


through the universal weak 
Fermi interaction. Then the transition scheme among these families is as 


follows (Fig. 1): 
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Analogously to the N-G scheme for the strong interactions, we assume another 
hypothesis than the charge independence for the fundamental weak inter- 
actions. Though it may be more reasonable to treat both the weak interaction 
between the baryon- and meson-families and the weak lepton processes on 
the same footing, we are forced to treat them differently because of the 
peculiar character of lepton processes. But 7,-charges of leptons are shown 
to produce a good selection rule. 

Now we shall go further to investigate the transition processes from the 
view-point of family mentioned above as follows: 


1) The strong interactions; [baryon] and/or [meson] + [baryon] and/or 
[meson]. 
Since these processes are strong and not concerned with the leptons, we 
do not enter into a detailed discussion here. The isotopic spin J, its 3rd 
component J; and ,-charge are conserved in this case. 


2) The weak interactions: [baryon] or [meson] [baryon] and/or[ meson]. 
This processes are governed by the rule 


(3) IADIÆE, MAN = 0 and.  (|An,;|21. 


For these processes it will be an interesting problem to examine whether 
these are due to direct interactions or to indirect ones. Here the indirect pro- 
cesses are meant to be the processes that occur through the intermediary of 
the universal weak Fermi interaction extended to baryons. 


3) The weak interactions: [baryon] and/or [meson] +[baryon] and/or 
[meson] + [lepton]. 

For this case the main task is to exclude the processes which are consistent 
with the conservation law, but are not observed experimentally: several at- 
tempts (57) to explain these difficulties have been published so far. Our at- 
tention will also be concentrated on such processes. 

We will here give a scheme to rule out the unwanted processes for leptons 
by the 7,-charge rule and the lepton-number conservation law instead of the 
charge independence hypothesis without introduction of the redundant neutral 
particles. Thus the following scheme is assumed 


(4) CATANIA MA NG 0", [Am |= 0 An =0 and |An;|=2 


for the process: [baryon] > [baryon] + [lepton], 


(5) (AL = AN = 07,77 An, awd An =0 and |Ay,| = 2 


for the processes: [meson]— [meson] + [lepton]. 
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4) The weak interactions: [lepton] + [lepton]. 
Only one process belonging to this case is known so far (u— e + v + v). 
For this process, the following assumption will be adequate 


(6) An = 0 and Nj 0 


It is interesting to investigate the interaction structure of this process in con- 
nection with the problem of the elementarity of the interactions. Compared 
with the case 3) this assumption seems somewhat unreasonable. Perhaps it 
may be due to the difference of the interaction structure of leptons from that 
of baryons. 


4. — Allowed lepton processes. 


In this section we would like to discuss the processes which are allowed 
by our model and consistent with the conventional conservation laws. As is 
referred in Sect. 3 the following characteristic features are recognized experi- 
mentally for the phenomena involving leptons: 1) The lepton processes involve 
at least one neutrino, 2) they take place always in pairs (u,v) or (e, v) 
and 3) in two body decay u and e appear only in an asymmetrical manner. 

The established lepton processes so far are as follows (°): 


(7) ND a ee Yc 
(8) PT - N +5, 
(9) DAC VIVI, 
(10) TV, 

(11) K},—ut+5, 
(12) K',—=n+ut+5, 
(13) KERI 


The processes except the decay of heavy mesons have been analyzed on the 
hypothesis that the elementary interactions are strong interactions (z-N) 
and that weak universal Fermi interaction occurs not only among nucleon 
and lepton, but also among leptons themselves, and that reasonable results 
have been obtained. As for the process (12), the accumulated experimental 
facts indicate that one of the two neutral decay-products is 7°. This process 


(°) See, for example, A. M. SHAPIRO: Rev. Mod. Phys., 28, 164 (1956). 
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might be taken as an alternative mode of decay of the same K particle decaying 
according to the reaction Ki > x*+7°, that is the K-meson decays virtually 
into xt+7° and then the charged pion decays into a u+v pair. Even if this 
turned out to be acceptable for the K,,, a similar explanation can not be 
applied to the K,, because the decay of a pion into an electron is highly forbid- 
den so that the K,,-decay should be much less probable than K,,, while expe- 
rimentally K,, are found as often as K,,. Therefore we will assume that bosons 
interact through virtual fermion pairs and that universal weak Fermi 
interactions exist not only between nucleons and leptons, but also between 
hyperons and leptons. Accordingly K,,- and K,,-decays occur through the 
baryon family loop shown as follows (1°): 


. ZT 7 a 
Z 
f F 
v vy 


where G is the strong coupling constant (1 ~ 10-1) and f is the weak coupling 
constant (~ 10-1). We get in this way for the transition probabilities (11), 


nmax 
0.044 1/ @ \fMic (eh ((1_. (9 
Cru 12898 (amie) ho Meu e] 4" pie 
(14) 0 


le so 
+ apr + aa 


where Q is the Q-value of the K reaction and M is the mass of the baryon 
and their ratios to the decay probability of the K_, are given by 


max 


MO CRC E CO: 1 
ORs _ n Arc w2 vito, \u2 | 11 pa je Bi! LT Sf 4 
(15) Keg 


OK OK 
0 E 
WE 3 WK, 


where v, and v, are the velocities of the two final pions of K_, respectively. 
These results are in good agreement with the experimental data. A similar 
consideration can be also applied to K,,-decay. Although the perturbation 
calculation may not be adequate for the strong interaction, these consideration 


(1) S. Fururicur, T. Kopama, 8. OGAWA, Y. SUGAHARA, A. WAKASA and M. YONE- 
ZAWA: Prog. Theor. Phys., 17, 89 (1957). 
(11) G. Costa and N. DALLAPORTA: Nuovo Cimento, 2, 519 (1955). 
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might be sufficient to give an indirect test concerning the elementarity of the 
weak Fermi interaction. 

If we assume the generalized weak Fermi interaction as the elementary 
one, the 7,-charge selection rule of our model allows the following lepton pro- 
cesses besides the transition processes (7)-(13). Here the processes are re- 
stricted to two or three body decays (12). 


1) The transition processes due to the generalized weak Fermi inter- 
action : 


[ A p+e-+5, 
(16) ta AP DE pr se; 
[ca 5 


and the charge conjugate of these processes. The analysis of these proces- 
ses has been done basing on the generalized weak Fermi interaction by 
Iwata et al. (1%), and they obtained the life times for the above processes 


(17) t(A°> p+e-+ 9) ~ (10-8 = 10-19) s, 


showing the detection of them as possible in future. And their detection will 
give a clue for the investigation of the interaction structure. 
2) The transition processes occurring through the combination of the: 


strong and the generalized weak Fermi interaction: 


| K°> xz +er+y, 


(18) | K°> 7 +ut+y, 


Though these are not yet firmly established, they are expected to be observed 
with frequency comparable to the K,, > r°+e+v and Ke T°+u+v because 
they can be considered as the neutral counter particle of the K. 


3) The transition processes through à combination of the strong and the 


(12) This restriction would be reasonable because the matrix elements of the other 
transition processes, containing the weak elementary interactions twice or more and 
having small phase volume, should be much smaller than the usual observed weak 


processes. 
(38) K. Iwata, S. Ocawa, H. Oxonoai, B. SAKITA and S. ONEDA: Prog. Theor. 


Phys., 18, 19 (1955). 
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generalized weak Fermi interaction: such as (14) 


| T > @e+v 
(19) 
| K+>e+y (5). 


These are not yet established too and they seem to be highly forbidden. So 
far as we take the strong and the generalized weak Fermi interactions as the 
elementary ones, it may be difficult to forbid these processes by means of the 
selection rules only. Perhaps it may require the consideration of the dynamical 
structure of mutual interactions. The investigation along this line has been 
done by several authors (1%) taking as the effective hamiltonian due to the 
strong and weak Fermi interaction 


(20) MP, or () VP) Pr Ga C. C. ? 


where M', the contribution from various intermediate states, is a function of 
momenta and masses of the intermediate particles and y, is y; or I. They 
showed the ratio of the 7-e-decay probability to that of 7-u-decay 


R = œ(r— e +v)/o(r -u+v)— 10", 


which is consistent with the experimental results. Such à consideration may 
be applied to K > e + v. 


5. — Forbidden lepton processes. 


In this section we shall discuss the lepton processes which are forbidden 
in our model, but allowed by the other conservation laws, which is one of the 
most interesting but delicate problems for lepton processes. Though this 
problem has been investigated by several authors, satisfactory results 
have not yet been obtained. The main lepton processes which should be for- 
bidden by the experimental indication are as follows (limited to two- and 
three-body decays): 


(14) To forbid these processes, SHIMAZU tried to assign the lepton number to the 
Bose particles (x, K). H. Sximazu: Soryushiron-Kenkyu (Mimeographed circular in 
Japanese), 13, 194 (1956). 

(5) Though this mode has been inferred by CRUSSARD, it has not definite evidence. 
J. CRUSSARD: Proc. 6th Rochester Conf., V-19. 

(16) See, for example, S. ONEDA and A. WAKASA: Nucl. Phys., 1, 445 (1956). 
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1) three-body decay of heavy meson 


Ker er Eur 
K++ rt+et + e- 
(21) Kt at4+ut+u- 


K+t> 7t4v +9 


—_——$—— T_—— 


Ki rt + utt 


(22) | Kt> rt+u-+et 


2) three-body decay of baryon 


Ap > N +v4+39 
D > N +v+ÿ 


(23) 
MESSE, 


3) two-body decay 


(ec 
(24) | K°— et +e 


Among these processes, Konopinski-Mahmoud, by the conservation of the lepton 
number, could prohibit (22) and u-> e+et+e-, but not (21), (23) and (24). 
For the reaction A°> N +v+%, we get the life-time 


(25) tA°>N+yv+9%) 107 8 


and its ratio to the main decay mode of A°, A°+> p+r- is 


o(A°>N+v+ 3) 


ZANE 
OI 


(26) 


459 
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so the processes (23) can be rejected by dynamical considerations, though they 
are allowed by conservation laws. 

On the other hand, for the reaction K+— r++v+9, we get the following 
life-time 


(27) tK>r+v+7) 21048, 


which is comparable with t(K,,— m+u+yv) and 7(K,— x°+e+v). And 
its ratio to the decay probablity of K_,, the main decay mode of the K-par- 
ticle is, | 


w(K — x +v +) 


w(K > n° + 7) Li 


(28) 


Thus, as for the processes (21), we cannot reject them dynamically, too. With 
our model, all these processes, (21), (23) and (24), can be ruled out by the 
n:-charge conservation rule. 


6. — Concluding remarks. 


In this paper we have discussed mainly the lepton processes based on the 
hypothesis that the Fermi interaction is an elementary weak one, and « uni- 
versality » exists among the interactions concerned with the fundamental 
particles. We have clarified how K-M’s attempt to prohibit the unwanted 
processes concerning leptons fails to explain the fact that processes 
such as K—7m+v+9 are not observed in nature. By introducing a new 
quantum number 7, for leptons, minimizing reasonable assumptions we have 
tried to forbid the unwanted processes and succeeded in prohibiting all the 
unwanted processes including those which were impossible to be forbidden 
by the K-M model. 

But there remain some defects in our model, that is; the reaction 7 > e+ v 
which is indicated experimentally to be highly forbidden could not be excluded 
by our selection rules, and we were also forced to put some unnatural as- 
sumption for the 7,-selection rule to get good coincidence with experimental facts. 

Nevertheless, as is noted in the introduction, whatever standpoint we may 
take, it is necessary to keep in mind that the weak interactions have the 
coupling constant of the same order of magnitude (~ 10-11) for both hyperon 
and lepton. So it may be still worthwhile to seek some scheme to explain all 
the processes concerned with fundamental particles on the same footing, though 
our model has some defects. 
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As for these defects, they may be due to the peculiar dynamical structure 
of interaction of leptons, or they may depend on the form of the interaction 
taken an as elementary one. So it seems to be important to enter into the 
details of dynamical structure of the interaction. 


* % * 


The author would like to thank Prof. K. SAKUMA, Drs. S. OGAWA, T. OUCHI 
and M. YOoNEZAWA for helpful discussions. 


Notes added in proof. 


1. — A similar analysis in which a Konopinski-Uhlenbeck type of weak interaction 
is assumed instead of the Fermi interaction taken in this paper has been performed 
by Y. Miyacui (Progr. Theor. Phys., in press). 


2. — Recently J. J. SAKURAI and J. WERLE independently have proposed a possible 
way to check directly the validity of time reversal invariance in the K,, decay process. 
It can be done by measuring the muon polarization perpendicular to the decay plane. 
If the polarization along this direction is observed the condition (A.5) is violated 
in this process. It is noted that the conditions for the coupling constants used by J. 
WERLE (Nucl. Phys., 4, 171 (1957)) are not correct. The correct oneis (A.5). The author 
would like to thank Dr. J. J. Sakurai and Dr. J. WERLE for sending him preprints. 


TRIASSICO) (9), 


Si sviluppa un modello di prova, modificazione del modello di Konopinski-Mahmoud 
riguardante i leptoni, per la classificazione dei processi cui partecipano leptoni. Per 
mezzo di questo modello si ottiene una buona regola di selezione, senza dover intro- 
durre la particella neutra sovrabbondante che sempre appare nei modelli basati sulla 
ipotesi dell’indipendenza della carica applicata ai leptoni. Si fa in breve riferimento 
alla relazione tra il modello N-G e il nostro schema e si discute al lume di tale schema 
la struttura dell’interazione. E per quanto si ottenga una buona regola di selezione 
per processi leptonici, si conclude che per una classificazione dei leptoni si deve sempre 
tener conto della loro struttura di interazione. 


(*) Traduzione a cura della Redazione 
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On the Threshold Behavior 
of the Negative to Positive Ratio in Pion Photoproduction. 


M. J. MORAVOSIK (*) 


Brookhaven National Laboratory - Upton, N.Y. 


(ricevuto il 31 Maggio 1957) 


Summary. — It is shown that the extrapolation used by BENEVENTANO 
et al. to obtain the negative to positive ratio at threshold is very sensitive 
to small alterations which can be made due to the approximate nature 
of the procedure. Such small alterations bring about very substantial 
changes in the energy dependence of the ratio near threshold and hence 
the extrapolation is quite ambiguous in the absence of extensive experi- 
mental data on the ratio in this energy region. Since such extensive set 
of data does not exist at the present time, it is concluded that there is 
no good reason to believe that the ratio at threshold is in disagreement 
with the theoretical predictions. 


1. — Introduction. 


The value of the negative to positive ratio in pion photoproduction and 
its relationship to other experimental quantities has been somewhat of a puzzle 
for some time. This quantity is an important clue to the S wave part of the 
pion-nucleon interaction, which is considerably less well known than the P_wave 
part. In fact, a consistent scheme which includes the interpretation of the 
threshold ratio in photoproduction, the S wave scattering phase shifts, the 
Panofsky ratio, and perhaps even other quantities, has been the aim of nu- 
merous investigators during the past few years. 


(*) Work carried out under the auspices of the U.S. Atomic Energy Commission. 
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One of the important papers in this field has been that of BENEVENTANO 
et al. (*). This paper, besides presenting some new experimental data, at- 
tempted an over-all analysis of photoproduction data in the energy region 
between threshold and about 240 MeV photon energy in the laboratory system. 
In particular, it derived a semi-empirical equation (reference (1) Eq. (21)) for 
the energy dependence of the ratio near threshold and used the coefficients 
of the angular distribution of positive photopions in conjunetion with four 
pieces of data on the ratio itself to extrapolate and determine the ratio at 
threshold. The resulting ratio turned out to be 1.85 or even higher. It was 
remarked that this value strongly disagrees with the value predicted by the 
threshold theorems and the dispersion approach, which claim that the result 
of the perturbation theory (about 1.4) should be correct. Similarly, the pion- 
nucleon coupling constant obtained by the usual extrapolating procedure also 
turns out peculiarly low when we utilize this high value-of the ratio. The 
value of the coupling constant has been determined with great apparent re- 
liability by several procedures, all agreeing very well, while the coupling 
constant obtained from this high threshold ratio is more than 20% lower than 
the commonly accepted value. On the other hand, it was also remarked that 
the high ratio can be correlated more easily with the latest value of the dif- 
ference in the slopes of the two S wave scattering phase shifts and with the 
latest value of the Panofsky ratio. The implication of BBCST is that the real 
discrepancy is between the theoretical prediction of the ratio and the « exper- 
imental » value of the ratio as obtained by BBCST, and that this « exper- 
imental » value is in reasonable agreement with other clues to the S wave pion- 
nucleon interaction. 

The purpose of this paper is not to solve this discrepancy and to give a 
unified scheme for the S wave interaction. Instead, the aim is to show that if 
one alters slightly the approximations used by BBCST one can obtain a quite 
different and considerably lower result for the threshold ratio. It is not claimed 
that the new ratio is a priori a more reliable value than the old one. It is 
claimed however, that our modification of the BBCST analysis is equally in 
agreement with the general principles used and with the experimental data 
which go into the numerical evaluation. . 

The important criterion for the choice of the extrapolating formula is of 
course the experimental material on the ratio itself. BBCST used four pieces. 
of data at 90° (see reference (1), Fig. 5). It has been pointed out (?) that 
apart from the quantitative inadequacy of these data they also fall at higher 
energies below, and at lower energies above, the value one obtains when inter- 


(1) M. BENEVENTANO, G. BERNARDINI, D. CARLSON-LEE, G. SroPPINI and L. Tau: 
Nuovo Cimento, 4, 323 (1956). This paper will be referred to as BBCST. 
(2) M. J. Moravesix: Phys. Rev., 105, 267 (1957). 
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polating the measured values of the ratio at angles other than 90°. In other 
words, the data used by BBCST give a steeper energy dependence of the 
ratio at 90° than one would deduce from using all data on the ratio. This 
contention is supported by the new preliminary data of DUDZIAK (*) on the 
energy dependence of the ratio near threshold. These data, which are con- 
siderably more extensive than those used by BBCST, fit much better into 
the angular distribution curves of the ratio (?) which have been measured 
previously. The new data seem to give a considerably gentler energy depend- 
ence and if extrapolated « by eye » to threshold they give a considerably lower 
ratio, a value around 1.6. The purpose of this paper is to show that such a 
low threshold value and gentle energy dependence (if confirmed by the final 
version of Dudziak’s data) is by no means incompatible with a semiempirical 
analysis which is of the same general type and of equal reliability as the one 
worked out by BBCST. 

Such a lower ratio, if yielded by the experiments and supported by a 
semi-phenomenological scheme, might easily be in agreement with the theore- 
tical prediction of 1.4. The difference, now only about 10%, could perhaps 
be explained by the deuterium correction which is yet to be given a really 
reliable treatment, but which could certainly amount to that much. 

Of course, on the other end we made things even worse, and the disagree- 
ment with the Panofsky ratio and the S wave scattering phase shifts now 
replaces the disagreement we just eliminated. It should be noted, however, 
that the experimental values of the Panofsky ratio and of the S wave phase 
shifts are still subject to large changes, and that the connection between these 
quantities and the threshold ratio might also be subject to modifications (4) 
Thus, without offering an over-all solution of the general problem of the 
S wave interaction, we will restrict ourselves to pointing out a way to a con- 
sistent picture of the special problem of the threshold value of the negative 
to positive ratio in photoproduction. 


2. — The analysis of Beneventano et al. 


Let us first review the assumptions and the derivation of the BBCST 
analysis. The relevant formula giving the energy dependence of the ratio 
near threshold is Eq. (21) of reference (1). We will follow the notation used 
in that paper. The derivation which leads to the above equation can be sum- 


(3) W. F. DUDZIAK: private communication. I am greatly indebted to Dr. DUDZIAK 
for sending me his preliminary data which aroused my interest in this problem. 
(4) See e.g. H. P. NoyEs: Phys. Rev., 101, 320 (1956). 
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marized as follows: 


(2.1) 


(2.2) 


(2.3) 


(2.4) 


yen 
where a,= 4 


(2.5) 


(2.6) 


(2.7) 


s = RP a 


o a + a; COS O0 + a; cos? 0 
ot at + a; cos 0 + at cos? 0 ? 


=f = rationat 90° , 


ag = (SF R)?+ a, 
af ——2(S R)K, 


0, = K*?*—ay Y — dp; 


p — ao _ WS + RY —a, 
ya (SR) Lg, 


2 8 = (ray + a) + (a+ a, 
2R = (ra + Ay)? — (G+ a»)? ; 


Ie pon Se 
O 
Vo V vo 


S(v, w) 


Si Va, (Vr, + 1) ; 


445 


where 7, is the threshold value of the ratio and » is the photon energy at 


threshold. 
Ne 
(2.8) (ra, + 42) + (do + a2)? = Ga (/% +1), 
Vo 5 Vo \} Ay 
(2.9) r= (2) 0 +ym—2(22) a+ vm) (1+ 2) +2, 


which is identical with Eq. (21) of BBCST. 

Next let us list all the assumptions which enter the above derivation. 
Eq. (2.1) implies that the angular distribution for photoproduction can be 
represented to a good approximation by a quadratic expression in cos 0. We 
know (5) that this might not be so, depending on what feature of the angular 


(5) M. J. Moravcsix: 


Phys. Rev., 104, 1451¥(1956). 
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distribution we are interested in. For instance, BBCST state that a, (which 
gives the cross-section at 90°) is constant in the energy region we are consi- 
dering, namely, from threshold to 260 MeV photon energy in the laboratory 
system. The more accurate quartic analysis (*) gives a 10% change in the 
same energy range. Since our analysis as well as that of BBCST utilizes only 
the ratio a,/a), which is small in any case, the assumption of Eq. (2.1) will 
probably not cause any drastic errors. Nevertheless this is one possible source 
of small corrections which, as we will see, have a very substantial influence 
on the conclusions. 

Eq. (2.2) is the definition of r which, granting Eq. (2.1), is precisely the 
ratio at 90°. 

Eq. (2.3) is based on several assumptions beside the one underlying (2.1). 


a) It is assumed that the quantities on the right hand side are real. 
In view of the fact that at 240 MeV the P,, phase shift is around 15°, this 
assumption might easily introduce an error up to 5% in the P wave para- 
meters on the right hand side. 


b) It is assumed that the so-called «spin flip» amplitude is negligible 
compared to the «no-spin-flip » amplitude. Using Table V and Eq. (A-5) of 
BBCST we can see that this assumption might introduce another error up to 
8% in the magnitude of a. 


c) Itis assumed that the S-D interference term can be neglected. Actually 
BBCST considered the S-D interference but, since they found that it could 
not resolve entirely what they considered was a discrepancy between theory 
and experiment, they omitted it in their subsequent analysis. The inclusion 
of such term might again contribute a change of about 10%. 


d) It is assumed that a, is the same for positive and negative pion, or 
in other words, that the P wave negative to positive ratio (a not too well 
defined quantity) is unity. This is clearly not so, as we know (?) from the 
behavior of the ratio at higher energies. The effect of the P wave ratio at 
lower energies is something which is not too easy to estimate, but another 10%, 
is not an unreasonable guess. 


We see, therefore, that Eq. (2.3) is contingent on several assumptions 
which hold only approximately. We will consider below what happens to 
Eq. (2.9) if we use slightly different values for the parameters, values which 
are equally possible within the approximations we made. 

Eq. (2.4) follows from the previous equations without further assumptions, 
and so does Eq. (2.5). 


(6) M. J. MoRAVOSIK: Phys. fev. 107, 600 (1957). 
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Eq. (2.6) is discussed in detail by BBCST and so we will not repeat the 
arguments which lead to this form of S. It is important to notice, however, 
that BBCST use the approximation in which o is neglected. It is pointed out 
by BBCST that the leading term in the expansion of o into powers of « is 
of the order of (u/M)?. We will consider below how the inclusion of such à 
term affects the energy dependence of Eq. (2.9). 

Eq. (2.7) follows from Eq. (2.5) without further assumptions, and Eq. (2.8) 
and (2.9) follow from the previous ones without further assumptions. 

Summarizing, we see that all assumptions are introduced in Eqs. (2.1), 
(2.3) and (2.8). 

Before we proceed to illustrate the effect of these approximations on 
Eq. (2.9) let us rewrite this equation so that all dependence on v and © is 
explicit. In doing so we will introduce the following notation 


(2.10) (vo) À = y 


| 


Ge l+vm =, 
and we know from kinematics that », = 0.944. Using these we can write 
(2.1) r = 0.944 By? — 1.944 Dy(1 + a,/a))? +1. 


Now let us write out the y dependence of a,. BBCST give the following em- 
pirical and approximate expression for it 


(2.13) Gola, = (— 0.62 + 0.10)y?/r? . 


One can easily convince oneself that it is a good kinematical approximation 
to write 


(2.14) y?/y? = 1 — 0.944y?. 
To approximate (1 + a,/a))? we can use the expression 
(2.15) [1 — 0.62(1 — 0.944y?)]* & 0.65 + 0.33 y? . 


Both of these kinematical approximations are accurate to 2% or so. In ad- 
dition, our conclusions do not depend at all on these approximations which 
are given only to present Eq. (2.9) in a more transparent form. This form is 


(2.16) r = — 0.642 Dy? + 0.944 Oy? — 1.264Dy +1. 


n 
+ 
A 


ds 
‘ 
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T This type of formula: was used by 
BBCST to obtain the curve in their 
Fig. 5. In our Fig. 1 we: show 
Eq. (2.16) for 7 = 1.60. It is marke 
by the solid line. 


Fig. 1. — The energy dependence of the 

negative to positive ratio r near thre- 

Ey Jab. MeV shold. All four curves are plotted using 

a à —— 720 240 60 ro=1.60. The solid curve is the form 

given in reference (1). The dotted line 

includes the correction of a,, the dashed line includes the correction of 0, while the 
[dash-dot line includes both corrections. For quantitative details see the text. 


3. — Possible modifications of the analysis. 


We will consider two ways of modifying Eq. (2.16). These two ways by 
no means exhaust the possibilities, but since our aim is mostly illustrative 
these two examples which, we believe, are particularly realistic, will suffice 
to prove the point we want to make. 


A) The first method is to make an allowance for the approximations 
made in connection with Eq. (2.3). As a consequence of these approximations, 
the first formula in Eq. (23.) is not necessarily 


(3.1) af =(SF R)—a,, 
but rather 
(3.2) ai = (SFR) —a*, 


where the a, is made up of a,, K, a,,, and of the effect of the P ratio. In 
fact, for the sake of simplicity, we will consider the example when 

(338) ast =a = 4, 

We will also assume that the energy dependence of a, is the same as that of a, 
which is quite plausible anyway. In this case the only change we have to 
make is in Eq. (2.15). Let us assume, for instance, that Eq. (2.15) changes to 


(3.4) 0.550 + 0.425 y2 , 
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which means that we change the value of a, by about 25% — quite a con- 
servative assumption in view of our estimates of the approximations involved 
and in view of the experimental error on a, itself, which is more than 15%. 
In this case we have 


(3.5) r = — 0.826 Dy? + 0.944 Dy? — 1.069Dy + 1. 


This formula is also shown in our Fig. 1 for the case of n= 1.60. It is marked 
by the dotted line. 


B) The second way of modifying our formula for r consists of the in- 
clusion of an appropriately small o in Eq. (2.6). In particular, the o we include 
will be, in accordance with BBCST’s Eq. (16), of the form 


PATATE PAESI, 
(3.6) o-a(ffosi, 


where 7= 0,(u/M)?, that is, of the order of (2- 5)%. 

The derivation of the expression for r must be modified somewhat in this 
case. Eqs. (2.1)-(2.5) hold unchanged, and so does Eq. (2.7). On the other 
hand, Eq. (2.8) reads now 


2 0\? 
(3.7) (ran + az)? + (A + aye = (=) (a ee 
Va 
which leads to 
MALE 2 22 PA LATE 
Lo n RIA RUI A ae ect e 
DS) ee Iv va V& | VAY | dI A ¢ ae si 


and finally we get 


1\}2 
0.944y — -|| — 
a) 


2 il a,\* 
— È — —— | 0.944y — -|| [1 + — IL 
2 ‘o 912 By — 7 (0 944 y A ai n) tte 


(3.9) 0 9720y — 2? 
à Peas Y — = 
VM 


where in the last three equations the threshold value of y was denoted by % - 
Using for (1 + &/a) the unmodified form, that is, Eq. (2.15), we obtain the 
effect of the inclusion of o alone. This is shown in Fig. 1 by the dashed line. 
The values of the parameters are 7, = 1.60 and 4/\/a, = 0.05. The combined 
effect of the two independent corrections we have considered is shown in Fig. 1 
by the dash-dot line. 
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4. — Conclusions. 


As we see from Fig. 1, the small and quite conservative corrections we 
made bring about a very substantial changes in the ratio. The reason for 
this sensitivity to small changes lies simply in the analytic form of the ex- 
pression for 7. The value of the function is the algebraic sum of several terms 
most of which are several times larger than the final result. 

It is also evident, however, from the above discussion that good experi- 
mental data will be able to determine the slope of the energy dependence and 
thereby tell not only the threshold value of the ratio but also the value and 
sign of the corrections which we discussed above. In the meantime, however, 
while we lack such data, there is no reason to believe that the threshold value 
. of the ratio is very large and that it is in disagreement with theoretical pre- 
dictions, even if indirect and uncertain indications from the Panofsky ratio 
and the S wave scattering phase shifts seem to point in that direction. 


RIASSUNTO (*) 


Si dimostra che l’estrapolazione usata da BENEVENTANO et al. per ottenere il rap- 
porto negativo a positivo alla soglia è molto sensibile a piccole modificazioni rese pos- 
sibili dalla approssimazione adottata nel procedimento. Tali piccole modificazioni por- 
tano con sè variazioni sostanziali nella dipendenza dall’energia del rapporto vicino alla 
soglia e pertanto l’estrapolazione risulta ambigua in assenza di estesi dati sperimentali 
sul rapporto in questo intervallo di energia. Poichè queste estese serie di dati per il 
momento non esistono, si conclude che non esistono buone ragioni di ritenere che il 
rapporto alla soglia sia in disaccordo con le predizioni teoriche. 


(*) Traduzione a cura della Redazione. 
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Summary. — Diffusion of slow neutrons in an extended uniform medium 
is treated by a method which ascounts for the thexmal motion of the 
medium. If neutron capture is absent, a single integral equation has to 
be solved, in order to determine the velocity distribution and the dif- 
fusion coefficient. A system of such equations is obtained for a capturing 
medium. For the case of a non-capturing monoatomic gaseous medium 
several approximations are discussed, which express the diffusion coef- 
ficient as a function of the atomic weight. 


1. — Introduction. 


Calculations of neutron diffusion coefficients are usually based upon in- 
vestigations of stationary neutron distributions in media which are macroscopi- 
cally isotropic and uniform, and free of neutron sources. For such cases Boltz- 
mann’s transport equation (+2), in the notation of DAVISON (!), reads as follows 


(1) vQ -grad N(r, vQ) + 000) + 27!(0)]N(r, 02) = 


= [for 0, v'Q')f(v'Q! > vQ) dv' dQ’. 


(1) B. Davison: Neutron Transport Theory (Oxford, 1957). 
(2) J. H. Tait: Rep. Prog. Phys., 19, 268 (1956). 
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Herein /,(v) and L,(v) denote the scattering and the capture mean free paths 
for neutrons having the velocity v. The neutron distribution function N de- 
pends upon the position vector r and the velocity vector vQ, with Q being 
the respective unit vector. dQ’ represents the element of solid angle, per- 
taining to the velocity vector v'@' before scattering. 

The scattering function f, which is normalized to unity, 


iù) f(0'L'> 12) dv dQ = 1, 


depends upon v, v and £-£' (the cosine of the scattering angle), and may 
be expressed by the Legendre series 


(2) f(v'Q! > 192) = (21-4 1)f(0 +0) P(Q-Q’) . 


SH 
Ms 


t=0 


The first two coefficients (fo and f,) are of special interest. f is used in the 
theory of the slowing down of neutrons if only their energy distribution has 
to be determined (**"), whereas f, is of importance for the determination of 
the diffusion coefficient. 

If the neutrons are fast, the atoms can be assumed to be free and at rest 
before encounter. In this case, and if scattering is elastic, f involves a 6 func- 
tion (1:2). With slow neutrons, however, the thermal motion of the atoms 
must be taken into account, whereby the function f is more or less smeared 
out (>). 

The mean cosine of the scattering angle and the transport mean free path 
for neutrons of the initial velocity v' are 


(3) b(v') = || [roa VQ) Q-Q' dv ao fi — v) dv, 


(4) lie(v') = 1.(0)/[1— B(v)] . 


For an infinite non-capturing medium the Maxwellian distribution, i.e. 
Nov exp [— fo], with #=Vm/2kT, represents a solution of Eq. (1). The 


3 


(3) G. PLACZEK: Phys. Rev., 69, 423 (1946). 
(4) E. P. WiGNER and J. E. WILKINS, Jr.: U.S. AHC Document 2275 (1948). 
(5) C. C. GROSJIEAN: Mededelingen Nat. Lab. Univ. Gent (1955). 

(5) E. R. COHEN: Proc. Int. Conf. on the Peaceful Uses of Atomic Energy (Geneva, 
1955), Vol. 5, p. 405. 

(6°) H. Hurwitz, Jr., M. S. NELKIN, and G. J. HABETLER: Nuclear Sci. and 
Engr., 1, 280 (1956). 

(7) G. JAFFÉ: Phys. Rev., 88, 603 (1952), equation (56). 
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principle of detailed balance then asserts that (4) 


(5) 1; *(v')f(v'Q! + vQ)v" exp [— Br] = 
= 15 "(v) f(v > v'Q')v exp [— B2v°], 
with relations of the same form holding for f,(v' + v). 


Multiplying both sides of Eq. (1) by dvdQ and integrating we obtain the 
equation of continuity 


(6) divj =—T!{r)n(r). 


n, j and t denote the neutron density, the neutron current, and the mean life 
time of the neutrons, respectively, i.e. 


(7) nr) = | [ f(r, o@)avan, 


(8) j(r) = [flex vR) dvd, 
(9) Tir) — n) [| pente, vQ) do dQ . 


In general 7 depends upon the velocity distribution, and therefore may 
vary from place to place, as indicated. For a Maxwellian neutron distribution 
Eq. (9) leads to the particular value 


(10) To = <l.>/0, 


where 5 is the average neutron velocity (8— 48kT/rm), and <l,> the aver- 
age capture mean free path. The latter average is defined in analogy with 
Maxwell’s mean free path in the kinetic theòry of gases (8°): 


foe) 


(11) le = BA fr exp [— f?0°]dv . 


0 


If 1,(v) is proportional to the neutron velocity, 7 is independent of the 
neutron velocity distribution, viz. 


(12) Gree, Lo). 
(8) J. JEANS: Kinetic Theory of Gases (Cambridge, 1952), $ 108. 
(9) B. Davison: I. c., p. 44. 
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(> is then equal to the capture mean free path for neutrons having the 
velocity %. 

In general the variations of the capture mean free path with the neutron 
velocity will be represented by a function ¢,(v) = <I,>/l,(v). Similarly we in- 
troduce the ratio g,(v) =1,/l,(v), where J, denotes some standard value or 
average of [,(v). The Maxwellian average <I,) defined in analogy with Eq. (11), 
could be substituted for 4; in the special case, discussed in Sect. 3, it will 


however be more convenient to take J, =1,(co) = V(M +1)/M-<l,>. 


2. — Diffusion in a non-capturing medium. 
In order to determine the diffusion coefficient of thermal neutrons in a 


non-capturing medium (/;'= 0), we look for solutions of Eq. (1) which have 
the special form (19:11) 


- (13) N = v* exp [— fv? KO, — C,[a/l, — uUi(v)]}. 
Herein ©; and €, are arbitrary constants; arccos uw is the angle between the 


velocity vector and the positive x-axis. The following equation for U,,(v) is 
obtained from (1): 


(14) 4) | Tale) — file Ud] = 2. 
0 
For a neutron distribution like (13) the proportionality 
(15) j=— D gradn 


holds. The above equations show that 


(16) D = 3x0, 
with 
: (17) We 29: [sto exp [— f?v?] dv. 


0 


(1°) L. Botrzmann: Siteber. Akad. Wien, 88, 835 (1883); Wiss. Abh., III (Leipzig, 
1909), p. 39. 

(1) Compare also 8. CHANDRASEKHAR: Radiative Transfer (Oxford, 1950), p. 14, 
Eq. (80). 
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Hence, in order to determine the diffusion coefficient D, it appears necessary 
to solve Eq. (14). For a specified kernel f, this could be done by numerical 
iteration. 

The following approximate solution of Eq. (14) is obtained if U, (v) is 
substituted for U,,(v') into the integrand 


(18) U,.(0) =fp,(0)[1 — b0)]} = 1, (0)/h . 


When this is introduced into Eqs. (17) and (16) an approximation for D follows 
which can be expressed in the conventional form 


(19) D = ih, 


if the average transport mean free path 7, is defined by 


(20) 1, = 268 Pieters exp [— f?v?] dv. 


3. — Diffusion in a non-capturing monoatomic gas. 


As an example the case of a monoatomic gas, where scattering of slow 
neutrons is isotropic in the center-of-mass system, will be treated in more 
detail. We assume furthermore that the scattering cross-section (o,) is constant 
and that capture is absent. These assumptions make the problem equivalent 
to that with smooth hard spherical molecules in the classical kinetic theory 
of gases (12). It can be expected that the results valid for this idealized case 
approximately apply also to some non-gaseous media (°). 

The only essential parameter of this problem is the ratio M= m,/m of 
the mass of the atoms to the neutron mass; x (Eq. (17)) is a function of this 
parameter. The following auxiliary quantities will be useful: 


0—B(M+1)(2VM), é=B(M—-1)/(2VM), v'Q"=Q—v'Q’, 
E = (60Q + Cv'Q')-Q"= 4 (BV M/o")(v? — 0° + 0°°/M). 
V 2kT/Mm-é is the projection of the velocity which the atom had before col- 


lision, upon the vector 0°", which represents the change of the neutron ve- 
locity. 


(2) S. CHAPMAN and T. G. CowzinG: The Mathematical Theory of Non- Uniform 
Gases (Cambridge, 1939). 
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In the present case simple expressions exist for the functions f(v'2' +02) 
and ¢,(v), namely 


(21) (7) f(v'Q! > oQ) = 620°[4/ Ma: Bo'v'p.(v')F 1e *, 
(22) (#712) ,(v) = [1 + (2M?02)-1] erf (VI fv) + exp [— MB?0?]/(V/aM Bo) , 


with erf (x) = (2//x) fedi. We have chosen h=1(c0)=1/n0,, with m 
0 
being the number of atoms per unit volume, so that y,(co) = 1. 


JAFFE (7) derived the expression (21) by starting with the original form 
of the Boltzmann equation for a binary gas mixture, and by carrying out 
some of the integrations. The same result can also be achieved by geometric 
reasoning, based upon the vector diagram for binary encounters (1°). This 
way has been followed (for M—1) by Hrpy and PAHL (4). Also a quantum 
mechanical deduction is known (8°). 

Fortunately the Legendre coefficients f,(v— vw) of (21) can be expressed 
in closed forms by the error function and elementary functions. We find 


(23) (4) fo(v' > 0) = 20820" p.(v')] + {Lert (00 — Cv’) + erf (00 + Cv')] + 
+ exp [620° — B?v"][erf (00° — Cv) F erf (00 + Cv) ]}, 


(24) io > 0) = 6% M[Bt0"p,(v')]}> - 
 { (B00? — plo — VM) [ert (80 — Co’) + erf (00 + Co')] + 
+ (B00"— Boo? — VM) exp [6? — B*w?]erf (00 — Cv) = erf (00 + Cv)] + 
+ (2/V2)P(0+0) exp [— (Go—Co')?]— (2/2) 8 10 —v' |exp [— (00+%0)?]}, 


where the upper signs refer to v< v', and the lower to v> v. For M=1 
both expressions become simpler and permit a reduction of the corresponding 
integral equations to second order differential equations (49:1915). 

f.(0' > v) can be approximated by simpler functions if v’ is either large 
or small. Equation (14) then shows that U,(0)= 0, 


Lee) 


(25) U (0) = Va + : [Onna — 20Cv?) exp [— 6?v?] dv, 
Zi 
(26) U (co) = f = Sal 


(13) S. CHAPMAN and T. G. CowLiNG: |. c., p. 55. 
(14) J. W. Hisy and M: PAHL: Zeits. f. Phys., 129, 517 (1951), Eq. (18): 
(5) C. L. PEKERIS: Proc. Nat. Acad. Sci USA, 41, 661 (1955). 
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This supplements the picture of U,(v), which is obtained with the approxim- 
ate solution (18). The figure below represents a few cases, and also shows 
the function 2(6v)* exp [— B?v?], by which the function U,,(v) has to be mul- 
tiplied and then integrated, according to Eq. (17), in order to obtain x(M). 

Some special values and approximations of D are known from the kinetic 
theory of gases (12). If M >1, scattering becomes isotropic in the laboratory 
system, and therefore D= $l,v. This means that x(co)=—1. If M=1, D is 
equal to the coefficient of self-diffusion of gases with smooth hard spherical 
molecules (121516), According to PEKERIS (15) the value of this coefficient is 
D = 1.01896:(3/32)71/21,%, so that 


9 
(27) x(1) = 33 V/2-1.01896 = 1.27325 . 


The behaviour of the functions occurring in Eq. (14) permits the conclusion 
that the approximation (19), viz. x=1,-/l,, should be satisfactory for M > 1. 


32(BV) ep 
Zane 


Fig. 1. — Graphs showing U,,(v) for some values of M. The curve for M=1 represents 
the data of PEKERIS (1). 


(Even for M=1 we find the value x= 1.29 which is only by 2% too high.) 
When terms of smaller order than M~ are neglected an asymptotic approxim- 
ation follows: 


(28) x=1+1/6M. 


Some remarks may be added concerning the one-group treatment of the 
diffusion of thermal neutrons. The neutron spectrum is assumed to be Max- 
wellian. If both the thermal motion of the atoms and the anisotropy of the 
scattering (in the laboratory system) are neglected, the already mentioned 


(16) F. B. PIDDUCK: Proc. London Math. Soc., 15, 89 (1916). 
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result x= 1, i.e. D — 140, follows, which is correct for very heavy atoms. No 
improvement is achieved if only the anisotropy of scattering is taken into 
account, while the atoms are still assumed to be at rest before encounter. 
This treatment is equivalent to the substitution of U,,(00) (Eq.(26)) for U(v) 
into (17). The figure explains why the formula so obtained, x=1— b(co), i.e. 


Lv 
1 — b(00)’ 


1 Il 
2 DECISE 
(29) 3lu(00) 5 = 5 


with b(co) = 2/3M, gives the correct limit x(00)= 1, but much too high values 
for small M (e.g. x(1)=3 instead of 1.273). In order to account for the 
thermal motion of the atoms we may substitute for J, b(co) or l,,(co), into (29) 
some reasonable averages of /,(v), b(v) or Z,(v), respectively. Equation (19), 


_ with (20), represents an example. A simpler but less satisfactory way would 


be to keep the denominator in (29), but to use Maxwell’s mean free path 
Cd LV MJ(M +1) instead of {,. The result, though in accord with (28), is. 
again too high for small M. 

The method of Chapman and Enskog consists in approximating U,,(v) by 
an odd polynomial, the coefficients of which are determined by a variational 
principle (17). The result obtained by a cubic polynomial, when specialized 
to the present case, can be written as 


9x) /M +1 13M? + 16M + 30 


20) “= 39) M -12M?+.16M + 30- 


This formula is very good for small M (it gives x(1)=1.271 instead of 
1.273), but not so good for large M (it leads to x(co)=0.96 instead of 
x(co) — 1). This is understandable as the polynomial approximation cannot be 
satisfactory for large M, when U,,(v) approaches a step function (see the figure). 

By a modification, which is insignificant for small M, formula (30) can 
be made to agree asymptotically with (28) and to assume the correct value (27) 
for M=1. A graphical comparison of the various approximations seems to 
support the hope that the modified formula 


(31) 


B) te Fi Mt 4-2.530M + 2.727 
“ M  M°-2.863M + 3.087’ 


is correct to within a few tenths of a percent for the whole range of M. 


(17) S. CHAPMAN and T. G. Cowzin@: 1. c. It should be noted that if the medium 


is a monoatomic gas the above equation (14) is a reduced form of a special case of 


Chapman and Cowling’s equation (5a) on p. 141, with D, corresponding to U,,. 
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4. — Diffusion with capture. 


When capture of neutrons occurs we may look for solutions which are 
proportional to exp [— «/L], where L is an unknown parameter, called the 
diffusion length. Obviously also for neutron distributions of this kind the 
proportionality (15) remains valid. With (6) the well known relation 


(32) dO =a bag 
follows. 

If capture is very weak the diffusion length is large (ZL > 1,). If moreover 
the energy exchange with the atoms is not too weak (i.e. if the atoms are not: 
too heavy) the wanted distribution is almost Maxwellian. More exactly, this 
distribution but little differs from a certain distribution of the type given by 
Eq. (13). The results of Sect. 2 are therefore still applicable and Eq. (16) is 
still valid. In this case t can be substituted for 7 in Eq. (32), so that 


(33) L? = td) 


(or approximately L? = 14,,<4.)). The ratio »=1,/L is therefore related to the 
ratio e=1,/<l.> by v? = 3e/x. 

For stronger capture equations (16) and (33) become invalid. For such 
cases a possible method will be indicated, which aims at the development 


(34) PENA, 
n=1 


where in particular A, = 3/x. According to what has been said before a good 
convergence of this development can be expected if capture is not too strong 
and if all the atoms are not too heavy. (It is easy to see that also the second 
condition is necessary. In a capturing medium with only extremely heavy 
atoms the exchange between different velocity groups of neutrons is so weak, 
that at great distance from the sources by capture an almost monokinetic 
group is selected, which approximately corresponds to some maximum value 
of i,(v). The corresponding diffusion length is not necessarily anything near 
the value 4x1, 1,».) 

In order to determine the coefficients A, etc., we start with an inverse 
series 


(35) ea Be, 


and try to solve the problem by a perturbation method, using the development 


co [hn] 


(36) Nec exp [— 2/L]-v? exp [— Bol DD) ¥" Wayn-2i(0) Ps (y) « 


n= <=) 
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Herein [4] stands for 4n or 3(n — 1), according to whether n is even or odd. 
A system of integral equations for the functions U,,(v) (n= 0, 1, 2,...3 1<n; 
I and n both even or both odd) follows from (1): 


(37) Ps(0) LAC — fre VU (doi 


1 era ICE) 
fa DIRSI] Oger i-i(0) Ste A +3 ces Un-,,141(0) — Gel) 3, BiUnsi,:(0) À 


These equations can be solved in succession. 

The first equation (n= 0, 2= 0) is homogeneous and is solved by Us, = 
= const. We. take Un = 1, whereby the next equation (n= 1, 1=1), beco- 
mes identical with (14). The equations with /= 0, n= 2, 4, 6, ..., require spe- 
cial attention because they are inhomogeneous and have the same kernel f, 
as the mentioned homogeneous equation. These equations therefore have 
non-trivial solutions only if the orthogonality conditions 


foe} 


(38) | È Un-11(0) — pe(v) DI -2;0(0)| 05 exp [— B2v°] do = 0 


0 


are satisfied. These conditions determine the coefficients B;. We find 
B, = x/3, 


(9) B, = - Bf Us1(v) — xp.(v) Uso(v) |v? exp [— B20?] dv, ete. 


Then we have A,= By'=3/x, A, — — B,/B° ete. 
No attempts have as yet been made to apply this method to any par- 
ticular case. 


RIASSUNTO 


La diffusione dei neutroni lenti in un mezzo esteso omogeneo viene trattata con un 
metodo che tiene conto del moto termico del mezzo. Se l’assorbimento dei neutroni 
è nullo, per determinare la distribuzione di velocità dei neutroni ed il coefficiente di 
diffusione deve essere risolta una sola equazione integrale. Un sistema di tali equa- 
zioni è ottenuto per un mezzo assorbente. Per un mezzo gasoso monoatomico non- 
assorbente vengono discusse alcune approssimazioni, che esprimono il coefficiente di 
diffusione in funzione del peso atomico. 
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A. SCHRODER (*) 


Istituto di Fisica dell'Università - Padova 
Istituto Nazionale di Fisica Nucleare - Sezione di Padova 


(ricevuto il 6 Agosto 1957) 


Summary. — The shell model of nuclei is considered, and the average 
potentials for protons and neutrons are calculated on a phenomenological 
basis. Different assumptions about the types of spin-orbit coupling and 
their magnitudes are treated in the one-particle picture. Their relative 
influence on the one-particle levels is discussed for almost the greatest 
part of the periodic system (above 4A 40). Calculations show that a spin- 
orbit interaction — A(%/m_c)?(1/r)(d/dr)V(»)-(Zs/%?) with only one fixed 
parameter À leads to consistent results in one-particle levels: magic 
numbers are reproduced. The diffuse boundary potentials V(r) derived 
for neutrons and protons separately are in agreement with those derivable 
from elastic scattering experiments in the lowest energy limit. 


1. — Introduction. 


The purpose of this paper is to investigate how far the simple suppositions 
on which the shell model for nuclei are based, can yield improved predictions 
as to the energy levels of atomic nuclei available, particularly as regards the 
field of the heavier nuclei (from A x 40 to A & 250), which can be linked 
up with the so-called one-particle structure (in the widest sense of the word). 
Of special interest here is the problem, of whether spin-orbit level spacing must 
be linked in its main part with a overall quality (average potential) of the 
nuclei. To this end the one-particle levels were calculated on the basis of 
several assumptions about the spin-orbit interaction for neutrons. 


(*) A greater part of this work was carried out while the author was at the Uni- 
versity of Zurich. 
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Therefore average central potentials felt by a considered nucleon were set 
up on a phenomenological basis for a sequence of only one ideal isotope for 
a given A, along the domain of stable nuclei in the valley of the total energy 
surface. The nuclear radii parameters found for the potential wells are in 
good agreement with known data. 

It has been demonstrated that a coupling law, which is directly connected 
to the average potential; leads to the right magic numbers and, where applicable, 
to the probably correct size of the energy jumps of the undisturbed one-particle 
ionization energies. 

The same is also obtained for the protons by using an average nuclear 
potential which, for the case of zero neutron excess, has the same shape and 
depth as that for neutrons (Coulomb-potential excluded). Its depth has to be 
monotonically increased with the increasing of the neutron excess. 

The results are listed in a form which is suggested by the sequence of the 
average potentials and their spin-orbit interaction (1). 


2. — The shell model. 


Although the concept of the shell model is well known the basic idea will 
be reviewed. 

The quantum mechanical treatment of an atomic nucleus requires the 
solution of a complicated wave mechanical many-body problem the crux of 
which is the calculation of the wave function (the limits of a non-relativistic 
treatment of the system are hereby reached). The difficulties in managing 
such a problem are made even greater by the fact that the Hamilton operator 
for a nucleus, because of the rather scarcely known nature of the nuclear forces, 
cannot be explicitly set up. Experimentally known data about nuclear forces 
and the nuclei themselves must be used in trying to construct the Hamilton 
operator of the system and therefore it is of very great importance to get 
simple visual interpretations of its constructing elements. 

The predictions which can be obtained from such a Hamilton operator, 
depend for their validity and absence of contradictory elements, on the degree 
to which the lowest approach agrees with the real circumstances. | 

In the shell model this lowest approach is established by the fact that 
the total wave function of the system, like those of a Fermi gas, is built up 
on a total anti-symmetrical product of single eigenfunctions of a Hamilton 


(1) A similar type of diagram was first used by K. BLEULER and CH. TERREAUX: 
Helv. Phys. Acta, 28, 245 (1955). 
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operator, whose potential is a well of finite domain and represents the average 
of all the force effects in the system. 

This procedure makes it possible (*) to define a particle within the correspon- 
ding system with simple quantum numbers as in the Hartree-Fock approach to 
the many particle problem (?). 

It is supposed that this average potential of all the forces acting within 
the system on a particle at the point r represents the average value which 
may be described by a spatial isotropic one-particle force (central field), and 
further acts strongly on the spin. This leads to particle states with the quantum 
numbers nljm for one nucleon. A high degeneracy occurs (because the energy 
of a configuration belongs only to nlj with the half integer quantum number 
j =l-+4 for orbit and spin) specially for energy levels with greater j. Levels 
which are energetically separated from others by large amounts lead to shell 
formation. 

From the radial part of the eigenfunctions associated with a certain }, 
a strong ring formed density distribution is produced (n gives the number 
of ring structures). 

À weak self-consisteney condition has however to be introduced, namely 
that the part of the potential which only depends on the position is not in 
constrast with its total radial density distribution. The potential is fixed in 
advance and guarantees the stability of the system. 

Although the particles in such a field are treated as formally independent, 
their behavior depends on the nature of the origin of their potentials, and 
their properties are closely connected one with the other. In particular, the 
excitation of one particle must considerably alter the energetic conditions of 
all the others and therefore its potential. Such an effect is neglected in the 
simplest approach of the Hartree-Fock theory. Another difference is that 
the individual interactions between the particles are not yet specified (two 
— and many — body forces) so that the calculation of the total binding energy 
of the system will depend on the special convention adopted. A direct fixing 
of the energy zero-point of the system would give rise to a strong self con- 
sistency condition and would base the model on a theory or experimental 
data about the individual interactions. 

The physical existence of an average potential is certainly linked with 
the fact that strong individual interactions between the nucleons may be con- 
sidered as perturbations, so that the wave-function of a nucleus, before and 
after the perturbation has been introduced, may be made up of the one-particle 


(*) This possibility is only approximative because the total linear momentum of 
the whole system must be conserved. The neglect of the proper center-of-mass motion 
results in this case in errors of order 1/A in the binding energies, levels, ete. 

(2) P. 0. L6wDIN: Phys. Rev., 97, 1475 (1955). 
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eigenfunctions of this potential or at least may be made up of products of 
stable core wave functions and one-particle eigenfunctions produced from the 
average potential given by the core. Strong admixtures from other levels are 
thereby not excluded. 

The shell model concept arises from modern experimental knowledge of 
the behaviour of the nuclei at low energy (*). It has, however, considerable 
theoretical difficulties on account of the strong interaction between nucleons, 
which the Hartree-Fock method does not seem to be sufficient to overcome. 
Recently, as a result of the knowledge obtained from the description of highly 
energetic scattering experiments with nucleons, methods have been developed 
which attempt to explain theoretically the correlation and saturation pheno- 
mena in the nuclear matter with two-body forces (4). In this way an attempt 
has been made to solve the problem in the sense of a hypothesis as to the 
very important part played by the Pauli-principle in the matter of the nucleus. 
But by this method the important influence of shell formation is not directly 
taken into account. 

Characteristic alterations in the meaning of the simple average potential 
in the shell model have been made (velocity dependence) by the generalization 
of known methods of the perturbation theory and by further specific develop- 
ment on the basis of allowed individual interactions. 

But till now the simple one-particle concept, which is still free of additional 
hypotheses, has not been sufficiently investigated quantitatively. 

The following shows that the results for the one-particle levels produced 
by a spin-orbit interaction in strict connection with the average potential 
obtained with the simplest model are good, and do not require to be sub- 
stantially modified by invoking a more complicated model. 


3. — The one-particle Schrodinger equation. 


According to the problem we are considering we have to treat the one- 
particle wave functions which are solutions of the equation: 


(1) CHE IVREA 

The most general form of a central force acting on a particle is given by a 
potential 

(2) Vir) = Vie) + fr) (sy. 

(3) M. G. MAYER: Phys. Rev., 75, 1969 (1949); O. HaxEL, J. H. D. JENSEN and 
H. E. Suess: Phys. Rev., 75, 1766 (1949); M. G. MayER and J. H. D. JENSEN: Zle- 


mentary Theory of Nucl. Shell Structure (New York). 
(4) H. A. Berne: Phys. Rev., 103, 1353 (1956) and references. 
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Invariance of (1) against rotation and reflection in the configuration space 
is assumed. Let us take 


(3) “H=T+YV,(r), where DRE 
and 
(4) H'= j(r):(ls). 


According to (2) it is sufficient after separation of the angular part, to repre- 
sent Æ as eigenvalues of the radial equation 


2 9 pae 9 
fa SÈ 041) > M goes 0 do); 


>) dr r dr pa du he 


with the conditions for bound states 


R(0): regular 


(6) so that fr ar: converges , 


U 


R(co): vanishes. 


4. — The empirical data for the determination of V,(r). 


The task now arises of establishing as exactly as possible the potential 
function V,(r). We represent the well V,(r) in the first place by a reasonable 
physical assumption and determine its parameter in the light of the experi- 
mentally well known properties of the stable nuclei. Let V,(r) be the average 
potential that acts on a particle situated at r while all the other particles 
are in their normal positions inside the nucleus, so that their energy levels 
are completely filled with nucleons. The gradient of this function — grad V,(r) 
is the effective force felt by the particle. By this definition, the potential V,(r) 
contains also an energy shift produced by all the other nucleons due to the 
change in radius after the evaporation of the considered particle. A static 
description of V,(r) is an acceptable approximation only if the influence of 
its (virtual) part occurring from the contraction on the eigenfunction of the 
considered particle is small. It seems to be so here. i 

As empirical data for the determination of V,(r) we chose the mean ioniz- 
ation (separation) energy for a nucleon and the suitably defined density radius 
of the nuclei proportional to A*. Both are in close connection with the satur- 
ation phenomenon of nuclear forces. We do not consider a special V,(r) alone, 


19 
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but always a selected sequence V,(7; À) and we try to determine them indi- 
rectly by the above mentioned empirical properties. 

These empirical properties can be represented for a selection of idealized 
nuclei, by a smooth mean value curve provided one does not take into account 
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Fig. 1. — The chosen mean value curve A=N+Z in the N, Z-plane (corresponding 


roughly to the valley on the energy surface for the most stable odd nuclei). —. Separation 
of a neutron pair from a nucleus with even Z and odd N. | . Separation of a proton 
pair from a nucleus with even N and odd Z (contains all cases which can be built 
up from known binding energies). O: Known isobar with greatest total binding energy 
(note the difference between the chosen curve for the most stable odd nuclei). 


their subjective differences. They are simple functions of the particle number A 
for the whole periodical system. Especially in the case of the strongly varying 
separation energy t suitable mean values t(A) must be built up. 

As we want to calculate typical one particle effects later on, only data from 
nuclei with odd A are taken into consideration. For the determination of the 
neutron separation energy we take nuclei with an odd N and an even Z and 
the opposite for the proton. 

We must notice here, that the resulting «even A » nucleus, does not have 
the same one-particle properties (assigments). This comes from the pairing- 
effect in even A nuclei. In order to remain in the connect with odd nuclei 
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we remove à pair of nucleons but form 


— 
Can] 
Sat 


DO > 
Go 5 Gi )- Bee Si z o) for the neutrons , 


7 | 
and respectively 


o al (N NO 
(8) Cala sl Be (5) — Bap E VE 3) for the protons (*), 


di 


NAVE : 
where"B | 4 is the esperimentally known total binding energy of the 


nucleus in question (5) (see Fig. 2). 
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Fig. 2. — The experimentally known values of the separation energy of one nucleon, 
given by the formulae 


I si 


and the determination of the mean values 7,(4) and 7,(A) by means of the chosen 

mean value curve A=N+Z. (The smooth curves 7, from an isotope with N-neutrons 

and 7, from an isotope with Z protons are extrapolated for some cases, in agreement 

with the eliminated effects). From every curve 7 only one point is taken into account 
in order to fix the shell model parameters. 


(*) I am indebted for this important idea to Dr. CH. TERREAUX (see also (1°)). 
(5) A. H. Wapstra: Physica, 21, 367, 385 (1955); I. R. HuIzENGA: Physica, 21, 
410 (1955). 
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In this way through successive ionization of a nucleus one also attains its 
binding energy B,,, without destroying the pure one-particle picture which 
is needed for the lowest approximation of the model. 

From these values 7, and t, we determine the mean value given by smooth 
functions T,(A) and 7,(A4) on the periodic system from a smooth curve for the 
mean values N and Z of N and Z (not necessarily integral) chosen from the 
flow of all stable odd nuclei in the (N, Z)-plane. 

For the total density distribution of the nuclear matter, of the idealized 
nuclei (belonging to N and Z) we make the following assumption in order 
to give a clear picture of the different radius definitions. 


(es fore Oar C2 
C+ t/2) — 
(9) o(r) = | pee SEE 
| 0 » C+Hi2 <r<o, 


which approximately corresponds to the relations derived from electron scat- 
tering experiments. 

If one considers the distribution of the charge only, then (8) must be 
somewhat generalized in view of the electrostatic potential that is to be cal- 
culated from it. (Experimentally such a distribution is only tested in the 
case of 1°”? Au). 


4 2 
0° |! tan ea | for 0<r<O'—1#//2, 


(9a) A oll + 6 PE — i Mt DOT 
0 » C0+t[2<r<0c0. 
By means of the moments 
(10) Je ter] etre àr With nm 02 vue, 
the «half density radius » | 
(11) c= | mario, 
J 


10 = 470! DL C3 + te cul 


465 


ONE-PARTICLE LEVELS WITH SPIN-ORBIT INTERACTION ETC. 469 


the «root mean square radius » 


2 _ via  Q(8/5+0/2) (6/0) + (8/80)(/0)\! 
(12) Teo = VI — 0 | i sa a/Hwoyr AAO | 


sa LE Ai 3 
MO SON 


5\4 
(13) y= (3) les , 


(14) To = RrA * , 
(15) r = CA, 
(16) PR AS. 


From (11) and (14) it follows that 
(17) Po MEC) 


Only in the case of a proportionality between ¢ and C, if 7, were constant, 
would 7, also be constant. Experimental evidence contradicts this (5). 


5. — The neutron potential. 


Now a potential V,(r) is built up which should be in agreement with the 
distribution of the density (9) by taking into account a finite range for the 
nuclear forces. (But no assumption about the forces has been made). 


(5) B. HAHN, D. G. RAVENHALL and R. HOrSTADTER: Phys. Rev., 101, 1131 (1956); 
R. HorsTADTER: Rev. Mod. Phys., 28, 214 (1956). 
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for O0<r<kR—a, 


R_-a<r<k, 


RÉ CD - 


(One could certainly give a more analytical description, however this more 
schematical form is most suitable for the applied calculation method without 


- Density distribution and potential. 
(Definition of the parameters). 


any particular auxiliary calcula- 
tory devices). «a» determines a 
diffuseness of the potential boun- 
dary. Let À be defined as fol- 
lows (see Fig. 3). 

(19) R=rA?+b=C4+d. 
The parameters 7,, b, a, V, must 
now be chosen so that the above 
average separation energies T,(A) 
can be reproduced by the solut- 
ions Æ of (5) with f(r) — 0 and com- 
pletely filled energy levels (with 


N(A) particles). This determination of the parameters is a rather sensitive one. 
By the application of the procedures described in Appendices I and II, 
this problem can be solved by means of a lengthy and extensive calculation. 
We attempt to take r,, b, a and V, to be constant. 
Results: the given values T,(A), N(A) and the conditions 


Py = COMM - 
D” = conste, 
V,= const, 


Q =const, 


where chosen. It therefore follows that 


r, =1.18-10-% em, 


b = 2.05-10-% em, 


V,— 44.5 MeV, 


= Vea ear 
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6. — The spin-orbit interaction. 


A satisfactory theoretical explanation of the high doublet splitting of the 
one particle energy levels which would lead to a specific shell structure for 
nuclei, and to the formation of the so-called magic particle numbers, has not 
yet been given. It is therefore important to establish if the operator H' con- 
tained in (1) in the form (4), like V,(r), possesses a significance which covers 
the periodic system of the nuclei. i.e. whether f(r) is connected with an overall 
property of the nucleus in question (the nuclear field), or if H’ undergoes 
notable modifications as a function of the particle number. (In the latter 
case H' could in a certain way, represent rather a caricature of an effect 
produced by an individual two body tensor force interaction between the nucleons 
in the shells in question similar, for instance, to the behavior of the eliminated 
pairing energies). This second hypothesis would however make the simplest 
approach to the shell model very doubtful and would make senseless the level 
representation given later as continuous functions of A. 

In order to test this, we again consider the case of neutrons in the nucleus. 
Making the assumption suggested by the determination of V,(r), namely, 
that it is not considerably modified by the addition of H’ to the one particle 
operator H (average separation energies are approximately equal with H’ 
added), the adaptation of V,(r) by means of t, would otherwise be useless. 
Next, the angular momentum j resulting from the spin-orbit coupling of a 
particle is to be diagonalized with respect to H+H' (which leads to the 
jj-coupling scheme). 

For the case 1 # 0 this will be reached by means of the well known complete 
reduced product representation of the rotating group 

D, xD, =D; ,,, + D 


j+t—}° 


The matrix of Clebsch-Gordan coefficients, which performs the reduction is 


1 Wise me Nm 
(20) (G) = =" I a 3 
VA+I \VI— m =—VI+m4+1 
For both directions of the angular momenta it follows that 
x He |! for j=1+4}, 
9 DE + 
Di Hr) 5 (FSET) » j=1—+ with #0, 


and the eigenfunctions of H are then: 


55 Cond ae Fe =) eal Y7'(9, @) | 
(22) m +1/2 (r, j °,œ) Co 0 Ris = 12 (9 ) va (9, 9) . 


ny=l1—- /2 n,v 
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From (21) it follows that the radial eigenfunctions Rr) for larger { depend 


closely on j. 
For the sake of simplicity the assumption is frequently made that 


(23) fir) = —|const |, 


where the minus sign is required in order to split the level with j = l + 1 down, 
according to the empirical facts. The centre point of the occupied doublet lies 
in this case exactly on the level formed without splitting, so that the deter- 
mination of V,(r) is not in any way influenced. 
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Fig. 4. — Energy levels corresponding to the shell model according to two different 
assumptions about the spin-orbit interaction for a single neutron. (Potential wells 


determined by neglecting spin-orbit splitting). 
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In the lower part of Fig. 4 this coupling is represented with the special value 
(24) IERI 


(corresponding to an approximate adaptation to the level splitting of expe- 
rimental values for A = 40 up to A = 50). 

We can see here that the magic numbers cannot be employed satisfactorily 
for the héavier nuclei: the splitting of levels for greater A is much too large. 
The assumptions (23) and (24) are not realistic. Closely linked and employed 
by many authors is an assumption which binds f(r) with V;(r) in a simple 
way and also satisfies the invariance condition for spatial inversion: 


VISCO: 
ZI F == I À | — — . 
(25) f(r) =—|const| = = Vil") 


This corresponds formally to the well known Thomas-term for the magnetic 
« doublet splitting » in the structure of atomic alcali spectra. H' is here pro- 
portional to 


~ s(grad V,(r) xp) . 


It must be noted that such a coupling term, because of the magnitude of the 
coupling energy in question, can no longer be taken into consideration by means 
of the application in the usual manner of an ordinary perturbation calculation. 
Therefore, for every coupling parameter involved, a new calculation of the 
whole level scheme for 1 4 0 is necessary. The conditions are such that the 
corresponding method, described in the Appendix: II, can suitably be used. 

Results: the result of the calculation is given for three different values 
of the dimensionless coupling parameter À in Fig. 5 (also for the case with 
good adaptation in Fig. 4 above). With 


2 ‘ for = 1 
mc} rdr AO) END 2 es 1d witha 0) 


A= 0465 ADR. A= 0.29 , 


and applying V,(r) from Sect. 5 (*). 


(*) The Compton wave length of the x-meson is introduced for dimensional reasons. 
With the Compton wave length of the nucleon the best adaptation gives 2 = 18.4. 


ie 
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The best agreement with the average of the experimental data can be seen 
for À — 0.38. Fig. 5 shows that by comparison over the domain of À already 
À can be relatively exactly determined only from known experimental level 
assignments and spacing in a few cases. The applied coupling law represents 
the experimental facts relative to the one particle attribute of nuclei with 
shell structure. Only in the field of the nuclei from A & 160 upwards, (the 


T4 
(MeV) 


6 


0,46 


À 


0,38 


À 


20 


à 


Yo | 
oe 
2 \ | 
) 
LE 
ae, 
d? 
ÿ 
© 


0, 20 40 60 80 100 120 140 160 180 200 220 240 


V,, = 44.5 MeV = const R = (1.184% + 2.05)- 10-13 em G= 2.83 Oem 


a 


Fig. 5. — Energy levels corresponding to the shell model with the spin-orbit coupling 


— A(h]m_c)?(1/r)(d/dr) V(r)-(Es)/#2 and three different parameters 2 (potential ‘wells 
determined without taking into account spin-orbit splitting). 
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lightest nuclei have not been taken into account because the structure of their 
core seems to be more complicated) is the agreement less good. But in the 
calculation made, it has not been taken into consideration that in this domain 
the more rapidly increasing coupling energy between spin and orbit alters 
the value of the mean separation energy (Fermi energy) somewhat to the 
value of 7,(4) respectively. 

Therefore a new determination for V,(r) is necessary. Further, a rounding 
off of the potential well certainly brings a notable raising of the level having 
a high orbital momentum. 

Taking these arguments into account a new adaptation of the potentials 
V(r) (now V,(r)) was set up with the direct employment of H’ with the same 
grad V;(r) on the surface and 2 = 0.38 where the «corners » of the potential 
boundary in the region of 1.2 10-!* cm were rounded off. 

Results: the given values now are 7,(A), N(A), grad V,(7) on the surface 
and the conditions b°= const, r}= const, the final parameter values for the 
neutron potential. 

For A > 60 are: 


po OS One 
DEA 00M 
VIa 46 MeV, 

EN = IAS ia 


These results are given in graphical form in Fig. 6, and in Table I, also for 
values below A = 60. 


7. The proton potential. 


If the given spin-orbit interaction has a universal significance, tais fact 
ought also to lead to reasonable results in the case of the protons: this will 
be demonstrated in the following. If in the first place we consider the po- 
tential V,(r) to which, according to its definition there is a small contribution 
of electric energy formed by the closing up of the charges in the nucleus when 
a nucleon is eliminated (this contribution appears also in the neutron potential), 
then this would have to be completed with the Coulomb potential of the proton 
charge distribution. If in order to do this by means of the electrostatic po- 
tential (the influence of the electric exchange effect ought not to be big, as 
a statistic estimate shows) it follows that the condition for which the limiting 
energies of the Z(A) protons must agree in the mean with 7,(A) can be ob- 
tained with V% only in the field of the nuclei with N(A) » ZA). In nuclei 
with an increasing excess of neutrons the protons would however find too 
little space within the nucleus. 
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Fig. 6. — Pure one-particle separation-energies for the most stable nuclei derived from 

the shell model. The energy jumps corresponding to neutrons and protons for equal A 

are of the same order of magnitude. A 68-stability condition between neutron and 

proton levels brings irregularities in the total energy vailey just in those places where 

the differences between corresponding neutron and proton levels are great (but then 

we must take into account pairing effects). The calculated potential parameters 
are drawn below. 
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The parameter V° — V° has therefore been treated as a variable in order 
to take this into account. For E° and @ the values obtained for the neutron 
potential were maintained. The Coulomb potential 


is to be added to V,(r). For spherical symmetry o,(r) from (9) or (9a) gives, as 
is known 


Tr 


@,(r)/e = 47 (feta da + «(ar dr) i 


0 


(27) 


within the charge distribution 


d> 0 


Dro = dn? | oxler)a? de , 


0 


outside the charge distribution and must be normalized by 


Z—1= ta oa) din 


0 


(28) 


For the charge distribution (9a) it follows: 


[ Dir) + O1(0'— 11/2) for 0<r< C4 I 
V(r) =e DC Or) DR CID e CRE RIT 
| Dir) DA) ioe Cat epee, 


where I and II refer to the given 
integration limits. Because of the 
simple form of (9a), only elementary 
integrals result. In these d= (1 — q)/q 
(0° = 0.) With the linear function 
LL JU) g (Aye 9 (250) 
= 0.72, OC’ =1.07A?- 10-13 em and 
v= 2.4-10-13 em were used to take 
into account at least approximat- 
ely the central depression of the 


5 10" Cb/em?) 


= 
(©) 


Fig. 7. — Charge distribution and potential. 
(Definition of the parameters). 


charge density affecting the Coulomb 
potential. 
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The depression within the nucleus has a distinct influence on the situation 
of the proton levels in heavier nuclei, as the calculations have shown. V.(r) can 
very well be substituted by a Coulomb potential created from a uniform density 
distribution with a slightly variable r,, between 1.11 and 1.13 (10-13 em). 

We can therefore simply take 


(270) A SA [(rA)/r for r > To A 
A | (375.4% — r2)/2 » T< A}. 
The parameter V) was then chosen so that the mean separation energy for 
the protons (taking into account the spin-orbit interaction (26) with V,(r) = 
= V,(r) +V.(r) and A — 0.38) was satisfied. 

While for the calculation of the levels in case | 0, the calculation pro- 
cedure, according to Appendix II, was very suitable, an ordinary perturbation 
calculation for the s-levels was carried out in order to take into account the 
« wine bottle » effect of the potential well. 

Results: the parameter values for the proton potential are given in Table I, 
and a graphical representation of the result in Fig. 6. It can be seen here 
that the same spin-orbit coupling in the proton case gives reasonable results 
when the potential V,(r) + V,(r) is set up in the given way. 


8. — Comparison with experimental facts. 


In order to obtain a direct comparison of nuclear spin and parity (measured 
from the ground and lowest excited energy levels, which might be interpreted 
as one particle levels in the widest sense of the word) with the calculated 
level schemes, these, have been drawn separately for neutrons and protons 
together with all measured assignments (assignments have been taken from 
the « Nuclear Data » (*)). Care was taken to plot only such excited levels which 
might be considered as certain. Doubtful assignments were excluded (see Fig. 3 
and Fig. 9). 

A glance at the measured values shows clearly the appearance of magic 
numbers, especially when spin and parity change (7). On the basis of the cal- 
culated pure one-particle ionization energies, we are here able to distinguish 
between magic and half-magic particle numbers. The magic numbers appear 


(*) I am indebted to Dr. CH. TERREAUX who put a compilation at my disposal. 
(7) This type of assignments diagram was first used by R. LATHAM: Proc. Phys. Soc., 
59, 979 (1947). 


a 


~ 


480 A.” SCHRODER 


in cases when the energetic change in the calculated ionization energy after 
the filling of a sub-shell is greater than about 2.5 (MeV) in the model. 
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Fig. 8. — Shell model level scheme for neutrons together with the measured values for 
spin and parity from odd A nuclei: Bi ground states; O excited states. 
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Fig. 9. — Shell model level scheme for protons together with the measured values for 
spin and parity from odd A nuclei: m ground states; D excited states. 
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Results: the following magic numbers result: 


Neutron numbers (theoretical) 


| M 20 28 40 50 82 126 
hM || 14 16 32038 56 64 90 100 114 118 124 136 


Neutron numbers (from known experimental data for spin and parity) 


M | 14 16 20 28 32 38 40 50 56 64 82 (90) 100 (114) (118) 124 126 


Proton numbers (theoretical) 


M 20 28 50 82 
hM || 14 16 38 40 58 64 76 80 92 


Proton numbers (from known experimental data for spin and parity) 


| M | 14 16 20 28 38 40 50 58 (64) 76 80 82 


The numbers in brackets () are those for which the evidence from Fig. 8 
or Fig. 9 for the appearance of a magic number is very small and uncertain. 

An extensive comparison between the measured values for spin and parity 
of the observed odd nuclei shows that nearly all the low lying energy levels 
may be assigned as one particle states of the given central potentials V(r). 
This is also valid for the highly degenerate half-filled sub-shells (the known 
break-down of the normal structure in the shells 7/2- and 9/2+ with [=j —1 
for neutrons and protons separately by their first occurring in the periodic 
system being excluded). 

From the point of view of a perturbation theory, in which individual inter- 
actions occur (deviations from the average central potential) not only the 
degeneration is annulled but the perturbations will lead to strong admixtures 
to the unperturbated wave function from energy states in the close neigh- 
bourhood of the last occupied state (from neutrons and protons) in the one- 
particle picture. 

Excited hole states occur after a sub-shell is filled, and arise from the 
displacement of a particle from the occupied state to the first unoccupied state. 
These holes therefore lie energetically below the separation energy. 

Only in a few cases where admixtures are expected to be small one can 
test the calculated levels directly (closed shells plus or minus one nucleon) (8). 


(8) Levels in the neighbourhood of 1 f 7/2 and Ca isotopes: R. H. NussBAUM: Rev. 
Mod. Phys., 28, 423 (1956); C. M. Braams: Phys. Rev., 108, 1310 (1956). — Levels 
of 2°7Pb: I. A. Harvey: Can. Journ. Phys., 31, 278 (1953); D. E. ALBURGER and 
A. W. Sunyar: Phys. Rev., 99, 695 (1955). — Levels of 2° Bi: M. J. PooLE: Phil. Mag., 
44, 1298 (1953); E. A. Error, D. Hicks, L. E. BEGHIAN and H. HALBAN: Phys. Rev., 
94, 144 (1954); R. B. Day: Phys. Rev., 102, 767 (1956). 
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TABLE II. 
#Ca (particle) Ca (particle) 207Pb (hole) 
20 21 20 29 82 125 
exp. theor. exp. theor. exp. theor. 
TEE 0 £ 7/2" -0 EROI, p 3/2 0 12200 p w/2-.0 
yes 1:95% pie) 2a 3-8 (1/2) 2.02): p.1/2; 02:02 5/25. 0569-1 f 5/2- 0.8 
LES TE D Loe 6 Po 9 )2 008940 DS /2 NE 
5204 op 6:7 13/2* 1,633 | i 13/2419 
10) 22:33 SA ete 28 O20 
4$c (particle) 4%$c (particle) 209Bi (particle) 
21 20 21 28 83 126 
apo Sette 712770 ELA PEA) 0 0072280 
SORT 0 ED S/o 14.1 2 eed LU np 3120 4.2012 (7/2.).0.9 10712 el-.4 
(5/2) >4 f 5/2 6.4 SCO eg Cd 3/2 RS 
2.6 DIRO ET SS 
| 
Measured and calculated levels for spherical nuclei with closed shells plus or minius one 
nucleon. 

(Note: for the given Ca and Sc isotopes the discrepancies between the calculated levels are 
great, but we are near the top of the wells and in this domain small variations of the shape 
have a great influence on the levels because centrifugal — and Coulomb — forces are not high. 
For ?Pb and *°°Bi this situation is completely reversed). 


In particular for those parts of the periodic system where marked deform- 
ations on the nuclei are observed, the given shell model spin-orbit interaction 
together with the central field potentials V(r) predicts precisely a high dege- 
neration in the unperturbed one-particle levels and a very large level density. 
If this is the case for neutrons together (N = 90--114; Z = 60--75), one 
can not longer expect that central field shell-structure in the perturbed case 
is a useful approximation (strong perturbations between the shells and the- 
refore polarization of the core). Here the simplifying Bohr-Mottelson collective 
particle movement description leads to good results (°) but even the shell 
model situation suggests this. 

A detailed analysis of the experimentally known shell structure for odd 
neutron nuclei related in a new way to the shell model is given recently by 
K. BLEULER and CH. TERREAUX (1°). 


(9) A. BonR: Dan. Mat. Fys. Medd., 26, No. 14 (1952); A. BoHR and B. R. Mor- 
TELSON: Dan. Mat. Fys. Medd., 27, No. 16 (1953); K. GorTFRIED: Phys. Rev., 103, 
1087 (1956); S. A. MoszrowSsKI: Phys. Rev., 103, 1328 (1956). 

(1°) K. BreuLER and CH. TERREAUX: Helv. Phys. Acta, 30, 183 (1957). 
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À comparison is given there with the level scheme calculated here for 
neutrons. 

Particularly for the special cases of light nuclei with a stable core and 
an extra nucleon the spin-orbit coupling determined by the author, somewhat 
prior to the preparation of their paper, leads to reliable results. 
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Fig. 10. — Measured and calculated levels for spherical nuclei with closed shells plus 


du 


or minus one nucleon. —: excited state; #: ground state; #: hole state. 


The nuclear radii obtained from the adaptation procedure of the potential 
wells are consistent with such produced by a total density distribution which 
is expected from the Stanford measurements of the nuclear charge distribution. 
Between the radii C of the half density and R—a/2 of the half potential a 
difference d = 0.93-10-! cm is obtained. 

No difficulty arises in taking account of the Coulomb potentials V(r) in 
the nuclei (V,(r), V,(r) and V,(r) are potentials, not potential energies per 
particle). 

A difference between the average neutron and proton nuclear potentials 
V,(r) and V,(r) obtained in the case N Z is closely connected with the anti- 


bui 
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= 


hehe 
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Symmetry of the total wave function and can give information about the 
character of the forces (ordinary and exchange) between nucleons in nuclei. 

The potential wells correspond in a good manner to the real part of the 
optical-model potential (1) for the limiting cases in which the energy of the 
incoming nucleons tends towards zero. This can be tested for example with 
the positions of the s- (size) resonances for neutron scattering. 


Potential 


Density 


NUCLEAR RADIUS (10cm) > 


20 40 60 80 100 120 140 160 180 200 220 240 A 


Fig. 11. — Nuclear radii consistent with the mean separation energies 7,(A) and 7, (4). 

(Potential wells determined by taking into account spin-orbit splitting). O: Stanford 

results from high energy electron scattering experiments (with the assumption of a 
Fermi charge density shape). 


For the potentials V,(r) the calculation gives 


38: at Aw 60, 


As: at Awl, 
and the experiments. 
38: at ASSE Di 


48: at A » 160, 


(11) H. FesHBACcH, C. Porter and V. WEISSKoPF: Phys. Rev., 96, 448 (1954). 
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(but it should be born in mind, that the determination of the s-levels near 
the top of the wells gives values for À which are somewhat too high). The 
spin-orbit coupling parameter de- 
VARA termined by the analysis of proton 
È IRR polarization data (12), possesses, 
uy a te within the limits of error, the same 
value as that reached here by 
other methods. 

It must here be remembered 
la°:2. 43:10 "cm const that when comparing this paper 
with others on the same subject, 
F7 ys Example: the nucleon mass M is not reduced 

4 Ke inside the nucleus; nevertheless 
it represents the experimentally 


Fig. 12. — Mean potentials calculated for the i. 


heavy nucleus 2°°Pb. (Note the Johnson-Teller integer spin satisfactorily (*). 
effect, caused by Coulomb potential on protons). The proof of self consisteney 

between the average potentials 
and the zero-order radial density distributions, will be given in a later paper, 
together with a study about the nuclear core. 
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I am very grateful to Prof. Dr. K. BLEULER for having encouraged me 
to undertake this research, also to Dr. CH. TERREAUX for many encouraging 
discussions, and to Prof. Dr. W. HEITLER for the helpful interest he has taken. 
I would like to thank the « Martha Selve-Gerdtzen-Stiftung » for the grant 
of a scholarship which has made it possible for me to carry on the greater part 
of the calculations. I am also much indebted to Professors A. ROSTAGNI, 
N. DALLAPORTA and A. KIND for their kind interest and the INFN for a sup- 
port to finish the work. 


(12) T. Eriksson: Nuovo Cimento, 2, 907 (1955); S. K6HLER: Nuovo Cimento, 2, 
911 (1955) (with different energy), see also L. VAN Hove: Physica, 22, 979 (1956). 

(*) A velocity dependence of the average potential V(r, p) represented in a one, 
particle picture starting approximation, is the result of an enlargement of the shell model 
discussed in Sect. 2. Here we have to symmetrize V(r,p) in order to obtain a her- 
mitian operator, but the results will depend on the method used, since there is not a 
unique solution. A result of the present work is, that we cannot prove the effect of 
a mass reduction in the now formal kinetic energy operator 7(r) [which is the first 
approximation of the velocity dependence of V(r)] on the low lying one-particle levels. 
The energetic differences between the normal occupied shells inside the potential wells 
would give better information about the effective mass problem, but there is so far 
no direct experimental information in nuclear spectroscopy. 
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APPENDIX I 


The adaptation of the potential well (schematically). 


For this purpose it is best to transform (22) by rR,,,(r) = y,,(r) with the 
new limiting conditions 7(0) = y(oo) = 0 to form a one-dimensional Schrô- 
dinger equation and then use the approximation procedure for the solution 
which is given in Appendix II. 

In order to adapt the parameters we make use of the properties of the 
Schrôdinger equation for which the relative eigenvalues e are invariable with 
espe ct to the following transformation: 


(29) De, 

(30) D no) with E=5: 
Irt follows that 

(31) re + ee) CN) ye) —o, 
with 

(32) = pe IRE 


and by means of (18), where in particular: 


| — 1 for OE 

(33) HD —- 1-8 2A ilo ea 
| 0 » 1]< & < col, 

with 

(34) y = a/R. 


The parameters responsible for the essential physical properties of the system 
appear together in X, whereas the behaviour at the potential limits is repre- 
sented by y. 

Potentials with the same X and y are therefore equivalent to their relative 
eigenvalues e and one can in this way, by calculating (31) for a suitable domain 
of X and y, adapt the parameter so that the conditions necessary for the mean 
value of the last fully occupied level (according to the Pauli principle) are 


= 
on 
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satisfied. However, for this purpose a three-dimensional representation e(X, y} 
is necessary, the calculation of which demands an excessive effort. An approx- 
imation procedure was therefore resorted to. 


US U 


Fig. 13. — Approximation method for the determination of the potential parameters. 


We next substitute all the parameters r, b, a°, V° with constant upper limits. 
r, b, a, V: they can be reached in a relatively easy way by means of expe- 


rimental data. We then have: 
[ RE =RLAR, 
(35) | PSP EA Ves 
| a =a +e, 


and the level scheme is calculated over a domain of A. If by means of the 
present relations one now sets up the suitably fulfilled condition that if 


LOL 
(36) pia ae 


Il 
wins const, then E(V, a) & EV°, «)— AVS, 


and consequently 


(37) R=(R°—a)+a=U+KkV, 
(38) e =(Y+T(U)/Y. 

In particular however 

(39) X=pVVR=pV(U+KY), 
(40) (ee we 


Ruy Dh Uk VE 


and this leads to the solution 
1 
(41) gu = (7-1), 


(42) V = (Xyluk} . 
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Whereas for the adapted values 


(43) x = po(VRR, 
al 
(44) y? == Re , 


‘ought to be valid. Inserting (43) and (44) in (41) and (42) (for X°—X, y°=y) 
leads to ; 


wae Ne Reo ae TP, 
(45 Ue — ARE ET o = Ul 3 — 
(45) U ae m) ( * ELA, 
0 40 __\# 
(46) A a 
La 


and according to (38) it follows for one energy level that 


fa 
(47) Ba (P— TU). 


When one sets up this levels in the right way with N = f(A) particles for 
the nuclei in the trough of the energy surface, its mean limiting curve may 
be made to agree with t,(A) for a certain parameter combination over the 
whole periodical system. Then the level scheme can be recalculated, and so on. 


APPENDIX TI à | 


‘The approximation method for the solution of the Schrédinger equation for 
a potential well. 


The solution for the Schrédinger equation 


hh? 
— ap 0) + (Mil) — mM = 0, 
with 
+ h? (+1) 
HONOR: 


are required for the potential (18) (*). 
For the s-states with /=0 it can be solved with analytical eigenfunctions. 
In the case of 1 40 we replace the potential which appears in the one dimen- 
sional Schrôdinger equation as W,(r) with a suitable approach potential, for 
(*) The spin-orbit potential is excluded here for simplicity, because the method 
for the solution remains the same. Only the parameters 4, v, w are then functions of n, I, 7. 


489 


490 A. SCHRODER 


which there are analytical, and therefore easily manageable eigenfunctions. 
In the case of the potentials used in this work, the following are particularly 


suitable : 
= 0h) D 
il = Wii) 7 
Wi 
Morse potential: 


2 exp ( ma 2) — exp | 2 ca = WE (r) . 


The parameters u;(n, 1), v(n, 1), w,(n,t) with à — 1 or 2 must be chosen so 
that the approximate eigenfunctions obtained with them, have the same 
turning point on the potential boundary as the exact eigenfunctions. They 
must satify the condition 


oscillator potential: 


— Us, 


Jon Lee Wit) vee) dr = 0 ’ 


with (W,— W,,): small in the domain of 7#(r). This adaptation was carried 
out graphically (*). The results for the approximate eigenvalues are then 


h 1/2u, il 
ioe == | 1 
n Ui + à M n 9)? 


‘Quy h /2u 1 
E® — 2 I. | 3 > 
Lari, SV (n- 3) uso, Ÿ M\" 2)? 


and the exact eigenfunction y,,() is analytically approximated by yf(r). 
Only the values n = 1, 2,3 are necessary for the radial quantum number. 


(*) See for a special case without spin-orbit interaction: K. BLEULER and CH. TER- 
REAUX: Helv. Phys. Acta, 28, 260 (1955). 


RIASSUNTO 


Si calcola fenomenologicamente, considerando il modello a shell, il potenziale medio: 
dei protoni e dei neutroni. Si discutono diversi tipi di accoppiamento spin-orbita e i 
relativi ordini di grandezza; si analizza la loro influenza sui livelli ad una particella, 
per la maggior parte del sistema periodico (da A & 40). I calcoli dimostrano che una 
interazione spin-orbita del tipo — A(h/m,c)?(1/r)(d/dr)V(r)-(ls/h?) è un parametro fisso, 
porta a risultati consistenti per la determinazione dei livelli ad una particella: sono 
riprodotti i numeri magici. Si ottengono inoltre dei potenziali diffusi per protoni e 
neutroni, in accordo con quelli derivabili dalla diffusione elastica elettrone-nucleo a 
bassa energia. 
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On the Creation Operators of Physical Particles. 


G. MOHAN 


Department of Physics, University of Nebraska - Lincoln, Nebraska 


(ricevuto il 28 Agosto 1957) 


Summary. — The conventional annihilation and creation operators of 
free particles are modified in view of the fact that all measurements re- 
quire non-vanishing duration of time. The modified operators create and 
annihilate physical particles. Asymptotic assumption or adiabatic hypo- 
thesis are not invoked so that the renormalized S-matrix can be obtained 
without their help. The analysis of certain functions that occur in the 
definitions of the modified operators is incomplete. 


1. — Introduction. 


The renormalization procedure in the theory of the S-matrix is a process 
of converting the transition probability amplitudes between non-physical bare- 
particle states to those between physical particle states. However, it is de- 
sirable to develop a theory which avoids introducing non-physical states even 
in the intermediate stages. It is possible to write down neat expressions for 
renormalized S-matrix elements directly, but their derivation is always based 
on either some sort of asymptotic assumption (21) or on a switching technique (*) 
which has not been completely established and which looks quite non-physical. 
Our approach to the problem is entirely different. We attempt to improve 
the technique of construction of the in- and out-states. The conventional 
annihilation and creation operators are modified on the basis of the fact that. 
no measurements are instantaneous. The modified operators refer to an ex- 


(1) H. LEHMANN, K. Symanzix and W. ZIMMERMANN: Nuovo Cimento, 1, 205 (1955). 

(2). K. NISHIJIMA: Prog. Theor. Phys., 17, 765 (1957). Other references are listed 
in this paper. 

(3) G. KALLÉN: Physica, 19, 850 (1953). 
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tended time which should correspond to the duration of measurements (*). It 
is shown that these operators are the «renormalized » creation annihilation 
operators. Starting from them one can directly obtain the renormalized S-matrix 
by the usual methods. 


2. — Instantaneous creation and annihilation operators. 


We shall confine our attention on a system with pseudoscalar interaction 
between fermion field Y(xv) and boson field (x). The method developed is 
applicable to other renormalizable theories also. The field equations and 
commutation relations are 


0? es 
(ea — 1) D(x) = dx? D(a) + g P(x}y: Pa) + A D3(æ) = — J(x), 
(1) È 


Gr + n) Va) = dm) — gys LD) = I2), 


[P(0), P(0)]p lose = ty? dle — x’) , 
OD(x') 


dx! 0. 


(2) 


P(x), 


=10(x—x’'). 


a =x"0 


Corresponding to every instant 7 a set of annihilation and creation operators 
for the bosons and Fermi particles and antiparticles can be constructed from 
the Heisenberg fields. For example, the creation operators are (*) 


ci) i fara & f(a) Pa) ; 
(3) az (pe) = — i [ata opt; æ) , 


vt 


pati | dn bap; ap Pa), 
where 


fa(æ) = 5A exp[ika];  kR=Vk+x, 


ug(p; 2) = (27) * VE up(p) exp [pe em pe V p? +m; cete. 


(+) (E. C.G. STUECKELBERG has studied S-matrix with non-precise epoche of initial 
and final states separated by finite time interval; Phys. Rev., 81, 130 (1951). 


0 of dg 
(5) È oi g(x) = 5,09 130” 
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These satisfy the typical commutation relations of creation and annihilation 
operators, for example, 


[a,(p'8'), a (pf)]}- = sy (p — p') 5 
La,(p'B’), a, (pB)]+ = (|). 


(4) 


It cannot be asserted in general that c*, a* and b* create physical part- 
icles when acting on the physical vacuum state |0). In order that the state 
vectors c*|0), a*|0> and b*|0> represent physical particles they must be 
eigenstates of — P,P" with eigen-values x?, m? and m? respectively, where P“ 
is the total energy-momentum operator and x and m are the masses of a 
physical boson and a physical fermion. For finite 7 this condition is not 
satisfied, but if appropriate limit to remote future or to remote past can be 
taken, then, indeed, the above mentioned states represent physical part- 
icles (4) (the so-called outgoing and ingoing particles). Such a limiting process 
will be proposed in this paper. 

In order to construct a complete SONORA set of multiple physical 
particles by the repeated application of the creation operators it is obviously 
not enough to show that the proper limits of c*|0», a*|0> and b*|0> are 
physical one-particle states. One also needs to establish commutation relations 
of the following type 


lim 2, "Le, (k), e (k)]_ = Ô(k — k'), 


(5) 
lim 25![6,.(k'), c.(k)]_ =0 ’ 


T'—o 


and similar ones for Fermi particles and antiparticle operators. The important 
point is that the limits to the remote future (or past) for each operator in à 
product have to be taken independently. 

For the present work, before taking the limit to remote future or past 
we modify the instantaneous creation annihilation operators themselves in a 
way so as to conform more with the actual measuring technique. Such pro- 
cedure yields the desired commutation relations also. 


3. — Commutators of instantaneous operators. 


Before introducing the modified creation and annihilation operators, we 
want to obtain in a convenient form the commutators between instantaneous 


(4) G. Mowan: Suppl. Nuovo Cimento, 5, 440 (1957). 


32 - Il Nuovo Cimento. 
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operators that refer to different instant parameters 7 and 7°. These relations 
will be used in the next section. 

By means of the Gauss theorem we can relate the corresponding operators 
belonging to the two different instant parameters t and t’. For example, 


— OF (k) + eî(k) = i farai fr(æ) Da] = 


0%, 


To 


or 


(6) cé(k) = cé(k) — i fate) 7 


Similarly we have 


(7) ot pp) = P_i [do Tea)glps a), 
and other relations. The product a,.(p'8')a*(pf) can now be expressed in a 
symmetrical fashion which will subsequently yield a simple expression for 
the anti-commutator of a,(p’p’) and a* (pf). 

a (pB) — à [ave Te)ug(ps0)| — 


To 


(8) a; (p'B')« (pb) = dx (p'p') 


T 


= a,(p'B')az(pR) —iar(p'B') ‘ d'2I(2)up(p; 0) + 


To 
sr 


a | dtyitp (p's y) Ly) aX(pB) - 
The product 1(y)a*(pB) in the last integral can ee in a different form: 
(9)  I(y)aë(pB) = Ty) a¥ (pB) + in) fase ram OSs 
= [I(y), aÿ (ph) — a (DOT) = 
i i dt T(x) Ug(p; æ)1(y) +iI(y) fi dt T(a) ug(p; 0) = 


y Ta 


T 


“iaia | die T {I(y) F(a)} up(p; 2) - 


Te 
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Substituting eq. (9) into eq. (8) we obtain 


(10) a (p'B')a (pp) = a (p'B')az(pP) — 


T 


— ia (p'B') dir La) (p 2) — iat (pp) | ty us (p'; y) L(y) — 


= | ty fair Le Tout: a) + 


+ à [ave (ps DCI) a) 
A similar expression can be dérived for the product a*(pBja, (p'B'). Adding 


these two expressions we obtain 


(11) [a(p'87), a7 (pB)]+ = [æ,(p'8"), a (pB)} + 
= i fav y[L(y), a (ph): — i fasti), a,(p'B')|+Ug(p5 2) - 
All other commutators between various annihilation and creation operators 


can be derived by this method. Some of these are given below: 


(12) [a,(p'B'), 4:(pP)]+ = 


Ù 
T T 


= à fareitp; « U( ap + à favi oa 


To To 


(3) [ex(k°), oF (be) |. = [02 (4°), (k)L + 


+i faut), UL + à Janet e. 
Expressions for other commutators have essentially the same structure. Sub- 
stituting for Z(æ) from eq. (1) into eqs. (11) and (13) we obtain 
(14) [a (p'B"), a (pb), = [a.(p'8"), a (pB)h + i [ate Np pal 


t' 


T 
a 


(15) [er (k), cé(k)]_ = [e,,(h'), &(k)]_ + i fare Gen (&) , 


Tel 
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where 
N ,8.p8(@) = Up(p'; {dm — gy, Da)}u,(p; 2), 
Gi, (0) = fe(0)fx(0){ dx? + 32 D(2)}. 


A significant result that follows from these calculations is the vanishing of 
the following anticommutators for arbitrary 7 and 7’: 


La, a, }+ =0, 
(16) [b;, 8.1: =0, 
[a,, b*), n= 0 x 


Now we are ready to introduce the modified operators and obtain their 
commutation properties by the help of the results of this section. 


4. — Modified creation and annihilation operators. 


It is generally believed that the creation operators of pyhsical particles 
are obtained by taking suitable limits in the definitions (3) to remote future 
or past. In any expression involving product of such operators, such as the 
left hand side of eqs. (5), the limit for each operator must be taken inde- 
pendently. One might expect that this has the effect of making the reference 
to the instant parameter diffuse and thus, at least partially, take account of 
the non-vanishing duration of the measuring processes. 

Our procedure is somewhat different. We first modify the definitions of 
the creation-annihilation operators themselves so as to be characterized by 
an extended time interval and then take the limit to remote past or future. 
The purpose of the limit is then only to avoid the insignificant transient effects. 
How exactly the time extension should be incorporated can be answered, in 
our opinion, only by a future well-developed quantum mechanical theory of 
measurements. Presently, we propose an extremely simple modification which 
undoubtedly can only be an idealization of completely realistic creation ope- 
rators (*): 


(17) a*(pp) = lim fac exp[ia,t] az (pl), 
May +7) 0 Ta — T1 


Ti 


(*) We shall in the following study the limit to remote future alone. The case of 
remote past is very similar. 
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To 


il 
b*(pp) = lim ——— {dt exp [ta,7] br (pp), 
(Ti +72) = 


Ta 


ai 
c*(k) = lim ———|drexp{ixr|c*(k), 
Hr, +7)>o0 Tz — Ti 
Ti 
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where x, and «, are real infinitesimal functions of the parameters of the theory 
and also of p and k respectively. These functions will be further specified 
presently. The interval (7, — t;) should be connected with the duration of 
The exponential factor in the integrands can be interpreted 


measurement. 
as altering the energy in the functions w(p; x), v(p; x) and f,(x) 


finitesimal amount. 


We now compute the anticommutator of a(p'f') and a*(pB). From eq. (1 
we have (*) 


(20) 


by real in- 


4) 


T1 


— exp [1(a,”° — Ta) + exp [i(oiti —04T2)] — exp [1(%,T2 —0T;)]) 


a*(pB)|+ = 
= lim 2 for [ae exp [i(a,t — ot) ][a,(p'B"), a*(pB)l- 
REA 1) 
ZE {4 exp [43(% — ot4)(T + T)] Sin Jan (t— 71) sin $oa(z Une Ti) 
(Ti + Te) > 061 04 (Ta lè Ca) 
[a,.(p'9); aB(pA) + — — | de po pol) >" 
Sn cai 


‘(exp [— i(a10° — at,)] — exp[— i (c10° — t,) + exp [1(a,4° — a,7,)] — 


Let M(x°) denote the quantity within the bracket in the integrand, An 
arbitrary matrix element of the anticommutator between energy momentum 


eigenstates will involve 


easily evaluated. The 


only through a factor 
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exp [— i(0°+a —%)4(t%1 + 7). 


an integral of the type i dix exp [1Qx|M(a#) which is 


4(t,+7,) dependence of the resulting expression is 
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Now if we let 4(7,1+7,) go to +co, then unless there is a factor singular at 
Q° =a,— «x, we can replace the above factor by 6,,,,_,;- Putting Qo =a,—«, 
and retaining only the lowest order terms in «, and «, we find that 


To 


fato exp [iQx] MU (a) = — 


Ty 


a 
(=)? 


(27)? 
12 


(21) d(Q)(a, +a1)(T — Ti) Omen © 


From eq. (21) we find that to the lowest order in «, and «| the matrix element 
of the anticommutator between any two energy-momentum eigenstates |P’y'> 
and |P"y"> is given by the following equation: 


(22) <P"y" |[a(p'B'), a*(pB)}+|P'y) =[a,(p'8"), az(pB)]+<P"y" |P'y> — 
a fini 5 ee ne Pap 
MT: — T1)? 
12Vp°p" 


-0(p + P'— p'"—P")-ug(p'\P"y"|9gysP(0)|P'y)u;(p) - 


; ; 
& + Xi) Op? + Pp 0_ po, a, a" 


So far we have made only one assumption; that is, «, is vanishingly small. 
Now the objective is to make it vanish in such a way that 


2, [a(p'8"), a*(pB)l: = 65 y dp — p’) - 


Consequently the last term on the right hand side of the eq. (22) must vanish. 
Further, since 0 < 2, < 1, the second term on the right hand side must not 
vanish. This means that for the applicability of this method dm, and simil- 
arly 0x? must necessarily be infinite. Following properties of «, can now be 
concluded: 


(23) a(p) = PM ; 


where 7, is independent of p, 


ES, 2 
(24) dal {1 = Moni M om} = 
and 
(25) 2, ng = 0. 


A similar consideration for the vanishing of the commutator [a(pf), c*(k)] 


00 
(cr) 
> 
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leads to the condition 


(26) 2,2, ng =0. 


The left hand sides of the eqs. (24), (25) and (26) are supposed to be defined 
through a limiting procedure in which 7, goes to zero and dm goes to infinity. 

With a better understanding of the magnitudes and relationships between 
the renormalization constants dm, dx?, 2,, 2 and 2, it is hoped that an explicit 
expression for «, and also for x, will be obtainable. As it is, full implication 
of the eqs. (25) and (26) need to be further studied. 

Similar relations for the function «, can be obtained from the commutators 
of boson creation annihilation operators. 

Thus we find that functions a, and «, can possibly be so chosen that the 
operators 2, *a*(pB), <;*a(pB), 2;*c*(k), ete., have the typical commutation 
properties of annihilation and creation operators. In the next section we show 
that they annihilate or create physical particles. 


5. — Physical particle states and S-matrix. 


Let us consider as an example the state vector 2;!a*(p8)|0). If it re- 
presents a physical Fermi particle with momentum p then we must have 


(27) Pres ta*(pB)|0)= p"z;'a*(p8)|0 
so that 
(28) — P#P 2, ‘a*(pB)|0> = m2, *a*(pB)|0> . 


From the definition of a*(pB) and noting that 


Wie) = exp[— iPx] (0) exp[iPx]|, 


it is easily seen that <P’y'|\z;!a*(pB)|0> vanishes unless P'= p and P'= 
= p°= Vp*?+ m*. This implies the eqs. (27) and (28). 

The renormalized S-matrix elements are constructed from the well known 
method of reduction formula (5). 


(5) See reference (1). 
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6. — Conclusion. 


We have introduced à set of modified creation and annihilation operators 
defined by the eqs. (17), (18) and (19). They differ from the conventional 
operators by incorporating a reference to extended time in their definitions, 
and thus there is a closer correspondence with actual measurements. It is 
shown, independently of asymptotic assumptions and provided certain consis- 
tency requirements can be met, that these obey the typical commutation re- 
lations of creation and annihilation operators, and they annihilate and create 
physical particles. There is no need for adiabatic hypotheses or asymptotic 
assumptions in the construction of multiple physical particle states or the 
renormalized S-matrix from these modified operators. However the analsyis 
of the functions «, and «, which go to define the modified operators is in- 
complete. We expect that a study of «, and « might be related with the 
renormalizability of the theory. 


IRAN SAW INE) Cea) 


In vista del fatto che tutte le misure richiedono un tempo non scomparente, si modi- 
ficano gli usuali operatori di distruzione e di creazione. Gli operatori modificati creano 
e distruggono particelle fisiche. Non si invocano ipotesi asintotiche o adiabatiche, 
cosicchè si può ottenere la matrice S rinormalizzata senza il loro ausilio. L’analisi di 
talune funzioni che compaiono nelle definizioni degli operatori modificati è incompleta. 


(*) Traduzione a cura della Redazione. 
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Decay of *’*Ku and the Unified Nuclear Model. 


S. K. BHATTACHERJEE, T. D. NAINAN (*), SHREE RAMAN and BALDEV SAHAI 


Tata Institute of Fundamental Research - Bombay 


(ricevuto il 29 Agosto 1957) 


Summary. — The radiations from !5?Eu (t;=13 yrs) were investigated 
using an intermediate image spectrometer, a magnetic thin lens spectro- 
meter and a coincidence scintillation spectrometer. Three §~-groups are 
present with end-points (1450 +25), (1 040 +20), (710 +20) keV; the highest 
energy group exhibits an «-shape. y-ray investigations showed that the 
following y-rays are present: 101 (5), 121 (88), 244 (30), 315 (4), 344 (100), 
410 (9), 442 (4), 710 (13), 770 (77), 888 (14), 970 (89), 1090 (96), 1127 (96) 
and 1410 (135) keV. K-conversion coefficients were estimated for most 
of these transitions and multipole orders were assigned. X-ray-y, Y-Y 
and f-y coincidence experiments allow us to construct the decay scheme 
of 152Ku. The difference in level structure in 152Sm (N=90) and Gd 
(N =88) is evident from the experimental results. The first three excited 
states at 121, 365 and 807 keV in 152Sm are identified as the three members 
of the ground state rotational band (K =0); the levels at 1090 and 1248 keV 
having 1—2+,and 3+ respectively are shown to be the quadrupole y-vibra- 
tional states corresponding to K=2. The highest excited state in 152Sm is 
at 1417 keV and is most likely fed entirely by L-capture; this is the only 
odd-parity state in ?Sm having the character /=1—. The levels in 
152Gd at 344 and 754 keV both with J=2+, reached by L--decay, are 
characteristic of the collective vibrations of the electric quadrupole type 
about a spherical equilibrium shape. The level at 1114 keV in Gd 
has the character J=3— and is interpreted as due to the octupole 
vibrations. The ground state of 15?Eu has been assigned the spin J=4 —; 
the branching ratios and the transition probabilities for the 8--groups 
and electron capture transitions have been calculated. The experimental 
results are consistent with the predictions of the Unified Nuclear Model. 


(*) Present address: Department of Physics, University of Indiana, Bloomington, 
Indiana, U.S.A. 
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1. — Introduction. 


The 13 year !Eu has been shown by several workers (14) to decay by 
electron capture to !*Sm and by B--emission to Gd. The recent works of 
NATHAN and WAGGONER (3) and of BHATTACHERJEE and SHREE RAMAN (8) 
have clearly brought out the importance of the study of this isotope from the 
point of view of the Unified Nuclear Model (5). A change in the nuclear struc- 
ture in the neighborhood of N = 90 has been indicated by anomalies in the 
isotopic shift (6), quadrupole moments (7) and Coulomb excitation energy 
levels (8). Several abnormalities (*) exist in the binding-energy data in the 
vicinity of the neutron number N = 90. MOTTELSON and NILSSON (1°) em- 
ployed an ellipsoidal potential with a spin-orbit force to calculate the indi- 
vidual-particle energy levels as a function of nuclear distortion. With these 
levels, the ground state equilibrium deformations were calculated for nuclei 
in the region N = 82 to N —126. A sharp increase in nuclear deformation 
was predicted in going from N = 88 to N = 90. One encounters in the decay 
of !5*Eu different types of nuclear spectra as one crosses N= 90 at °%Sm to 
N — 88 at 1Gd. In !*Sm well-developed low-lying rotational states exist 
corresponding to collective excitations of the strongly deformed nucleus; on 
the other hand, in !52Gd the low-lying excited states can be classified as vibra- 
tional levels corresponding to the shape oscillations of a spherical nucleus (°). 
The recent works (**), as mentioned earlier, on this isotope were carried out 
with the mixed activities of 13 year !*Eu and 16 year Eu as obtained 
from the neutron irradiations of natural europium. Because of the similarity 
in the decays of these two isotopes it was found difficult to make proper iso- 
topic assignments of some low-intensity y-rays. 

The present paper deals with the results of our work on the 13 year !52Eu 
as obtained from the thermal neutron irradiations of 92 percent electro- 
magnetically enriched isotope of 51Eu. In this work attempt has been made 
to determine the relative intensities of the different y-rays and the internal 
conversion lines; this enables us to estimate the K-conversion coefficients of 


) R. E. SLATTERY, D. C. Lu and M. L. WIEDENBECK: Phys. Rev., 99, 1615 (1955). 
) L. Gropzins: Bull. Am. Phys. Soc., 1, 163 (1956). 
(3) O. NATHAN and M. A. WAGGonER: Nuclear Phys., 2, 548 (1957). 
) S. K. BHATTACHERJEE and S. RAMAN: Nuclear Phys., 4, 44 (1957). 
) K. ALDER, A. Bour, T. Huus, B. MortELSON and A. WINTHER: Rev. Mod. 
, 28, 523 (1956). 
(6) P. Brix and H. KOPFERMANN: Phys. Rev., 85, 1050 (1952). 
(7) P. Brix: Zeus. Phys., 132; 579 (1952). 
(8) N. P. HEYDENBERG and G. M. TeMMER: Phys. Rev., 100, 150 (1955). 
(9) W. H. Jounson and A. O. NieR: Phys. Rev., 105, 1014 (1957). 
ga) 


502 


: 
f 
| 
| 
| 
f 
' 


Cane 


DECAY OF l?Eu AND THE UNIFIED NUCLEAR MODEL 


503 


most of the y-rays and thus assign the multipole orders. The relative transition 
probabilities of electron capture and B--emission were calculated from the 
measured branching ratios and the experimental results were then compared 
with the predictions of the Unified Nuclear Model (5). 


2. — Experimental procedure. 


Reactor. 


Electromagnetically enriched isotope of *!Ku (92 per cent), obtained from 
the Stable Isotope Division, Oak Ridge, U.S.A., was bombarded by thermal 
neutrons for à. period of 3 weeks at a flux of 10! n/cm?/s in the Oak Ridge 


of the sample from the 
reactor. The whole exper- 
imental procedure con- 
sists 
measurements : 
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1) 


‘tal 


of the following 


Low energy elec- 
tronspectrum with 
the Siegbahn-Slatis 
intermediate ima- 
ge spectrometer. 


High energy elec- 
tron spectrum with 
a thin magnetic 
lens spectrometer. 


y-rays energy and 
relative intensity 
measurements with 
a calibrated 2 in. 
thick x 13 in. dia- 
meter NalI(T1) crys- 
spectrometer. 


X-ray-y, y-y and 
B-y coincidences 
with a coincidence 
scintillation spec- 
trometer. 


Experiments were started after a month from the date of removal 


20 


He (Gauss-cm) 


0 337 674 1011 1348 1685 2022 2359 


Fig. 1. — Low energy electron spectrum of ?Eu taken 

with the intermediate image spectrometer. Energies, 

isotope assignments, and relative intensities are sum- 
marized in Table I. 
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2°1. Low energy electron spectrum. — The Siegbahn-Slatis Intermediate Image 
spectrometer of the Institute was employed for this measurement. The source 
was deposited on 150 ug/em? Al-foil, the Geiger-Müller counter window con- 
sisted of 75 ug/em? collodion film. The momentum resolution was about 
4.5 per cent. The spectrum is reproduced in Fig. 1. For the low-energy 
electron lines the KL differences in Sm and Gd allow identification of the 
product nucleus. 

In Table I these results are summarized along with those obtained with 
the Lens spectrometer. The isotopic assignments are made mainly on the 
basis of our coincidence results and our assignments agree well with those 
of Cork et al. (11). The Z-line of the 101 keV y-rays falls at the low energy 


TABLE I. — Energies and relative intensities of internal conversion lines (*). 
Elec- Relative Iso- Corres- 
tron Inter- Re tope ponding K/L Spectrometer 
energy pretation sic assign- | y-energy employed 
(keV) i ment (keV) 
33.1 | Auger-KLL 700 Sm — cs Intermediate image 
39.2 | Auger-KLM 110 Sm — — » 
54.2 K 220 Sm 101 ms 2.3 ) 
74.5 K 2355 Sm 121 w 1.25 » 
95 L 95 Sm 101 » 
114 L 1900 Sm 121 » 
160 (+) K 16 Sm 207 » 
197 K 83 Sm 244 rw 3.2 Thin lens 
218 (+) K 2 Sm 265 » 
236 L 26 Sm 244 » 
243 (+) Il, INS 5 Sm 244 » 
268 K 4 Sm 315 » 
294 K 100 Gd 344 a 3.5 » 
336 L 29 Gd 344 » 
362 K 4.5 Gd 412 » 
396 K 6 Sm 443 ÿ 
656 K 4.5 Sm 703 » 
69914) K 1 — — » 
724 K 2.5 Gd 774 » 
926 K 5 Sm 973 » 
1047 K 6 Sm 1094 » 
1080 K 4 Sm 1127 » | 
170) K 2 Sm 1417 » 
(*) The energies are generally believed to be correct to within 1.5 % and the intensities to 
within 15% except when otherwise stated. 
(*) Since the lines are weak and broad, rather large errors must be assigned to their mea- 
surements. 


(11) J. M. Cork, M. K. Brice, R. G. HELMER and D. E. Sarason: Bull. Am. Phys. 
Soc., 2, No. 1, E316 (1957). 


504 


Obi; 


ee 


pm ee n 


DECAY OF !52Eu AND THE UNIFIED NUCLEAR MODEL 505 


tail of the L-line of 121 keV y-rays because of poor resolution, this renders 
the intensity measurements of both these lines rather inaccurate. A conver- 
sion line corresponding to a y-ray of 104keV has been reported by Cork 
et al. (4) and they have assigned it to !5*2Gd; our coincidence measurements, 
however, favour its origin in 1?Sm. 


2°2. High energy electron spectrum. — Measurements were carried out with 
the thin Magnetic Lens Spectrometer (2?) with a momentum resolution of 
about 2.4 per cent. Sources, were deposited (about 3mm in diameter) on 
0.6 mg/cm? rubber hydrochloride film with the source properly grounded by 


0.9 


Ho (Gauss-cm) 


0.3 


Ho (Gauss-cm ) 


1300 1600 2300 2800 3300 3600 


Fig. 2. — High energy electron spectrum of !52Eu with the magnetic lens spectrometer. 
Energies, isotopic assignments and relative intensities are summarized in Table I. 


aquadag lines. The Geiger-Müller counter window consisted of 1.6 mg/cm? 
mica foil. The spectrum is shown in Fig. 2, and the results are summarized 
in Table I. Besides the well resolved low-energy lines, some strong conversion 


(2) T. D. Nainan, H. G. DEWARE and A. MuKERJI: Proc. Ind. Acad. Sc., 44, 
111 (1956). 
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lines are present in the higher energy range. These higher energy lines are 
situated on a large continuous background and hence intensity measurements 
of these lines involve 
rather large errors due to 
this tailing. 

The total f-spectrum 
was analysed by means 
of a Kurie plot and was 
found to be complex 
(Fig. 3). The upper ener- 
gy portion of the uncor- 
rected Kurie plot has a 
slight curvature concave 
towards the energy axis 
characteristic of a unique 
first forbidden spectrum. 
When the unique correct- 
ion factor S,=—(W?—1)+ 
+(W,—W)? is applied, 
the Kurie plot is straight- 
ened out, thereby giv- 
ing the end-point of the 
highest energy f_-group 
as (1450425) keV. The 
zi corrected Kurie plot starts 
Reema, | "a 1450 Lev to deviate from a linear 

DITTE plot near 1000 keV; when 


Fig. 3. — Kurie analysis of the B--spectrum. e---@---@ Le ing meet SEL bra 
Gross Kurie plot. o...0...0 Corrected for the unique forbidden spectrum’ is 
shape for the highest energy component. The correction Subtracted from the total 
factor employed: Sy= (W2—1)+(W,— W)?. A...A..4 spectrum, a second Kurie 
‘ The 710 and 1050 keV 8° -groups. +...+...+ Thelowest plot is obtained which 
energy B--component at 710 keV. 


\710 keV 
1 2 3 4 


gives the end-pointenergy 
of 1040 keV for the se- 
cond £--group. This second Kurie plot starts to deviate from around 
700 keV and subtracting from the second plot, we obtain the third 6--group 
at end-point energy of 710 keV. The energy differences between the three 
B--groups indicate that the 770 and 410 keV y-rays are associated with !52Gd 
along with the 344 keV y-ray. Our determination of the maximum energies 
of three f--groups agrees very well with that of CORK et al. (11). 

The area of the total B--spectrum was analysed into its three components 
by means of the Kurie plot analysis. The branching ratios of these three groups 
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are shown in Table II together with the results as obtained from the measured 


y-ray intensities. The disagreement for the highest energy group is a little 
too large. However, considering the uncertainties involved in extracting the 


TABLE II. — fB-group maximum energies and branching ratios. 
B-energy i Branching ratio Branching ratio from | 
(keV) (from Kurie analysis) y-intensities | 
| 
1450 27 18 
1040 9 8 | 
710 64 74 | 


branching ratios from Kurie analysis, the intensities as derived from the y-ray 
measurements are believed to be more reliable and hence these intensity values 
will be utilized in evaluating the transition probabilities in the 8--branching 
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Fig. 4. — y-spectrum of !5*Eu in 2in. thick x 1?in. diam. Nal(Tl) crystal. The energies 
and relative intensities are summarized in Table III. Energies are in keV. 
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of Ku. It has been established by HILL and SHEPHARD (1) that !52Eu does 
not emit positrons, as is also evidenced by the absence of any annihilation 
radiation in the observed y-spectrum (Fig. 4). 


23. Energies and relative intensities of the y-rays. Singles spectrum ana- 
lysis. — The single-channel spectrometer employed a 2 in. thick x1?in. dia- 
meter NaI(T1) crystal and a Dumont 6292 photomultiplier. The peak effi- 
ciency curve for this crystal at a source-to-crystal distance of lin. was de- 
determined (Fig. 1 of ref. (‘)) for various monoenergetic y sources. The re- 
solution was about 9 per cent for the 661 keV line of !87Cs. The pulse height 
distribution generated by the !52Eu y-rays is shown in Fig. 4, and the y-ray 
energies and their relative intensities are summarized in Table III. 


TABLE III. — Observed y-ray energies and their relative unconverted intensities (*). 


y-ray Energies Relative Intensities 
121 74 
244 24 
344 100 
410 (+) 9 
710 (+) 10 
770 62 
870 (+) 11 
970 71 
1090 154 
1410 108 


(*) The energies are believed to be correct to 2% and 
intensities within 15 % unless otherwise stated. From the 
considerations of thermal neutron absorption cross-sections 
of 1Eu and ‘Eu, the isotopic enrichment ratios (92% 
of 1Eu), and their lifetimes one can conclude that the 
interference due to the radiations from ‘Eu present is 
practically negligible. 

(*) These lines are not well resolved in the single 
spectrum and hence their energies and relative intensity 
measurements are subject to larger errors. 


204. X-ray-y and y-y coincidences. — The coincidence experiments were 
carried out with a Coincidence Scintillation Spectrometer which is described 
elsewhere (4). It utilizes the « fast-slow » principle, the fast coincidence re- 
solving time was 6-10-8 s. For the y-y coincidences both the detecting crystals 
were 2in. thick x1?in. diameter and 6292 photomultipliers were employed. 


(13) J. M. Hrcz and L. R. SHEPHARD: Proc. Phys. Soc., A 68, 126 (1950). 
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All the coincidence experiments were carried out with the counter axes oriented 
at 90° with a view to avoid spurious coincidences; an anti-scattering shield 
was used in the usual ray. The ratio of the accidental coincidence rate to the 
genuine coincidence rate in the region of the true coincidence peaks was usually 
small. The results are summarized in Table IV and a few representative co- 
incidence spectra are shown in Fig. 5, 6 and 7. 


400 


X-X COINCIDENCE 


= 
à 
a 
da) 344 
264 SINGLES 
100 
COINCIDENCE 
AT 
0 L 
0 22:5 45 67:5 
PULSE HEIGHT (Vv) 
ig. 5. — K X-ray coincidences with the low energy y-rays in 152Eu. 
ig. 5. 3 y y-ray 


One of the features of our coincidence results is the presence of weak 
y-rays of energies (in keV) & 108, 315, 410 and 442; all these lines are also 
shown up as weak conversion lines in the 6--spectrum (Fig. 2). These low- 
energy coincidence peaks cannot be explained as spurious coincidences in view 
of the geometry used in these experiments. The relative intensities of these 
lines could be roughly estimated from the coincidence data; the energies quoted 
above are appropriate to the observed conversion lines (Table I). Within the 
limits of experimental errors the energies of the coincidences peaks as shown 
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in Fig. 5, 6 and 7 agree well with the above-mentioned energies for these y-rays. 
Such an observation was also made in our earlier work (*) with the mixed 
activity. In coincidence with the 121 keV y-ray, the relative intensity of the 


h 
| 
i 
| 
| 


400 


300 


C.P.M. 


200 


K X RAYS 
244 


100 


È 
9 15 30 45 60 75 D 9 18 27 36 
PULSE HEIGHT CV) 


Fig. 6. — (a), (b), (c) show the coincidences of !5°Eu y-rays with 121 keV photopeak. 

(d) represents a portion of the spectrum in coincidence with 244 keV photopeak. Note 

that the unresolved peaks in the single spectrum at 710 and 870 keV are completely 
resolved here. 


1090 keV y-ray is reduced by 50 per cent as compared to that of 970 keV, 
if the relative intensity ratio of these two y-rays from the singles spectrum 
is taken as standard. This establishes the presence of the y-ray of energy 
1127 keV which is of equal intensity as that of 1090 keV; these two lines were 
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not resolved in the singles spectrum, the former is in coincidence with the 121 keV 
y-ray, while the latter goes directly to the ground state from the 1090 keV 
level in 15°Sm. 


30 
330 770 
- 
20 
L 
4 
240 a 
= 
Q 
G af 
10 
PULSE HEIGHT (V) 
0 i dread 
0 22:5 45 57:5 


Fig. 7. — The y-spectrum in coincidence with 344 keV photopeak. Note that the un- 
resolved line at 410 keV is completely resolved here. No appreciable coincidences 
were observed for y energies higher than 770 keV. 


It is significant to note that no KX-ray or 121 keV coincidences were ob- 
served with the 1417 keV y-ray. In the earlier works (*') with the mixed 
activity the coincidences with the 121 keV y-ray attributed to the 1410 keV 
y-ray are most probably due to !5*Eu. This fact has been verified (1) with an 
enriched activity (95%) of 15iu and we have found coincidences between 
123 keV and the higher energy y-rays present in this isotope. 


2°4.B B-y coincidences. — An anthrancene crystal of jin. thickness and 
14 in. diameter was used to detect the f-particles. This crystal was covered 
on the top by 150 ug/em? Al-reflector to ensure efficient light collection; the 
resolution obtained was about 18 per cent for the 624 keV Æ-conversion line 
of the ?’Cs y-ray. The full photopeak at 344 keV was accepted in the y-ray 
detecting gate channel and the coincidences B-spectrum was obtained by 


(14) Work on 16 yr !54Eu decay is in progress and will be reported in due course. 


all 


121, 244, 710, 
870, 970, 1090 
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EELS IWS = SUEY Lies eis si and ey Oa 
o 9 C oc | 
radiation | radiations Remarks 
(keV) (keV) 
39 K-X-rays | K-X-rays, 121, Strong X-ray coincidences indicate that in the 


capture branch (425m) some low-lying states 
de-excite by highly converted y-rays. No coin- 
cidences were observed with 344 keV indicating 
that this y-ray is in bulk in the B_-branch. 


of equal intensity. 


branch is also confirmed. 


Weak coincidences ata ~ 108 keV confirm our 
earlier observation that the 121 keV y line is com- 
posed of a y-ray of energy close to 108 keV and 
must be present in the capture branch. Similar 
arguments hold for weak 340 keV coincidences. 
No appreciable coincidences were observed near 
the 770 keV photo peak. The relative intensity of 
1090 keV photons per 121 keV photons is about 
50% of that in singles spectrum showing that the 
1090 keV peak is composed of two y-rays, each 


The unresolved peaks near 770 keV in nol 
spectrum are completely resolved in coincidence | 
spectrum showing that the 244 keV y-ray is in 
coincidence with the 710 keV and 870 keV y-rays 
but not with the 770 keV y-rays. 
y-ray 


The existence 
= 340 keV in the capture | 


The 410 keV unresolved peak in the singles 
spectrum is clearly resolved. Negligible coinci- 
dences in the A-X-ray region indicate that | 
344 keV y-ray in bulk must be present in the 
No appreciable coincidences were | 
observed i in the energy ranges higher than 770 keV. 


The X-rays, 121 and 244 keV coincidences 
are presumably due to 870 keV y-rays present | 
in the composite peak at 770keVin the gatechan- | 
No appreciable higher energy coincidences | 
were observed after the intense 344 keV peak. 


121 K-X-rays, | 
n 108, 244, 
340, 450, 970, 
1110 
244 | K-X-rays, 121, | 
> 240, 340, 450, | 
| 710 and 870 
| 
of a weak 
344 108, 240, 330, 
410 and 770 
B- branch. 
770 K-X-rays, 121, | 
240, 344 | 
| nel. 
970 K-X-rays, 121, 
a 240 
1090 K-X-rays, 121, 
240, & 330 | 
ISO |-__ —__ + 
1410 No coincidences | 


state in 152Sm. si 


coincidence with the 121 keV Y-rays. a 


This y-ray presumably leads to the 121 keV 


Coincidences with 121 keV Yy-Trays indicate 
that it is indeed the 1127 keV y-ray which is in 


The 1410 keV y-rays did not show any ap- 
preciable coincidences either with K-X-rays or | 


with the 121 keV photopeak. 
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scanning the anthrancene crystal spectrometer. Fig. 8 (a) represents the 
Spectrum and one clearly sees two B-groups of end-point energies at 1500 
and 750 keV in coincidence with the 344 keV y-ray. In Fig. 8 (5) is shown 
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Fig. 8. — (a) B--spectrum of !92Eu in coincidence with 344 keV y-rays, the two 8--com- 

ponents at 750 and 1500 keV end points are clearly discernible. (b) B--spectrum of 

152Ru in coincidence with 770 keV y-rays. The only coincident £--group has the end 

point energy of ~ 750 keV. The resolution of the anthracene crystal spectrometer 
for the 624 keV K-conversion line of ’Cs is & 18%. 


the coincident f-spectrum when the full 770 keV photopeak was in the gate 
channel; the only f-group coincident with the 770 keV y-rays had the end- 
point energy at ~ 750 keV. In all these B-y coincidences, the y-y background 
was separately determined and only the genuine {}-y coincidences are ploted 
in Fig. 8. 

In order to ascertain the y-rays that are present in the 8--branch of !°2Eu, 
a band of B-particles of energy between 400 and 440 keV was selected as the 
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gate radiation and the coincident y-spectrum is reproduced in Fig. 9. The 
two y-rays at 344 and 770 keV are the ones that are strongly in coincidence, 
no y-ray of energy higher than 770 keV shows any appreciable coincidence. 
The weak peak at about 410 keV is not well-resolved. The coincident peak at 
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0 25 50 WS 


Fig. 9.— The y-spectrum of 5?Eu in coincidence with a band of B--partieles of energies 
between 400 and 440 keV. This band was particularly chosen as not to include the 
conversion electrons due to intense low energy y-rays. No appreciable B-y coincidences 
were observed with the y-rays of energies higher than 770 keV. In all B-y coincidences 
the y-y background was observed by absorbing all f--particles by means of a 3/16 in. 
perspex absorber. The spectra shown are the genuine 8-y coincidences. 


= 150 keV is presumably due to the back scattering of the 344 keV Y-TAYS, 
since this experiment was done at 180°-geometry. The result of the ®-y co- 
incidences are summarized in Table V. 
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TABLE V. — Summary of B-y coincidence measurements. 
| Gating Information 
| radiation about the coin- Remarks 
(keV) cident radiations 
| 344 keV B--groups of The 344 keV y-ray is in coincidence with 
. photopeak > 1500 and two _-groups as is clearly shown in Fig. 8 (a). 
| | w 750 keV 
| . 
| end-points. 
| 770 keV B--group of The only B--group which is in coincidence 
| . . pa 
| > 750 keV with the 770 keV y-ray is the one of 710 keV 
end-point. end-point energy (Fig. 8 (b)). 


6--particles in 


344 and 770 keV 


This experiment clearly demonstrates (Fig. 9) 


the energy photopeaks that the strong y-rays that are present in the 
band 400 = B--branch of !5*Eu are 344 and 770 keV. No 
— 440 keV appreciable =z~-y coincidences were observed 
accepted in for any y-rays higher in energy than 770 keV. 
anthracene The slight coincidences near & 150 keV might 
crystal be due to back scattering of 344 keV y-rays 


since this experiment was performed with 
| counter axes oriented at 180°. 


3. — Multipole order assignment to the y-rays. 


From a knowledge of the relative intensities of the A-conversion lines and 
the y-rays, it is possible to estimate the relative magnitude of the K-con- 
version coefficients of the different y-rays. The relative values then can be 
converted to absolute conversion coefficients if the multipole order of one of 
the transitions is known in an independent way. The {}-spectrum analysis 
and the 8-y coincidence experiments establish the fact that the 344 keV y-ray 
arises due to the de-excitation of the first excited state in Gd and the 
maximum f--group of 1450 keV terminates at this level; so it is concluded 
that no ground state B--transition takes place in !52Gd, hence the 344 keV y-ray 
intensity represents the total number of B-disintegrations from !?Eu. The 
K-conversion coefficient «, of the 344 keV y-ray can then be determined 
from the following relation, 


I ie 
or let Ts 1) (1 pr aid | 


The conversion contribution due to shells higher than Z is neglected in this 
expression. 1,=1,+1,+1,, 1, and I, represent the intensities of the 
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K and L conversion lines and Æ/L is the measured ratio. From the lens spectro- 
meter measurements (Fig. 2) it is found that the total y-continuum intensity 
is 3390 as compared to 100 for the K-conversion line of 344 keV. Hence 
from the other measured data (Table I) ie. I,=29 and K/L= 100/29 we 
get the value of the K-conversion coefficient for the 344 keV y-ray as 
a, = 0.032 + 0.005. The theoretical K-conversion coefficient for a pure H2 
transition is 0.033 and is in excellent agreement with the measured value. 
The K-conversion coefficient of 344 keV can now be used as a reference for 
calculating the absolute AK-conversion coefficients of the other y-rays. Table VI 
summarizes the experimental conversion coefficients calculated in this way; 
the theoretical K-conversion coefficients for the different multipole orders are 
obtained from RosE’s (5) tables. 


TaBLE VI. — K conversion coefficients and multipole order assignments. 


| 
Iso- | Lu IE ao 
De intensity| verted OK a i È Pata Proposed 
energy assign- of K-con- inten- Experi- | Multipole 
(keV) CRE version sity mental El | #2 Wal Order. 
lines (d) | 
——si = — = + sà | 
1014) Sm 2 DUREE 1-32 (RR 1.6 E2 
121% Sm | 2900 || Bk SCONTO cs | 6.2-10=1 _ | E2 | 
244 Sm 83 MIS OMS SOA 102188 00 A2 Om E2 | 
315°| Sm 4 4 | 3.2-10-2] 1.2-10-2 | 4.1-10-2| 8.0-10-2 | (#1 or E2)° 
344 Gd 100 LOO MGS. 2510-44 2.0810 A 8 MO 6:60 a4 E2 
410 © Gd 4.5 9 Or Om" NO A0) 20102 EA E2 | 
442 © Sm 6 4 | 45:10?) 4.9-10-3| 1.5-10-2| 2.7-10-2| (#2+ M1) | 
710 Sm abel | E 1.0-10-2| 1.8-10-3| 4.6-10-3 | 9.0-105 | (#2+ M1) 
770 Gd 2.5 UNI MOTOS ANG 103 12100 E1 
870 Sm — 14 e —. — | == — 
970 + Sm | 5 89 BARON 9.0 1054 182.2 LOS 3:97 1078 E2 
| 1090 ¢ Sm 6 96 OCTO AIO RIZZO | 3.0.-107=8 E2 
1276 Sm 4 JOM Se O see O LOSS RIGO enc ol Oma H2 
|} 1417 | Sm 2 | 135 (4410-447104) 10-10 | 16-101 F1 
(a) The conversion coefficient values quoted for this line are taken from F. K. McGowan: 
Phys. Rev., 93, 163 (1954) for Z= 61. 
(b) Quoted from ref. (*). 
(c) The relative intensities of these lines were estimated from the coincidence data. 
(d) The unconverted intensities are normalized with respect to the intensity of the 344 keV y-ray. 
(e) The assignments shown inside parentheses are to be considered as doubtful. | 


(5) M. E. Rose, G. H. GOoERTZEL and C. L. Perry: ORNL-report 1023 (1951). 
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4. — Level scheme. 


On the basis of the results of the present investigation a decay scheme 
is proposed for the 13 year Ku (Fig. 10). The intensities of the various 
radiations together with the calculated transition probabilities are listed in 
Table VII. The proposed decay scheme explains satisfactorily the origin of 
most of the observed y-rays as well as the various coincidence data. 


Radiation | Abundance | log fot | log f,t 
Bi = 710+20 keV i” 13 | 10.7 
By = 1040+25 keV | | Di 
Bs = 1450+25 keV | 3 | | 12.3 
EC to the level at 1417 keV | 29 EZIO) 
) 1 248 keV | 20 | 8.7 
) 1191 keV | 3 | 9.8 
» 1090 keV 29 9.1 
» 807 keV 1 11.0 
| 
(*) Caleulated per 100 disintegrations of **°Eu. | 
| (+) Calculated from the nomograph of MoszKowsKI: Phys. Rev., 82, 35 (1951). 
(°) Calculated on the basis of maximum decay energy corresponding to K-binding energy 
of Sm i.e. 46.8 keV. 


152Sm-levels: From the earlier works (**) and also from Coulomb exci- 
tation experiments (8), it has been established that the 121 keV transition 
arises due to the de-excitation of the first excited state in 15?%Sm, and it is of 
#2 character. The rotational character of this level is also evidenced from 
the half-life determination (16) of this state which shows an enhanced £2 trans- 
ition probability as compared to the single particle estimate. The second 
state at 365 keV de-excites by 244 keV y-ray, the measured Æ-conversion co- 
efficient and the K/Z ratio establish it as one of £2 type. The level at 870 keV 
de-excites by a weak y-ray of energy 442 keV, observed both as a weak con- 
version line and also in the coincidence spectra with 121 and 244 keV y-rays. 

The levels at 1090, 1191 and 1248 keV are proposed on the basis of the 
y-y coincidence data, the energy fits the internal conversion lines and the 
observes intensities of the different y-rays involved. The 1090 keV level 
de-excites by a y-ray of 1090 keV going directly to the ground state and 
970 keV y-ray terminating at the 121 keV level; observed conversion coefficients 


(16) A. W. SunyarR: Phys. Rev., 98, 653 (1955). 


517 


518 S. K. BHATTACHERJEE, T. D. NAINAN, S. RAMAN and B. SAHAI 


characterize them as predominantly Æ2 types. The 1248 keV level also de- 
excites by two y-rays; one of energy 1127 keV terminates at the 121 keV 
level and the other of 883 keV ends up at the 365 keV level. The 883 keV 
y-intensity is low and no observable internal conversion line corresponding to 
this y-ray has been observed; the 1127 keV y-ray is mostly £2 in character. 
The observed intensities of the three radiations of 244, 970 and 1127 keV 
terminating at the 121 keV level account very well (Table VI) for the ob- 
served intensity of the 121 keV y-ray after correcting for internal conversion of 
all the y-rays. This leads us to conclude that the direct electron capture to 
this level is negligible. Similar arguments hold for the 365 keV level which 
de-excites only by the 244 keV y-ray whose total intensity accounts rather well 
for the observed intensities of the 725 keV y-ray (to be identified with the 
unresolved 710 keV y-peak as observed in the singles spectrum, Fig. 4, and 
the broad conversion line of weak intensity corresponding to a y-energy of 
= 703 keV) and the 883 keV y-ray, both terminating at this level. It is also 
seen that no appreciable electron capture transition takes place directly to 
this level. 

The 1191 keV level is postultaed in order to explain the existence of some 
weaker low-energy y-rays that were observed both in the electron spectrum 
and in the coincidence experiments. It has been pointed out earlier that no 
appreciable K X-ray or 121 keV coincidences were observed with the 1417 keV 
y-ray. The observed intensity of the 1417 keV y-ray is about the same as 
due to that of 121 keV correcting for conversion and hence, according to 
the arguments put forward before, the 1417 keV y-ray cannot terminate 
at the 121 keV level. The intensity consideration and the coincidence results 
force us to conclude that the 1417 keV y-ray arises due to the de-excitation 
of the uppermost level at 1417 keV and it is predominantly fed by Z capture. 
This shows that the total electron capture decay energy of !5°Eu is more than 
1417 keV and less than 1464 keV. 


152Gd-levels: The Kurie analysis of the £-spectrum shows three groups at 
maximum energies 1450, 1040 and 710keV. The 6-y coincidence experi- 
ments establish the level scheme in %?Gd. No appreciable (-transition takes 
place to the ground state. No 6-y coincidences were observed for any y-ray 
of energy higher than 770 keV (Fig. 9). 

The highest energy f-group of 1450 keV has an «-shape corresponding to 
a unique first forbidden transition of the type AJ = 2, yes, the calculated 
log f,t = 12.3. Since the 344 keV level is a 2+ state, it is most likely that 
the ground state spin and parity of !5*Eu have the value 4 —. This assign- 
ment satisfactorily explains the existence of the 9.2 h isomer of '**Ku, situated 
at = 50 keV above, which emits a 1880 keV 6-group to the ground state of 


Prat 


152Gd with an allowed shape and log ft-value of 7.5. 
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The calculated log f,t = 12.1 for the 1040 keV 8-group also indicates that 
the level at 754 keV has spin /=2+. Because of the Æ1 character of the 
770 keV y-ray as inferred from the measured K-conversion coefficient and in 
conformity with the angular correlation work (17), J — 3 — has been assigned 
to the 1114 keV level. The 710 keV B-group going to this level corresponds 
to a transition of the type AI = 1, no; the calculated log ft-value is too high 
for an allowed transition. The total £-decay energy of 5?Eu is 1794 keV. 


152 

EU (13%) 
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S116, 37 — 
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Fig. 10. — The proposed decay scheme for 13 yr 5?Eu. The transitions shown as dotted 
lines may be regarded as doubtful assignments. All energies are in keV. 


From the proposed decay scheme, the ratio of total electron capture to 


B-emission is given by 


HC e Tai + 1,090 + Tia? 
DI Isa | 


(7) L. Gropzins and H. KenpALL: Bull. Am. Phys. Soc., 1, No. 4, 164B11 (1956). 
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correcting for internal conversion, we get from Table VI 
EC[/P & 4.8. 


The relative intensities and the transition probabilities of the different 
radiations from !5Ru decay are summarized in Table VII. 


5. — Discussion. 


The decay of 5?Eu is unique in the sense that the level structure of 15*Sn, 
as reached by the electron capture decay, and that of !52Gd, reached by B-decay, 
reveal quite different types of nuclear spectra in the region of N = 90 and 
N=88. The large change in the nuclear structure is evident. This sharp 
change in the nuclear structure has been first observed in the study of the 
systematics (18) of even-even nuclei. The nucleus !*Sm with N= 90 shows 
well-developed low-lying rotational states characteristic of the collective exci- 
tations of a nucleus of spheroidal shape. The energy interval rule for such 
rotational levels, corresponding to the ground state rotational band (A = 0) 
is that the first, second and third excited states must have their energies in 
the ratio H,: H,:H, = 3:10:21 and of characters 2+, 4+ and 6+; the ground 
state character of an even-nucleus is always 0+. The levels in Sm at 
121, 365 and 807 keV conform very well to the interval rule and the measured 
K-conversion coefficient of the de-exciting y-rays (pure H2 type) satisfy the 
required spin and parity assignments. On the other hand, in !5Gd, the two 
low-lying levels at 344 and 754 keV are characteristic of the collective vibrations 
of the electric quadrupole type about a spherical equilibrium shape (5). The 
shape oscillations of a spherical nucleus arise due to an oscillating electric 
multipole moment and are classified according to their multipole order, the 
lowest frequencies of collective vibration are in most cases expected to be of 
quadrupole type. Theory (5) predicts that in the limit of harmonic oscillator 
approximation the energies of the first and second vibrational states must 
be in the ratio 1:2, the spin of the first vibrational state is always 1 — 2 +; 
there lies a degenerate triplet of character 0+, 2+, 4+ at the second vibra- 
tional level. In !52Gd we indeed observe that the energy ratio of these two 
states = 2.1; 1—2+ character of the 344keV level is established by the 
K-conversion coefficient of the corresponding y-ray, and 2+ assignment to 
the 754 keV level is also in agreement with the prediction of the theory. In 
the decay of vibrational states, M1 radiation is forbidden even when A1 = 0; 
the 410 keV y-ray corresponds to a 2+ + 24 transition and the estimated 
K-conversion coefficient indicates an H2 character. 


(18) G. ScHARFF-GOLDHABER and J. WENESER: Phys. Rev., 98, 212 (1955). 
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The level at 1114 keV in !5Gd has been assigned a spin /=3. Re- 
cently, relatively low-lying odd-parity states of character 1 —, 3 —, 5 — have 
been observed in a few even-even nuclei in the heavy elements (1%). In the 
9.2 h isomer decay of !5°Eu a state (#2) in 428m at 960 keV has been iden- 
tified as J — 1 — and Æ—0. Such odd-parity collective vibrational states are 
recently (5) interpreted as due to octupole vibrations of even-even nuclei, 
where the one phonon excitation has the character J = 3 —. 

The Unified Nuclear Model (5) predicts for spheroidal nuclei two vibra- 
tional bands for even-even nuclei, each having one quantum of vibrational 
excitation. The so-called f-vibrational band (K= 0) will have rotational 
band members (0, 0+), (0, 2+), (0,4+).... The y-vibrational band (K = 2) 
will have members (2, 2+), (2,3+), (2,4+).... The rotational moment of 
inertia for these upper bands is about the same as that of the ground state 
band. The vibrational quantum energy may be of the order of 1 MeV for the 
heavy elements. The de-excitation of these levels to the members of the 
ground state rotational band will mainly proceed by £2 transitions even though 
M1 transitions are allowed. The intensity rules for the relative transition 
probabilities to the various members of a nuclear rotational band has been 
given by ALAGA (21) and others. 

The respective ratios of the electric quadrupole transitions from the 
1090 and 1248 keV levels in 28m to the states of the ground state rotational 
band provide means of evaluating the K quantum numbers for these 
levels. After dividing out the fifth-power energy dependence we obtain the 
ratio of the reduced transition probabilities. The branching ratio of the 
1090 keV level to the ground state (1—0+, K=0) and 121 keV (= 24, 
K=0) yields B(E2; 2, 2 — 0, 0): B(E2; 2, 2 >0, 2)=0.6, assuming pure H2 
radiations for both transitions. The ratio calculated from the theory (21) 
is 0.7 thus establishing K=2 for the 1090 keV level. The branching ratio 
of the 1248 keV level to the 2+ and 4+ states of the ground state rota- 
tional band yields B(H2; 2, 3 — 0, 2): B(H2; 2, 3 — 0, 4) = 2.1, the theoretically 
calculated ratio is 2.5; the agreement seems to indicate that the two y-rays 
involved are probably predominantly of #2 character and the 1248 keV level 
has also K=2. Thus we may conclude that the 1090 keV 2+ state and 
1248 keV 3+ state are quadrupole vibrational states corresponding to K=2 
(y-vibrations). 

The 1191 keV level de-excites by the 101 keV transition to the 1090 keV 
level. The measured Æ-conversion coefficient of this y-ray is ~ 1.32 and 


(19) F. S. SrePHENS, F. Asaro and I. PERLMAN: Phys. Rev., 96, 1568 (1954). 

(2) L. Gropzins: Bull. Am. Phys. Soc., 1, No. 7, 329J7 (1956). 

(21) G. Anaca, K. ALDER, A. BoHR and B. R. MortELSON: Kgl. Danske Videnskab. 
Selskab. Mat. fys. Medd., 29, No. 9 (1955). 
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K/L & 2.3; this y-ray is predominantly of #2 type since for the M1 transition 
K/L & 6 

value from 0 to 4 and even parity. No spin assignments could, however, 
be made on the basis of our experimental results. 

The assignment 1 = 1 — as the spin and parity of the 1417 keV level in: 
152Sm is consistent with the H1 nature of the 1417 keV y-ray. This is the 
only odd-parity state discernible in ’Sm. Theory (5) predicts that in an 
even-even nucleus the coupling between the quadrupole deformation and the 
octupole mode may give rise to a lowest odd-parity excitation /=1—. The 
only difficulty that one encounters in such an interpretation is that one expects’ 
a competing H1 transition of energy 1296 keV of relatively observable intensity 
to the 121 keV 2+ state. Because of this ambiguity no A quantum number 
has been assigned to this level. In order to ascertain the K-value of the ?’Eu 
ground state, the @-transition probabilities to the different levels in ?Gd have 
to be known rather precisely; the experimental errors in the 6-branching ratios 
are rather high (Table II) and this prevents us from assigning any Æ-value 
to !5:Eu. One fails to notice the X-forbiddenness of some electron capture and 
most of the f}-transitions in this decay. The electron capture transition to 
I=4+ in 18m is a first-forbidden transition (AJ = 0, yes), the higher order 
forbiddenness as revealed by the corresponding log ft-value is most likely due 
to K-forbiddenness. The same situation exists for all the three L8-transitions: 

Detailed studies of this isotope with a high resolution magnetic spectro- 
meter and the experimental determination of the spin of Ku would be of 
great interest. 
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Note added in the proof. 


After this paper was sent for publication, the detailed work of Cork et al. (Phys. 
v., 107, 1621 (1957)) has been published; our results are, in essence, in agreement 
with theirs. Further, the ground state spin of '?EKu has recently been measured to 
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This indicates that the level 1191 keV may have any spin- 
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be T=3. (M. ABRAHAM et al.: Phys. Rev., 108, 58 (1957)) by paramagnetic resonance 
method. This fact implies that the highest energy @-group of 1450 keV terminating 


at 2+ state at 344 keV in 1°Gd is most probably a transition of the type AI=1, yes, 
and it exhibits a forbidden spectral shape. 


RES SON Ox (ES) 


Le radiazioni del ?Eu (4 —13 anni) sono state studiate usando uno spettrometro 
a immagine intermedia, uno spettrometro a lente magnetica sottile e uno spettrometro 
a coincidenze a scintillazione. Si constatano tre gruppi f terminanti a 1450-25, 
1040+20, 710 +20 keV; il gruppo di energia massima ha forma «. Opportune ricerche, 
mostrarono la presenza delle seguenti radiazioni y: 101 (5), 121 (88), 244 (30), 315 (4), 
344 (100), 410 (9), 442 (4), 710 (13), 770 (77), 888 (14), 970 (89), 1190 (96), 1027 (96) 
e 1410 (135) kV. I coefficienti di conversione AK sono stati stimati per la maggior 
parte di queste transizioni e sono stati assegnati gli ordini dei multipoli. Esperienze di 
coincidenza raggi X-y, y-y e B-y ci permettono di costruire lo schema di decadimento 
del ?Eu. La differenza nella struttura dei livelli del 152Sm (N=90) e 152Gd (N=88) 
risulta evidente dai risultati sperimentali. I primi tre stati eccitati del 152Sm a 121, 
365 e 807 keV si identificano come le tre parti della banda rotazionale dello stato fon- 
damentale (K=0); i livelli a 1090 e 1248 keV con I=24+, e 3+ rispettivamente, si 
dimostrano essere gli stati vibrazionali y di quadrupolo corrispondenti a K=2. Lo 
stato eccitato più elevato del !52Sm si trova a 1417 keV ed è probabilmente alimentato 
completamente per cattura L; è questo l’unico stato con parità dispari del #?$m avente 
carattere /=1 —. I livelli del 152Gd a 344 e 754 keV, ambi con I=2+, raggiunti col 
decadimento 8 sono caratteristici delle vibrazioni collettive di tipo quadrupolo elet- 
trico intorno a una configurazione d’equilibrio sferica. Il livello a 1114 keV del !52Gd 
ha caratteristica 1—3 — e si ritiene dovuto alle vibrazioni da ottupolo. Allo stato fon- 
damentale del 2Eu è stato assegnato lo spin /=4 —; si sono calcolati i rapporti e le 
probabilità di transizione per i gruppi Be le transizioni per cattura elettronica. I risul- 
tati sperimentali si accordano con le previsioni del modello nucleare unificato. 


(*) Traduzione a cura della Redazione. 
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on the Spin Kinematics of Electrons in a Uniform Magnetic Field. 
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(ricevuto il 20 Settembre 1957) 


Summary. — The influence of a uniform magnetic field on the spin orien- 
tation of electrons, or other spin } particles, in a beam is computed by 
assuming that the Dirac-Pauli equation correctly represents the electrons 
having anomalous magnetic moment. It is shown that, when the field 
is orthogonal to the electron beam, the anomalous magnetic moment 
causes the spin to turn about the field direction faster than in the normal 
condition. This effect depends on the energy of the electrons and becomes 
particularly large for energies which are high compared with the rest 
mass: for about 250 MeV energies a half circle path in the field would 
suffice to convert the longitudinal polarization state of the beam into 
the transverse one. The calculations are performed by making use of 
the exact solution of the Dirac-Pauli equation. 


1. — It is interesting for many purposes to evaluate the effect of the ano- 
malous magnetic moment on the spin kinematics of electrons, or of other Dirac 
particles, moving in a uniform magnetic field. This effect is at present of 
special interest due to the fact that polarized particles can now be injected 
directly into an accelerating machine. 

MENDLOWITZ and CASE (1) have recently derived the angular velocity of 
the spin precession of an electron having anomalous magnetic moment. Their 
starting point is the Dirac Hamiltonian with the additional term of Pauli 
type aufo-H; subsequently, however, they have analyzed the problem in 


(1) H. MenpLowITZ and K. M. Case: Phys. Rev., 97, 33 (1955). 
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the FOLDY-WOUTHUYSEN (?) representation obtaining the new Hamiltonian up 
to terms of the first order in w,H. By means of this method, the authors have 
found that, when the constant homogeneous magnetic field is perpendicular 


to the electron beam, 
a 
(el Ds = W, (1 va #5 ’ 


where ©, is the relativistic cyclotron frequency, a = «/2z is the second order 
radiative correction to the spin magnetic moment, f = v/c. The result con- 
tained in the (*) deserves to be verified not only because of the importance 
it may have in the experiment of LOUISELL, CRANE and Prpp (?) for the 
direct measurement of the magnitude of the gyromagnetic ratio of the electron, 
but also for the large effect that it might produce when high energy particles 
are considered. In fact, for an electron energy of about 50 MeV, one should 
expect a 10 percent effect on the angular velocity; in other words, a few 
revolutions of the electron beam path should, in contrast to what is generally 
believed, suffice to obtain an observable anomalous spin rotation. Of course 
this statement is correct under the hypothesis that particles with anomalous 
magnetic moment are properly described by the Dirac equations in conjunction 
with the Pauli term awfo-H. This assumption, of course, leaves open the 
question whether it is correct to interpret the anomalous magnetic moment 
as an intrinsic property of the particle or such moment should be attributed 
to radiative correction; in the latter case the Pauli term may not suffice to 
describe the behaviour of the particle, whatever its energy. Bearing this limi- 
tation in mind, the spin rotation has been directly calculated in the present 
work by using the exact solutions of the Dirac-Pauli (‘) equation in a uniform 
magnetic field. In Sect. 2 it is shown that the generally used definition (5) 
of «spin direction » is again convenient using the Dirac-Pauli equation. 

The solutions are derived in Sect. 3 according to the method used in a 
preceding work (6). In Sect. 4 the spin direction for the electron beam emerging 
from a magnetic field is explicitly obtained as a function of the spin direction 
of the incident beam. The result agrees with that obtained by MENDLOWITZ 
and CASE by using the Foldy-Wouthuysen representation. 


(2) L. L. Forpy and S. A. WouTHUYSEN: Phys. Rev., 78, 29 (1950) see also K. M. 
Case: Phys. Rev., 95, 1323 (1954); furthermore H. MENDLOWITZ, in his thesis (unpub- 
lished), has treated, in the Foldy-Wouthuysen representation, an electron with an 
anomalous moment in a constant magnetic field, but we have not seen this work. 
(3) W. H. LurseLL, R. W. Pipp and H. R. CRANE: Phys. Rev. 94, 7 (1954). 
(4) W. Pauti: Hand. der Phys., 24, 232; for a discussion of the Dirac-Pauli equation 
in a uniform magnetic field see also: P. CALDIROLA: Rend. Ist. Lomb., 77, 1 (1944). 
(5) See for example: N. F. Morr and H. 8. W. Massey: The Theory of Atomic Col- 
lisions (Oxford, 1949). 
(6) M. CARRASSI: Nuovo Cimento, 5, 955 (1957). Referred to as (I) in the text. 
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2. — Let us consider the relativistie equation of Dirac-Pauli: 


(1) (3 Pr. + ime +23 3 oP, JF=0, 


uv 


were — 1,2} 3,4; do = (iPa, 5), a and f being the usual Dirac matrices; 
0,,= (8,8, — B, Fa =(0A4,/00, — 0A,/0x,) with v= (r, ict) and A,= (A, id) 


7, = — th (0/0a,) + aN À is a constant of the dimensions of a magnetic 
moment so that the magnetic moment of the electron is w= u,+ À (u, is the 
Bohr magneton). If we express the electromagnetic tensor #,, in terms of 


the field intensities and write ~ = y exp [— (i/h)wt], the equation (1) becomes: 
(2) {w+ eD + can + Pm c? + ABlo-H + i(a-E)} py = 0, 


where y is now the spacial wave function with four components: y, (j =1, 2, 3, 4) 
and w= (p’c?+ mie). 

Since our object is to study the spin kinematics in a magnetic field, the 
Spin direction is taken to be defined, as we did in (I), in the usual way through 
the ratio: 

(3) Ps _ te exp [io], 

Ya 2 
where y and © are the polar angles of the «spin direction » in the co-ordinate 
system in which the electron is at rest. 

This definition is the most convenient for all the problems concerning 
polarization states of the electrons, as TOLHOEK and DE GROOT (’) have clearly 
shown; furthermore, for a fast electron it is also the only method available 
at present by which the spin might be observed in practice since, as we already 
know (5), a decelerating electrostatic field does not change the ratio y,/y,4. 
This last property, indeed, might not be satisfied if we use equation (2) to 
describe the motion of the electron with anomalous magnetic moment. Actually, 
equation (2) contains a new term of interaction depending on the electric field; 
however, this term ijfa:-E, effects a mixture of the components in the same 
way as the term ca x and hence causes no difficulty. Indeed if we confine 
ourselves to the case of an electron moving parallel to the 2 axis, there being 
also an electrostatic field in this direction, we cally show that y, and y, both 
satisfy the same differential equation. 


Ò À 0D 
(4) (w + eD — mec?)p + hc Gira) 
E CA À 0D 


(?) H. A. ToLHoEK and S. R. DE Groot: Physica, 17, 17 (1951); H. A. TOLHOEK: 
Rev. Mod. Phys., 28, 277 (1956). 
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Thus, since the equation (4) is homogeneous and y, and Ya both satisfy 
the analogous boundary conditions, it follows that the ratio W3/Y, does not 
change. The same result can be deduced calculating explicitly the solutions 
of equation (2) as ToLHOEK and DE GROOT have done for a Dirac particle (8). 

This remark is also of interest since even when considering electrons with 
anomalous magnetic moment, any longitudinal electric field does not change 
the spin direction. 


3. — If we consider the influence of the magnetic field alone, the equation (2) 
becomes: 


(5) Star + img +5 BoHly =0. 


To solve the equation (5), we confine 
ourselves to the following case: H=0 for 
æ< 0 and H=H,=const. for «> 0 (see the 
figure). 

We proceed in a way analogous to that 
followed in I) in order to obtain the phys- 
ically interesting solutions. Hence, putting 
A, =— Hy, A =the, A, —0, the set. of 
equations (5) can be reduced to the system ;<0 
of two second order differential equations 


2 € LE 
(ai - ni) + ROH — (1 + 20020 pit 


+20, (202+ im)ay = 0, 


eee iti a, AH 
0 + moe? + H ? arie nee 


(0,1 


Among the four unknown functions there exist the relations: 


C; . 
| Ta [oz “a eV 
(7) | cà A i 
| a L(Te + dt) Va — Ta] . 


(8) E. T. WHITtTAKER and G. N. Watson: Modern Analysis (Cambridge, 1952). 
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Since 7, commutes with 7, and z,, it commutes also with the hamiltonian 
and therefore p, is a constant of motion. Then we put p,=0 and so we 
treat the problem as independent of the ¢-coordinate. This assumption means 
that the propagation vector of the plane wave in the region 4 < 0 is perpen- 
dicular to the z-axis. The system (6) is reduced to two independent diffe- 
rential equations 


(6) SEO moe — (Ta + Ty) — (i 2H sù | o. (|) UE 


e e 


while equations (7) become 


| Vi “Hi (yg — Uy) Va, 
(0e) = 
C; 1 
| dicon (x + Wty) Ps 
System (6'), putting w = eH/2c, can be written 
w? + 22H? n ooh = Ò ; Ci 
(8) Ti mg + REV? — 2ilio ly È 2x) 
— w? (a? + y?) — 2 (10 ale | o. (i = 0° 
\ € Pa 


Using the same substitution as in (6,1) and also É=2(r+(7/20)) Volt, 
we obtain 


; d? Il eH? wa 
(8') fe Dr Le | 7 hoë? [ro à | on 


Putting finally 


ia i MH? 2w4H 
| tent pe? pec ? 
(9) 
F 1 VA? | 2waAH 
FEV, pee pee ’ 


where » = p?/4fiw is the number of De Broglie wave lengths contained in a 
half circle of radius r= cp/eH, we can write the two equations 


| dos ne be cd» 
1 O e ie 
(10) | d'y, Papi 1 9 oA 

rae if) | 


We observe that (10) differs from equations (7,1) obtained using the normal 
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Dirac equation, only for the values of the constants y’ and »’. Therefore, all 
terms of interaction with the anomalous magnetic moment are contained in 
these constants. The solutions are again the parabolic cylinder functions and 
we put g, = AD, _,(&) and y, =BD,.(éÉ). Using (7'), (6,1) and the known re- 
currence formulae of the D, functions, we find the solution of the equation (5) 
in the case «> 0 


AV ho S 
(wje) + me + ANTO 0 Pra 
= AV hw An. i 
(11) y=| Co + me— UV exp |} (024 +n). 
ADy_, 
BD, 


We now have to match, in x = 0, the solution (11) with that in the region 
“<0 (see (9,I)). From this condition the following system is obtained 


(By — B;)p, + (By + B.)p,= y: WV ho BY Dy (6) 


a) (A,— A; )Pe — i(Ay + Asp, = — 2iy_Vho AD, (Eo) 
A, un A, =" AD, 2,(E) 
B; ala Bb, En BD,.(£0) ’ 
where 
Seve a EN ie) aN ee ak 
“muri vali) 4 ai w+ moe) ~ 


From (12) we can easily derive the ratio B,/A,, which specifies the pola- 
rization state of the reflected beam in terms of the initial data for the incident 
plane wave. 

Let « be the angle between the momentum of the reflected beam and the 
momentum of the incident beam. Then, as shown in the figure, «=a — 2%, 
where o, = 627 p,/Dz.- 

Hence we have 


(13) exp [ta] = — ne 
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Using (13), from system (12) we deduce 
jt NT 
14624 Lr) 1e 2, 2x) 
Px Pu V+ v o Pa Pe v 
~~ = à 
fe PT D 14i( 2, Pr) 
Da  PeY+1 ¥ Pa Pa y 


where we have used the definition v = p?/4hm and we have put 


TES 16%, 3 
14 — = — exp |ta 
(14) ven 6 p [ta] 


Dy(§o) R- D,(&o) 


Tape Re ee = == . 
Vy" D,'_1(Ë0) Vy! Dy _1(&o) 


Now, putting 


14/» 3 ly 
a nd pee? VEE: 


Pa Pa V+ v Px Px ‘a 
we have 
li “ea RES È 
1+iM 1—in 14 MN RI DI 
D EN MN Se He 
1+MN 9 
Therefore (14) can be written in the useful form 
JE BN, 2 
(15) ani [i(a + d)], 
where 6 is defined by the relation: 
ly 13/v 
n z + 
DEN rel yo 15») 
(16) tg = = = = = : 
2 Ji oa MN y! ibs lay y v' 
p—p,\y_ || —R - -R+]+ p al _ R+R- 
y V+ v V+ Vv 


4. — It may be noted that (15) and (16) are exact formulae because till now 
no approximations have been made though we have considered particular 
conditions. However it will be convenient now to remember that v=108--1015, 
in every case of interest; hence »’ and »” according to (9), are of the same order 
of magnitude. In addition we may observe that putting 1 — au, = a(eh/2m, ¢), 
where a is the percentage correction to the magnetic moment, we can write 


RTE ae PTT UE CPE A NT Mes EU a AU, g' : re, NI at PE 
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and therefore (9) gives 


1 2 2 
| is Ep ae (È) (25) 
Vv 2v MC? 4v) \m,c? 
(17) UÀ 2 2 
EEE) 
v 2y M,C? Ay} \moc? 


Furthermore if a is of order 107%, it is clear that the last term in the two 
relations (17) is quite negligible both compared to unity and to the first order 
term. With this in mind and remembering the definitions of y: and y_, we 
obtain 


With these simplifications, (16) becomes 


Ô (p./p)(R” — R*) 
di 5-1 EP +R) RFR 


The functions R* and À can be evaluated asymptotically for large values 
of the parameters »’ and »”, as we have shown in the Appendix. Hence, we 
obtain: 


& Dy($0) cos (y' — Bo) 
REN = i. = ctgy' cos fo + sin f 
anne) at: È si 
(19) n 
i a cf, e 
CET NT 
where 
7 Il 1 
È TV > Pu a 207 | 1 py ES 9 Py | q si | Py 
On RE A im)! 2p’? 2 p\' ame) 7 2p 


By putting (19) into (18) and keeping terms up to the second order in 
lp = Bo, we find: 


Ô cotg y° — cotg y” 
2 14 cotgy'cotgy” 


tg 


I ! Py W 
tte li T4 | 


D mel 


= (v° y’) 


From the expressions of »’ and »” and assuming p,/p=f,=% we obtain 


w 
(20) 0 — à aia (7% — 2a) . 
0 
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Hence from (15) finally we have 


Ul 
J 


ilo | a(14 a = 


Me 


It appears from this result that the spin direction of the reflected beam is 
rotated around the direction of the magnetic field by an angle larger than 
the rotation angle of the momentum. If «=z, one gets 


B, Dine ae w 
ARSA Si moe | * 


The last result, corresponding to the particular condition p,= 0, can be 
obtained directly, as we have shown in the Appendix, without using the 
asymptotic formula for the D, function. In this particular case the effect of 
the anomalous magnetic moment on the spin precession about the field di- 
rection might become particularly large for sufficiently large energies. Indeed, 
the rotation of the spin exceeds the «normal precession», by an angle 
Aa = ra(w/moc?) and therefore for an energy w= 250 MeV, one gets 
Ac = 7/2. Thus, at this energy, an incident longitudinal polarized beam is 
«reflected » by the field with a transversal polarization. 

The angular velocity of the spin precession can be easily obtained from (21); 
in fact, if r= cp/eH is the radius of the classical half circle described by the 
electrons in the magnetic field, the time spent in traveling the path ar is: 


alt Mo x Mo € 


t= - = : 
Vip vip eH 


Hence the angular velocity is 


A _ «(1 + a(w/moc?)) VIF eH (1 fe RR eee eee 
pag t ie ITA pe?) DEN 4/1 83) 


where w, is the angular velocity of the normal spin precession and therefore 
our result agrees with that obtained by MEDLOWITZ and CASE. 


The author wishes to thank Professor A. BORSELLINO for encouragement 
and helpful advice received in working out this paper. 
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APPENDIX 


For the calculation of the R*+ and R functions, it is necessary to find 
explicit expressions for the D,,(£,) and D,:(é,) functions for large values of the 
parameters v and v”. In the mathematical literature (*) there are, under 
suitable hypotheses, several asymptotic expressions for the confluent hyper- 
geometrie functions for large value of the parameters on which these depend. 
In our case, the question is more complicated, because the variable 
E, = p,|Vhw = (p,/p)2V% contains also the magnitude of the parameter itself. 
We refer to a detailed work of N. SCHWID (?°), where two independent solu- 
tions of the Weber’s equation are obtained assuming n to be very large. These 
solutions w,(z) and w,(z) are related to the usual parabolic cylinder functions 
as follows: 


| P(A) green (€) = Lt + DD (ie) + Dati, 
(LA) | 
| P(G +2) Berea) = illa + DID) — Dio) 


(n is real). Further, by using the relation: 


(2A) D,(z) = 
= leo LD (iz) + D_,-,(—iz)] + isin = DE Nat i, 


the expressions of w,(z) and w,(z) and the asymptotic expressions of the /'(n) 
functions, from (1A) and (2A) one obtains: 


e-"l2prl20# J PRA (saree Tei | 
(3A) DO) peo Oh a Dj 0s | -(n 3) | + O[n 1, 


The relation (3A) holds good within the semicircle of radius |2|— 
—[2(2n +1)]} contained in the half plane R(z) > 0 excluding the neighbourhood 
of the circumference. For the purpose of our calculation, we must consider 
that <= è and that n will be either equal to v' or y'— 1 for the functions in- 
volved in Rt, or, alternatively, to v" or v’—1 for the functions ni in: ie: 
In order to ensure the validity of (3A) and noting that §,= p)2v/v, we 
shall assume, although this hypothesis is not strictly necessary, joe Pulp LA 
(i.e. we shall keep terms up to the second order in p,/p —sinx). This hypo- 
thesis only means that we limit ourselves to consider those particular cases 


(9) A. ERDELYI, W. Magnus, F. OBERHETTINGER and F. TrICOMI: Higher Tras- 
cendental Functions, vol. I and II (New York, 1953); W. Magnus and F. OBERHET- 
TINGER: Special Functions of Mathematical Physics (New York, 1949). 

(10) N. Scuwip: Trans. American Math. Soc., 37, 339 (1935). 
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in which the angle a is not wider than a few degrees; furthermore, it allows 
considerable simplifications in the final formulae. Particular care shall be 
given to the simplifications of the arguments of the trigonometric functions 
contained in (3A). In fact it is permitted to neglect the terms of the order 
of 1/v’ (or 1/v"), but unity is not negligible compared to »v' (or v"). Making 
use of (17), one obtains: 


(4A) 


From (4A) and (3A) one thus finds: 


(1 — (&/49')) 2 

1/2 és PL (y Le 1) ; ay’ 3 Lo a w | all Dy 
(1 — (8/47)? sin| ? ( | 

4/2 et a a ay" Dy a w i Dy 

Bett) = ane 2G Ot Sia) a 

1/2 ev D2 (y — 1) 
(1 — (£5/4v°))? 


1/2 e v'12y!v'12 TV Dp, a w 1 i 
Dio). = 1 COS 2 7 dea z| 


Dy_1($0) D 


(5A) 


Dy (So) = sin 


Since for large values of n, (1+(A/n))"= e‘, we may put in (5A) 


L 
2 


ES poll E 
(pyr prete) 


x 
; till 2 NE 
and (y" — 1)" v2 = pe |). 
Vv Vv 


Using these together with the (5A) and the definitions of Rt and R , we at 
once get (19). 
Whenep,= 0; p—#»., (18), becomes 


6h 0) =O) 
Si RO) RON 


(18") tg 


Considering that 


Tan 
D,.(0) = T(A— m)/2) ’ 


534 


THE INFLUENCE OF THE ANOMALOUS MAGNETIC MOMENT ON THE SPIN ECC. 535 


and using the well known relations satisfied by the Euler functions as well 


as the asymptotic expression 


T(n) = en t4/2x , 


while bearing in mind the definitions of R*, R-, one immediately finds 


" ! 


R*(0) = cotg = , Lo cotg 


from which (20) is directly obtained in the particular case in which «=a 
(o == 0). 


RIASSUNTO 


Viene calcolato l’effetto di un campo magnetico uniforme sulla orientazione dello 


spin di un fascio di elettroni, o di altre particelle di spin 3, supponendo che il moto 
degli elettroni, dotati di momento magnetico anomalo, sia correttamente descritto 
dall’equazione di Dirac-Pauli. Si dimostra così che, quando il fascio di elettroni è nor- 


male alla direzione del campo, l’effetto dovuto al momento magnetico anomalo è di 
far precedere lo spin intorno alla direzione del campo più velocemente che nel caso 
normale. Tale effetto dipende dall’energia degli elettroni e diviene particolarmente 
vistoso per energie elevate rispetto a quella di riposo: per energie di circa 250 MeV, 
basta il percorso di una semicirconferenza nel campo per mutare lo stato di polarizza- 
zione del fascio da longitudinale in traversale. I calcoli sono stati eseguiti facendo uso 


delle soluzioni esatte dell’equazione di Dirac-Pauli. 
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An Experimental Test of the Fermi-Teller « Z-Law » (*)(*). 


J. C. SENS, R. A. SwANSON (*), V. L. TELEGDI 
and D. D. YOVANOVITCH 


Enrico Fermi Institute for Nuclear Studies, University of Chicago 
Chicago, Ill. 


(ricevuto il 2 Ottobre 1957) 


Summary. — Frrmi and TELLER predicted, in their classic paper on the 
slowing-down process of negative mesons, the relative probabilities with 
which the latter should be captured by the different atoms of a chemical 
compound. They concluded that these atomic capture probabilities should 
be proportional to the nuclear charges, Z;, of the atoms considered. The 
experiments here described were carried out to check this prediction. 
They are based on the fact that negative muons reach rapidly the K-orbits 
around the nuclei of the atoms of the material in which they are slowed 
down, and disappear from these orbits (either by decay or absorption) 
at rates which are characteristic for the nuclei in question. Thus the 
time distribution of decay electrons emerging from a chemical compound 
in which muons are brought to rest presents in general as many expo- 
nential components as there are constituents; their intercepts at time zero 
(entry of the muons into the material) are directly proportional to n,p,, 
where n;=fraction of component of atomic number Z;, p;= capture prob- 
ability into the mesic K-shell of atom of same constituent. Neglecting 
the possibility of transfer of muons between neighbouring atoms, these 
p; are equivalent to the capture probabilities predicted by Fermi and 
Teller. The present experiments show that at least in insulators (A1,O,, 
P,0;, Si0,, CCL, KOH, KHFE,, p-dichlorobenzene) the predictions of 
Fermi and Teller do not hold. No weighting of n, with Z; seems to arise, 
the lighter element being even preferred in some cases. As an explanation, 
it is pointed out that the energy loss to lattice vibrations, not considered 
by Fermi and Teller, may be an important mechanism in insulators for 
negative particles of nearly zero kinetic energy, i.e. at the point of being 
captured. 
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(*) A preliminary account of this work has been presented at the 1957 Denver 
Meeting of the APS (Bull. Am. Phys. Soc., II, 2, No. 6, 320 (1957)). 

(+) This research was supported in part by the U.S. Atomic Energy Commission 
and the Office of Naval Research. 

(*) National Science Foundation Fellow, 1956-57. 
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AN EXPERIMENTAL TEST OF THE FERMI-TELLER « Z-LAW » Doi 


1. - Introduction. 


In the final section of their classic paper (1) on the capture of negative 
mesons in matter, FERMI and TELLER proposed to predict the relative proba- 
bilities with which mesons are captured by the different kinds of atoms when 
the slowing-down takes place in a chemical compound. On the basis of some 
estimates which they themselves considered crude, these authors concluded 
that the atomic capture probabilities should be proportional to the respective 
nuclear charges, Z;, of the atoms in question. 

In many experiments, negative mesons are necessarily brought to rest in 
chemical compounds and mixtures thereof (e.g. in nuclear emulsions, bubble 
chamber liquids, radiochemical targets, etc.). The interpretation of such ex- 
periments depends often critically on the relative capture probabilities referred 
to above. Lacking direct experimental evidence on this point, many inves- 
tigators have in the past relied on the theoretical conjecture of FERMI and 
TELLER, which has then become known as the « Z-law ». 

No experiments specifically designed to test the validity of this «law » in 
a variety of substances appear to have been carried out so far, while sundry 
contradictory evidence has been obtained as a by-product of various experi- 
ments involving negative mesons (24). Among these, the work of STEARNS 
and STEARNS (4) is most significant and shall be discussed in detail later. 

The present investigation, made possible by the availability of high in- 
tensity, high purity u--beams at the Chicago synchrocyclotron (5), aims to 
provide reliable empirical data on the relative capture probabilities in a number 
of chemical compounds. It leads to the conclusion that the « Z-law » does 
not hold in practice, the distribution of captures among the various consti- 
tuents being approximated more closely by the simple unweighted atomic 
ratios of the latter in the compound. While this conclusion is based exclusively 
on experiments performed with muons, it is reasonable to anticipate on the 
basis of our general understanding of the slowing-down processes that it will 
apply unchanged to the capture of negative pions. 


2. - Experimental method. 


Negative muons captured by an atom reach the mesic K-orbit in a time 
very short compared to their lifetimes in condensed matter, and disappear 


(1) E. FERMI and E. TELLER: Phys. Rev., 72, 399 (1947). 

(2) W. K. H. PANOFSKy, R. L. AAMODT and J. HADLEY: Phys. Rev., 81, 565 (1951). 
(3) A. FAFARMAN and M. H. SHamos: Phys. Rev., 100, 874 (1955). 

(4) M. B. STEARNS and M. STEARNS: Phys. Rev., 105, 1573 (1957). 

(5) N. P. CAMPBELL and R. A. Swanson: unpublished report. 
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from this orbit at a rate, 4,, that is characteristics of the nuclear charge Z. 
Indeed, A, — A, + A,, where A, is (to an excellent approximation) the muon 
decay rate in vacuo, and /, the rate at which the muon is captured by the 
nucleus Z from the K-orbit. A, is proportional to Z4, or rather Zi, (°). 

Decay electrons will appear at a rate proportional to 4/43, with an ex- 
ponential time dependence characterized by A,, time zero being that of the 
arrival in the Æ-shell. When muons are captured by the different consti- 
tuents Z; of a chemical compound, the time dependence of the electron decay 
rate will be a superposition of various exponentials, each with a different 
A,(Z;). Indicating the relative capture probability into the X-shell of con- 
stituent atom Z, by p,, it is clear that the magnitude of any component of 
the observed electron decay rate at ¢= 0 will be proportional to (4,/A,)Ap;n;, 
i.e. simply to p,n;, where n; is the fraction of atoms Z, in the compound. 
Barring the possibility (7) that a muon be transferred from one atom to another 
while cascading from its first bound orbit around a specific atom to its ultimate 
K-orbit, these p; will represent directly the atomic capture probabilities sought 
by FERMI and TELLER. They will in any case represent the quantities of 
most immediate practical interest, being the relative capture probabilities into 
the orbits from which nuclear capture occurs. 

To implement this method in practice, one has to determine accurately 
the time dependence of electron decays from negative muons stopped in a 
chemical compound. In principle it is not necessary to know a priori the 
A,’s of the various elements contained in the compound, as these parameters 
can also be obtained from the analysis of the electron decay curves. It is 
however statistically advantageous to work with compounds consisting of 
elements which: (a) have been, or can be, investigated in the form of ele- 
mental targets, (b) have reasonably different disappearance rates, i.e. differ 
sufficiently in atomic number. It is obvious that binary compounds yield data 
that can be analysed less laboriously and more reliably than those from others. 

The equipment used to determine the electron decay curves was the one 
repeatedly referred to in reports (*°) on muon work done in this laboratory, 
and need not be described here in detail. It consists of a scintillation counter 
telescope defining the entry of muons onto the target under investigation, a 
scintillation counter telescope defining the exit of decay electrons from said 
target, and of a time-to-pulse height converter circuit fed by the coincidences 
from these two telescopes. The output of the converter is connected to a 
100-channel pulse height analyzer. In the present experiment, the total useful 


JA 

T. B. Day and P. Morrison: Phys. Rev. (in print). 
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range of the converter was 70 channels, corresponding to 5.15 us. The con- 
verter was repeatedly calibrated and checked for linearity during the experi- 
ment by means of a precision variable delay line (*). 

For each of the compounds for which data are reported here, at least 3-104 
electron counts within the useful range indicated above were collected. The 
counts from each successive 5 channels of the analyzer were totalled, and the 
data so obtained were first analysed graphically (+). The results so obtained 
were found to be somewhat subject to the uncertainties connected with the 
method, and in particular sensitive to the assumptions made about the « back- 
ground » counts (these arise experimentally mainly from random coincidences). 
The data were hereafter subjected to least-squares fits to forms of the type 
CO, exp [— A,(1)t] + C exp [— A,(2)t] + B, assuming the A,’s as known and 
weighting the points with the inverses of their variances. This analysis was 
carried out by means of punched cards on an IBM-650 computer. The relevant: 
quantities n,p; were derived from the C,’s by applying corrections (a) for the 
lumping of data into groups of five channels, (b) for the known displacement: 
of the «earliest » group of channels with respect to the true zero of time. In 
all cases the backgrounds B yielded by the least-squares analysis had magni- 
tudes that were reasonable for the counting rates involved, as judged from our 
experience with elemental targets (?). 


3. - Results. 


We indicate underneath information relevant to specific targets and runs. 
The results are summarized in Table I, while Fig. 1 through 5 represent 
standard semilogarithmic plots of the electron decay rates for a number of 
targets. 


(a) SiO,. Target: 6 g/em? of Vitreosil, a commercial form of high purity 
fused quartz. Total electron counts: 3.18-104, B= 6-10? counts/us. 


(b) ALO;. Target: 10 g/em? of refractary brick. Total electron counts: 
14-104, B= 4.6-10° counts/us. 


(c) P,O;. Target: 6 g/em? C.P. grade, anhydrous material, sealed in a 
thin polyethylene bag. Total electron counts: 4.9-104, B= 1.4-10* counts/us. 


(d) CCI. Target: 6 g/cm? of liquid C.P. grade material, contained in 
a thin plastic bag. Total electron counts: 7.3-104, B = 1.4-10? counts/us. 


(*) Manufactured by Advance Electronics Co., Passaic, N.J. 
(+) These results were quoted in our preliminary communication, Bull. Am. Phys. 
Soc., II, 2, No. 6, 320 (1957). 
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TABLE I. — Summary of results on reltaive capture rates. 
Compound Ratio This expt. (*) HET (6) Atomic (*) 
| 
P,0; O/P OTe ss 1.3 2.5 
CEL, C1/C 4.1 4.8 1102 4.0 
ALLO; O/Al 2.3 + .2 .92 15 
SiO, O/Si 35H LI 2.0 
KHF, F/K 1.7 + 4 .95 2.0 
KOH 0/K DOTE A 42 1.0 
cié cl (1) C/C1 2.3 +2 1.06 3.0 
» » (8) C/C1 PARSO 1.06 3.0 
(*) The errors assigned to these ratios are larger than those from statistics alone. Where the 
decay rates A, are known from independent experiments (Al, O, P, Si, C) allowance has been 
made for the uncertainty in these rates; where the decay rates of one or more constituents 
| were derived from the same data, the error has been additionally amplified by an estimated 
| factor. For F, A,=— 6.31:105s 71 has been adopted, contrary to Ref. (°). 
| (+) The ratios indicated here are those conjectured by Fermi and Teller, i.e. correspond 
| to the assumption that p;= Z;. 
(*) The ratios indicated in this column are the stoechiometric ratios in the compounds, 
| i.e. correspond to the assumption that all p;= 1. 
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Fig. 3. — Time distribution of electrons from 
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(e) KOH. Target: 6 g/cm? of pellets 
of C.P. grade material, contained in thin 
polyethylene bag. Total electron counts: 
3.1:10*, B= 7.7-10? counts/us. 


(7)- KHF,. Target: 6 g/cm? of C.P. 
grade material, contained in thin poly- 
ethylene bag. Total electron counts: 
4.6-104, B= 8.7-10? counts/us. 
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contained in a thin pyrex beaker. Total electron counts: 2.6-104, B= 
= 6-10? counts/us. 


(h) p-dichlorobenzene. Target: same as (9), but maintained in liquid 
state at about 70 °C. Total electron counts: 2.6-104, B= 6-10? counts/us. 


4. - Interpretation of results and conclusion. 


It is clear from Table I that the observed capture ratios depart in all cases 
from the Fermi-Teller prediction and follow more closely the simple atomic 
ratios, unweighted by the atomic numbers, indicated in the last column. In 
some cases the constituent of lesser Z appears even to be favored. 

Before venturing into hypotheses about the reasons for this discrepancy, 
two rather obvious facts need consideration. First, not all the compounds 
investigated by us are equally reliable sources of information. Targets (a) 
through (d) consist of truly binary compounds, while the remainder contain 
hydrogen as a third element. The presence of hydrogen is not noticeable in 
the electron decay curves of these ternary compounds, where a component 
of A, essentially equal to A, should appear, if protons were capable of cap- 
turing and retaining a muon in these substances. On the other hand, as anti- 
cipated by FERMI and TELLER (1) and suggested by PANOFSKY’s observations (2), 
‘the small neutral system formed by a hydrogen mesic atom would easily dif- 
fuse through matter and ultimately transfer the muon to a nucleus of higher 
charge. If we extrapolate the conclusion that the constituents of a compound 
~ capture mesons simply in their ratio to hydrogen, then in a compound like 
KOH the ratio of mesons found in the K-shells of K and O could be appreciably 
affected by this transfer mechanism. However, even under the unlikely as- 
sumption that in this compound (as contrasted with the true binary cases) 
the K atoms would have a 2--4 times larger relative capture probability than 
the O atoms, and that the mesic hydrogen atoms would transfer the mesons 
overwhelmingly to oxygen, the observed capture ratio could not be attained. 
It might be that practically no captures occur in hydrogen that is ionically 
bonded, as the protons have practically no electrons «of their own » to con- 
tribute effectively to the last stages of the stopping process (1). 

Second, the case of ALO, needs clarification. This substance has been 
. investigated by STEARNS and STEARNS (*) from the same point of view as in the 
present investigation, but by a different technique. These authors concluded 
that the « Z-law » is well borne out in this substance, in contrast with our 
results on this and other compounds. Their experiment consisted in com- 


(0) We are greatly indebted to Prof. A. TURKEVICH for pointing this out to us. 


nN 


12 


ai ns ee aes LES à - a Seek een. 
ede CARA Ac x À 
a SA paw a 
AN EXPERIMENTAL TEST OF THE FERMI-TELLER « Z-LAW » 543 


paring mesic X-ray spectra (and absolute yields) from two targets: (a) ALO,, 
(b) a macroscopic mixture of Al and O consisting of the same gram-amounts 
as the constituents of the. oxide sample. The Al was in the form of 1 mm thick 
sheets, the oxygen in. the form of water filling the interstices between these 
sheets. The spectra from these samples, involving K-X-rays from O and 
L-X-rays from Al, resulted essentially identical. From this we would conclude 
that the macroscopic law of stopping governing the distribution of captures 
among the adjacent sheets of water and aluminum is essentially the same 
as the microscopic law of stopping governing the distribution of captures among 
the Al and O atoms. The Stearnses assumed that the relevant macroscopic 
stopping powers are proportional to Z, and argued therefrom that the « Z-law » 
holds microscopically. While their assumption is reasonably justified at high 
velocities, the velocity range in which the stopping powers of Al and H,O 
have to be compared is that of a muon that has barely enough energy to 
traverse one of the sheets. A rough estimate of this energy is 4 MeV. O is 
much more effective in stopping than Al at this energy, the ratio of stopping 
powers per gram being (1) about 1.6 and increasing in favor of O as the 
particle velocity decreases. Thus the Stearns’ data on ALO, are in reasonable, 
if not exact, agreement with ours. 

The Stearnses also report that in CaS the « Z-law » is fulfilled to within 
+ 20%. The ratio of atomic numbers being 20/16, this conclusion is perfectly 
compatible with the captures occurring in this compound too in the simple 
atomic ratio. 

: The prediction of FERMI and TELLER was based on the behavior of the 
energy loss of negative mesons near zero kinetic energy, which these authors 
showed to be proportional to Z. In seeking the cause of the discrepancy 
between their prediction and experiment, it is perhaps worth pointing out 
that FERMI and TELLER gave a reliable estimate for the energy loss only for 
the case of metals. They indeed pointed out that in insulators certain diff: 
culties arise at low meson velocities, these being due to the existence of à 
Brillouin gap in such materials. It might be that in insulators the energy loss 
near zero energy arises from mechanisms other than those considered by FERMI 
and TELLER, e.g. by inelastic collisions with the lattice. The dependence of 
such losses on Z, if any, is altogether unknown. Another theoretical possi- 
bility to explain our observations is that the mesons could be captured near 
zero kinetic energy not by any individual atom, but into a molecular orbital 
around a whole grouping of atoms. From such an orbital they would then 
transfer to individual atoms on a statistical basis. We hope to return to these 
points in more detail in a future publication. 

Lastly, we wish to discuss the relevance of the « transfer mechanism » 


(1) Stopping powers calculated from tables of proton ranges given in UCRL-2301. 
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suggested by DAY and MoRRISON (7) to explain the anomalously high Auger 
rates observed in the K- and L-series of mesic X-rays from light elements ({). 
One could think that the mesons are initially captured according to the 
« Z-law » and get later transferred to the lighter constituents of the compounds. 
This possibility appears, however, rather improbable as the heavier consti- 
tuents, both theoretically and experimentally, exhibit a much smaller fractional 
« éxternal » Auger rate than the light ones. Thus the Day-Morrison mechanism, 
important anyhow only for the lightest elements, would operate in the reverse 
direction. 


RIASSUNTO (*) 


Fermi e TELLER hanno previsto, nel loro lavoro classico sul processo di rallenta- 
mento dei mesoni negativi, le relative probabilità della cattura di quest’ ultimi da parte 
dei differenti atomi di un composto chimico, concludendo che queste probabilità di 
cattura atomica debbano essere proporzionali alle cariche nucleari Z; degli atomi con- 
siderati. Gli esperimenti qui descritti furono eseguiti allo scopo di verificare tali pre- 
visioni: sono basati sul fatto che i muoni negativi raggiungono rapidamente le orbite 4 
attorno ai nuclei degli atomi della materia nella quale sono rallentati e scompaiono 
da queste orbite (per decadimento o per assorbimento) in misura caratteristica per i 
nuclei in questione. Pertanto, la distribuzione nel tempo degli elettroni di decadimento 
emergenti da un composto chimico nel quale si arrestano dei muoni presenta in gene- 
rale tante componenti esponenziali quanti sono i suoi costituenti; le loro intercette al 
tempo zero (ingresso dei muoni nella materia) sono direttamente proporzionali à n,p;, 
dove n, è la frazione rappresentata dalla componente di numero atomico Z;, e p; è la 
probabilità di cattura nello strato mesico À dell’atomo dello stesso costituente. Trascu- 
rando la possibilità del passaggio dei muoni tra atomi vicini, le p; equivalgono alle pro- 
babilità di cattura previste da Fermi e Teller. I presenti esperimenti mostrano che, 
almeno negli isolanti (A1,0;, P,0;, SiO,, CCl,, KOH, KHF,, p-diclorobenzene), le pre- 
visioni di Fermi e Teller falliscono. Non sembra intervenire nessuna ponderazione di 7; 
da parte di Z;; in alcuni casi, anzi, l'elemento più leggero sembra preferito. Come pos- 
sibile spiegazione si fa rilevare che la cessione di energia alle vibrazioni del reticolo, 
non considerata da Fermi e Teller, può essere negli isolanti un importante meccanismo 
per la cattura di particelle negative dotate di energia cinetica quasi nulla. 


(*) Traduzione a cura della Redazione. 
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(ricevuto il 17 Ottobre 1957) 


Summary. — The direct exchange-energy Æ of two atoms in the quasi- 
crystalline picture of liquid argon is calculated. The result is 2.,,/E = 1.4. 
It is Shown that the atomic motions in the liquid, enhanced with respect 
to the solid by the mutual correlations, can explain the differences expe- 
rimentally determined in the structure of the two states. 


1. — In order to investigate the physical consistency of the assumption that 
the elementary diffusion process in a non-associated liquid, near the melting 
point, consists in an exchange of position between two molecules, we have 
evaluated the energy associated with this exchange (1). We have attributed 
to the liquid the f.c.c. structure which, in the case of argon that we consider, 
is also the structure of the solid. Obviously, the numerical result obtained de- 
pends essentially on this hypothesis, even if mitigated by the use of a liquid 
region with an extension of only few interatomic distances. 

The idea of attributing to the liquid state a close-packed structure might 
be suggested by the fact that, at temperature not much higher than the 
melting point, X-ray analysis shows a nearest neighbours’ number of about 
10 for every liquid examined, which does not present molecular association, 
independently of the structure of the corresponding solid (?). This is consistent 
with a f.c.c. structure in which large and correlated molecular motions 
decrease the order-degree and the number of nearest neighbours. This 
point will be discussed again in Sect. 4. 


(1) A similar calculation has been made for copper in the solid state, by F. SEITZ 
and H. B. Huntineron: Phys. Rev., 61, 315 (1942). The activation energy found was 
too large with respect to the experimental value. 

(2) N. S. Ginericu: Rev. Mod. Phys., 15, 90 (1943). 
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2.— Fora f.c.c. lattice with the cubic 
unitary cell edge of length a = 5.716 A, 
determined from the density of liquid 
argon (d—1.420 at the vapour pres- 
sure p=0.8 atm. and temperature 
T= 84.6°K), near the melting point 
(Th = 83.8°K), we have calculated the 
difference, H, between the energy of 
the saddle-point configuration for the 
rotation of the atom pair indicated by 
1 and 2 in Fig. 1, and the energy of 
the equilibrium configuration, given by 
the perfect lattice. Among the possible 
rotations of the atom pair 1, 2, we 
have considered the rotation in the 
(0, 0, 1) plane through the cell-center, and about the [0, 0,1] axis which 
bisects perpendicularly the line 1,2, because the corresponding saddle-point 
configuration has the least energy. 


The nearest neighbours’ distance in the f.¢. e. lattice is (4/2/2)a = 4.04 A. 

In the saddle-point configuration, which is reached after a rotation of the pair 
of 90°, we shall let the displacements, in the [1, 1, 0] direction, of atoms 1 and 2 from 
their interstitial positions, situated respectively in the middle point of the edge 7,10 
and in the centre of the cubie cell, be ha(4/ 2/2): here 1 is a parameter to be determined. 
Let similarly w(a/2) be the displacements of atoms 3, 4, 5, 6 from the equilibrium 
positions in the rotation plane and »(a/2) those of atoms 7, 8, 9, 10 in the [0, 0, 1] di- 
rection. 

We consider also the displacements of A,, A,,....4,, and B,, B,, ... By, in the di- 
rections along the lines joining their equilibrium positions with the positions of 3, 4, ... 10 
in the saddle-point configuration. We write the displacements of the atoms 4;, which 
are all equal, as za(4/2/2), and the displacements of the atoms B,, which are only 
approximately equal, as (z+dx)a(4/2/2). da is small with respect to the other para- 
meters and may be expressed through w, » and a. | 

We assume an interaction potential between two argon-atoms at distance r, of 
the Lennard-Jones type: 


; rama) 


with: 


e = 165.3-10-* erg, o = 3.4054, 


determined by A. MicneLs from the second virial coefficient (?). 

We describe by means of (1) the interactions of atoms 1, 2,..., 10 (Fig. 1) with 
the others and in particular with the A,’s and B,’s: the interactions of each 4, and B, 
with the set of its other 11 first neighbours is described instead by means of the averaged 


(3) A. Micners, Hus. WiskER and H. K. WIJKER: Physica, 15, 627 (1949). 
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potential 11@(r) used in the Lennard-Jones and Devonshire cell-theory. If 7 is the: 
displacement of any A; or B, from its equilibrium position, this averaged potential is 


bo 


( 


) 11 G(r) = 11% 


(es eee 
toa nue (( a) ( ii "tf 


with 


In the saddle-point configuration the mutual distance between atoms 1 and 2 is 
r, = (1— 2A)a(4/2/2). The distance between atoms 1 and 3, and those of the same 
type i indicated by the symbol <1 — 3)), are: r, = [(14+w+4)?+4?}}(a/2) and, with 
a similar meaning of the symbol 


ARTT de 24713, o a = 
2 
C3 — A;): ra = [(1—w)?4 ne CT A5): 7, = (24 pato 
2 2 2 
> | a 
{3 — 49: r= (14 4a be | ha og ST Lew) UE > + (x + dr)a — 
x | 2 1/2 
3— D: re=[(1+ 2)?+ (+), | mo =m, ru (a+ dla, 
| 2 


In this configuration, we have also for each B; 


dV dd 
(3) = —1i | 
dr r=Ty dr r=(7+d7) uns 
For the analogous relation available in the case of A;, eq. (3) becomes 
TUA do 2 
(4) | :) ar = —u(F 3 Nada. 
RE dr? Ja Van 2 


From eq. (4) we have 


dsl fes i lois 4{È 
T= 2 (» = u) 7" iy 19" 9 (2 - fl . 2 5 


to the first order in «/2—2z and with a parabolic approximation to ®(r). 


The energy difference Æ between saddle-point and equilibrium con- 
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figuration if one considers only changes in interaction energies between im- 
mediate neighbours is obtained determining the minimal value of 


DI 
(0) E(A, pv, 7) = V(r) + 4V (13) + 4V (15) — 23V ( VE) - 
: E E TE © 
+ 167 (r,) + 2V (rs) + 8V(rs) + 16V(r,) — 427 (a C+ 
> 
alg 12V(r;) + 16V (rs) — 28V CE + 160@(r,;) De 1607 (a + 


ae 1060(r,,) — 1067 (« a). (9) 


The result is 


E = 796-10-* erg dor 4007, 


LE 1019 
P= 0257, 
0:02 


(Of Course the value of Æ might be lowered if one considers fluctuations 
in the density.) 


3. — A comparison of the H-value, calculated in this way with the experi- 
mental H-value determined through the Arrhenius form of the diffusion co- 
efficient 

D = A-exp[— E/kT] 


is not possible because of the scarcity of experimental data (*). 
Although the elementary processes for momentum transmission and for 
molecular diffusion might be completely different, we observe that E. /E— 0.67. 


vis 


with E, = 529-10- erg derived from the experimental values (5) of the vis- 
cosity coefficient in the form 


n= B-exp (#,,/kT] . 


(*) The coefficients of ©(r) in this expression are equal to the number of links 
between the A; and B; and their first neighbours in the equilibrium configuration. In 
this way we consider also the fact that both B, and B, interact with 7 and 8, and 
analogously for the symmetricals. 

(4) In the work of J. CorBerr and J. H. Wane: Journ. Chem. Phys., 25, 422 (1956); 
only one point of the curve D = D(T) is determined. 

(5) N. S. RupeNKo and H. H. ScHuBNIKOV: Phys. Zeit. Sowiet., 6, 470 (1934). 
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From experimental data concerning the liquids, one has the general empirical 
relation HE,,,,/H~3. In the case of Argon EÆ,.,/atom — 1103-10-%erg at 
T = 87.3 °K and if we use the calculated Æ-value, we have the ratio E,,,,/E=1.4. 

We may say that this calculation indicates that the direct-exchange of 
2 atoms might be the basic mechanism for the molecular diffusion in the 
liquid, even though, obviously, it is not possible to exclude a vacancy motion 
or more complex processes involving a number of molecules higher than two. 
Against these more complex processes one can only say that they have a 
higher activation energy. Let us now examine in more detail this statement. 
For the high density of the solid, the greater contribution to the saddle-point 
energy in the exchange of two-molecules is given by the repulsive terms. If 
one takes into account exchange, for example through 4 atoms rings, one 
obtains repulsive terms with much smaller values than in the previous case 
and hence a smaller saddle-point energy (6). On the contrary, owing to the low 
density of the liquid, the larger contribution to the activation energy is not 
given by the few repulsive terms but by the many attractive terms. It follows 
that when one considers exchange through rings with 3 or more atoms, one 
has negligible decrease in the values of each repulsive term, while the increased 
number of attractive terms rises definitively the activation energy. We have 
verified this with a rough evaluation for the case of a three-atom ring (1, 2, 9 
in Fig. 1) wheeling in the (1, 1, 1) plane, with the possibility of a small shift 
from it. Perhaps more definite information about the presence of these mo- 
lecular ring motions will be supplied by experimental techniques such as 
neutron spectrometry. 


4. — It is possible to believe that the essential difference between the liquid 
and solid states consists in the correlation of molecular motions, which is very 
strong in liquids, but negligible in solids, as is confirmed by the good results 
of the Einstein-model. 

The order-degree, described by the radial distribution function g(r), is 
smaller in liquids than in solids. This molecular disorder in liquids might 
be attributed to the larger displacements of the molecules from their equi- 
librium position, caused by the mutual correlation in their motions: the direct 
exchange of the position of two molecules is the most evident example of this 
correlation. We recall that X-rays analysis in argon at 7= 84.4 °K (7), gives 
the g(r) function. According to the upper limit o of the integration interval 
(0, 0), chosen in order to evaluate the area under the first peak (at r= 3.79 A) 
of the 4z7?g(r) function, one finds a first neighbours’ number lying between 
10.2 and 10.9. If we assume o = 4.69 A, the position of the first minimum, 


(6) C. ZENER: Ac. Orys., 8, 346 (1950); (Evaluation for copper in the solid state). 
(?) A. EISENSTEIN and N. S. GINGRICH: Phys. Rev., 62 (1942) 627. 
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any molecule placed at a distance larger than this from another fixed mole- 
cule can no longer be considered as its first neighbour. Now, in the model 
chosen 50 atoms are involved, essentially, for each elementary process. In 
the saddle-point configuration one has (in angstrom): 


IR 
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If one assigns to each atom of this group the first neighbours’ number found 
for this configuration, one has five f.n. for the atoms 1 and 2, nine for the 
atoms 3, 4, ..., 10, eleven for the atoms 11, 12, ..., 18 and twelve for the re- 
mainder. The averaged f.n. number for the 50 
atoms of the process comes out equal to 11. Natur- 
ally, if the model has to be consistent, the elementary 
processes of diffusion must be a prevailing cause of 
the diminution of the first neighbours’ number with 
respect to the solid. On the other hand, such a 
diminution is not sufficient to obtain the exper- 
imental co-ordination number 2, if we wish to ex- 
clude a frequency of these processes which is so 


SE x t î 
o large that there is appreciable interference amongst 
Fig. 2. them. Therefore it is necessary to consider also the 


contribution to the decrease of given by the large 
oscillations, without exchange, about the respective equilibrium positions. 

In order to discuss this point we use the characteristic methods of the 
Lennard-Jones and Devonshire cell-theory (Fig. 2). We consider an atom that 
during an oscillation shifts by r from the center of the spheric cell with a radius 
a(V2/2) (Fig. 2). On the cell surface are uniformly distributed the 12 first 
neighbours. 

If rLa(V2/2) > 0, all the points of a spheric «calotte» with surface 

ES GAIA IL fesa ; 


4/2 ar 


have a distance greater than o from the oscillating atom. The first neighbours’ 
number for the atom so shifted from the sphere center is lowered to 
2(r) = 12(1 — (o/2xa?)). The averaged first neighbours number #(r) for the 
atom in oscillation, placed in the field of force described by eq. (2), results 
r=Q— ay 3/2) is 
zi | exp Bs 12@(r)/kT] dr + |< ) exp De 12D(r JKT] dr 
a(n) = —"* DR ae eS ee 
Jesp [— 122D(r)kT]dr 
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From this result it follows that the first neighbours’ number of the atoms 
11, 12, ..., B; and A; might be lowered by 0.6. In this way the averaged f.n. 
numbers for the 50 atoms of the process results to be 10.6. 

With a more refined cell-model, for example with no hard walls, it is 
certainly possible that this number can be reduced. In this way the number 
of elementary processes per unit volume necessary for obtaining the mean 
experimental number of first neighbours can be sufficiently low to give no 
appreciable interference between these processes. 

In any case, from this estimate, it is clear that in order to explain the 
number of first neighbours in the liquid state, in particular for Argon, it is 
sufficient to consider the effect of the molecular motions. 


The author would like to thank Prof. J. DE BOER for his kind hospitality 
at the Instituut voor theoretisque Physica of Amsterdam and Prof. G. CARERI 
for helpful discussions. 

He also wishes to thank Prof. J. F. Koksma and Prof. VAN WIJNGAARTEN 
of the Mathematisch Centrum of Amsterdam, who have carried out the mi- 
nimum value calculation for the function H(A, 4, r, 7), contained in this note. 


RIASSUNTO 


Si è calcolata l'energia £ per lo scambio diretto di due atomi, nel modello quasi- 
cristallino dell’argon liquido. Il risultato è E,,,,/£ ~ 1.4. Si mostra come i moti degli 
atomi nel liquido, esaltati rispetto al solido dalle correlazioni mutue, possono spiegare 
le differenze determinate sperimentalmente nella struttura dei due stati. 
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Summary. — Water loaded emulsions have been used to study cut-off 
energy, spectrum and flux of primary particles at geomagnetic latitude 
and longitude 47 °S, 135 °W. Results agree with Stormer cut-off energy 
for the latitude. Most interactions in the water loaded emulsions were 
shown to occur in light nuclei. 


1. - Introduction. 


A number of balloon flights (1) were made during 1954-1956 in Victoria, 
Australia, at approximate geomagnetic latitude 47°S, dip 68°, and longi- 
tude 135° W (geographic longitude 145° E). A study of the fragmentation 
events of the heavy primaries recorded during these flights is of considerable 
interest in view of the recent work of WADDINGTON (?) and SIMPSON et al. (3) 
on the effective geomagnetic co-ordinates for cosmic rays. 600 um Ilford G5 
glass-backed emulsions were used and a large proportion of these were swollen 
to approximately three times the normal thickness (k —3) with water and 
exposed while wet. Details of the method will be described elsewhere but the 
main advantages for this particular investigation are that the paths of part- 
icles passing through the water-loaded emulsions are three times as long as 
in a single normal emulsion and that the chance of a collision with light nuclei 
including hydrogen is considerably increased. Expansions were limited to 
about three times since with greater expansions minimum ionization tracks 
could not be efficiently observed. 

Sixty collisions of heavy primaries with nuclei in the emulsions have been 
analysed, 55 in the swollen emulsions, 4 in dry emulsions and 1 in the cello- 


(1) J. E. LABY, Y. K. Lim and V. D. HoPPER: Nuovo Cimento, 5, 249 (1957). 

(2) C. J. WADDINGTON: Nuovo Cimento, 3, 930 (1956). 

(8) J. A. Simpson, K. B. FENTON, J. Katzman and D. C. Rose: Phys. Rev., 102, 
1648 (1956). 
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phane packing between two emulsions. The charge and energy of à primary 
particle were inferred from its d-ray density combined with the grain count, 
angular distributions and nature of the disintegration fragments. In some cases 
scattering measurements could also be reliably made. From the zenith angle 
of the primary particle a correction was made to deduce the energy at the . 
top of the atmosphere. Seven heavy primary particles were arrested in these 
emulsions and on the assumption that these particles had not changed charge 
an estimate was also made of their energy at the top of the atmosphere. From 
these results it is inferred that the cut-off energy at this latitude and longitude 
is 0.6 to 0.8 GeV/nucleon. 


2. - Energy determination. 


The interaction of an energetic nucleus with a hydrogen nucleus in the 
emulsion is the same as that due to an energetic proton bombarding the heavy 
nucleus at rest, but viewed from a different frame of reference. There are 
important advantages from a study of the former interaction, namely, the 
products of the disintegration can be more readily identified and the energy 
of the primary particle can be estimated from the angular spread of the frag- 
ments. If the energy of the primary particle is above 1 GeV/nucleon there is 
a possibility of mesons being created but these have a wider angular distri- 
bution than the fragments of the primary and can usually be separated from 
them. 

In the case of a collision between a heavy primary particle with an atom 
other than hydrogen in the emulsion, the study is more complicated in that 
fragments are also emitted in the disintegration of the target nucleus; further- 
more, the mesons created may be confused with the knock-on protons from 
the target nucleus. It is however in general possible in these events to sepa- 
rate most of the fragments of the incident nucleus from the mesons and the 
fast protons from the target nucleus. 

The charged fragments consist mainly of protons and «-particles and one 
common method of estimating the momentum p per nucleon of the primary 
particle is to use the formulae (*”), 


(1) (x-particle) 0==10:050)/p; 
(2) (proton) fiala oe 
where 6,,., the root mean square angle, is measured in radians and pe in 


(4) M. F. KAPLON, B. Peters, G. T. REYNoOLDS and D. Ritson: Phys. Rev., 85, 


F. HANNI: Helv. Phys. Acta, 28, 345 (1955). 


) 
1 
(5) K. GorTSsTEIN: Phil. Mag., 45, 347 (1954). 
) 
) Y. EISENBERG: Phys. Rev., 96, 1378 (1954). 
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GeV/nucleon. The formulae are based on the assumption that the «-particles 


and protons observed as fragments of the primary particle are evaporation. 


particles emitted isotropically in the centre of mass system of the primary 
particle with mean kinetic energies of 10 and 12 MeV respectively. 

It became evident during the analysis that where both formulae are ap- 
plicable they do not give identical results. This might be expected as a fraction 
of the protons are emitted in high energy disintegrations with much greater 
energy than can be attributed to a pure evaporation process. An analysis 
of high energy proton and neutron induced stars in normal emulsions by 


PARSONS et al. (8) shows that one-third of the protons with above minimum: 


ionization have energies above 50 MeV and cannot be attributed to a pure 
evaporation process. The energy distribution of protons emitted in the dis- 


integration of a light nucleus with production of mesons in normal emulsions 


obtained by these workers, after correcting for the presence of x-particles and 
heavier nuclei using fragmentation probabilities derived later, is given in 
Table I. 


TABLE I. — Energy distribution of protons from high energy stars induced by protons 
and neutrons. 
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It is possible to estimate from this distribution the angular distribution 
of protons emitted in the fragmentation of a light nucleus bombarding a sta- 
tionary hydrogen nucleus on the as- 


Abb sumption of isotropic emission in the’ 

centre of mass system of the primary 
12F particle. The result for a primary 
i energy of 1 GeV/nucleon is given in 


Fig. 1. It can be seen that a peak 
occurs at 2° to 6°. The root mean 
Square angle occurs at 8.9° which is 
2.1 times the value expected from 
Eq. (2). Agreement with Eq. (2) may 
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Fig. 1. — Angular distribution of protons 
De from fragmentation of a light nucleus with 
0 5 50 =“ energy 1 GeV incident on hydrogen. 


(5) R. W. Parsons, F. A. BRISBOUT and V. D. HoPPER: Phys. Rev., 95, 193 (1954). 
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TABLE II. — Summary of events recorded on flights. 
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TABLE II (Continued). 
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be obtained by neglecting the outer third of the protons in computing 6 


rms ? 


as these are probably due to the non-evaporation protons. For collisions of 
primary particles with atoms heavier than hydrogen there is the possibility of 
still greater deviation from the formulae due to additional collisions. A similar 
method was applied to these events and the angle obtained by this method 


was combined with the 0 


primary particle. 


rms 


for «-particles to obtain the energy of the 
The justification of this method is commented on later. 
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The grain densities and in some cases relative scattering of the secondary 
particles were taken into account in the energy estimation when the primary 
energy is low. 

To correct for the ionization loss in the matter traversed before collision, 
use has been made of the stopping power given by BRADNER and BISHOP (?) 
for glass, the extrapolation of the stopping powers given by WILKINS (1°) for 
dry emulsion and water and of the energy-range relation given by Rossi (1) 
for air. The stopping powers for glass, normal emulsion and 3 x expanded 
emulsion, defined as ratio of range increment in g/cm? relative to air are res- 
pectively 0.96, 0.75 and 0.87. The corrections have been made individually 
for each primary particle and are listed in Table II. 


3. - Charge determination. 


It will be shown that in the loaded emulsions used, most of the interactions 
took place with the light nuclei of the target. This has the effect of limiting 
the target prongs to a small number and so reducing the possibility of target 
prongs being confused with the fragments of the primary particle. Thus it 
was possible to identify the primary nucleus from its fragments with reason- 
able certainty particularly in the several events where collisions occurred with 
hydrogen. These events served as the basis of ¢-ray calibration from which 
other particles could be identified. 

Results for the 60 primary particles undergoing fragmentation and 7 pri- 
mary particles that ended in the emulsions are shown in Table II. The stack 
flown on 26 March 1956 was area-scanned and every star in 116 cm? of 2.6 X 
loaded emulsion and 39 em? of normal emulsion was examined carefully to 
determine whether it had been produced by a particle of charge greater than, 
or equal to 2. The remaining stacks were scanned for heavy primary dis- 
integrations, mainly by the method of line scanning. 


4. - Distribution of target prong number. 


The distribution of target prong number N, for events in the loaded emul- 
sions (mean k= 2.6) are compared in Fig. 2 with the distribution obtained 
for normal emulsions. While the loaded emulsion data give a much smaller 
percentage of events with large prong numbers, the distributions at lower 


(9) H. BRADNER and A. S. BisHop: Phys. Rev., 76, 587 (1949). 
(109) J. J. WILKINS: A.E.R.E. Report G/R 664 (1951). 
(1) B. Rossi: High Energy Particles (New York, 1952), p. 40. 
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prong numbers are strikingly simula: 
This would be expected if in the normal 
emulsions the events in heavy target 
nuclei tended to occur with large prong 
numbers and in the loaded emulsions 
most of the events took place in light 
target nuclei. In particular the sudden 
decrease of frequency at N,=7 for 
both the normal and loaded emulsion 
distributions suggests that few light 
target nucleus events occur with 
N, > 7. The distributions are also com- 
pared in Table III. 

If it is assumed that events with 
N, >6 arise entirely from interactions 
in the heavy nuclei of the emulsion, 
one would expect the number of such 
events in a given area to remain unaf- 
fected by loading. The numbers in the 
table are therefore normalized so that 
the number of these events is equal 


10 7 Fig. 2. — Distribution of target number: 
(a) normal emulsion, Noon and KAPLON (!?); 
5 (b) normal emulsion, FOWLER et al. (13); 
dmn n Vp (c) loaded emulsion, present work. 
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TABLE III. — Distribution of target prong number. 


Observed frequency Normalized 
Geomagnetic | Primary of stars frequeney of stars 
Emulsion latitude charge with N, with N, 
20 3 er É 

MO eo GE Gu RO AE 0 = EC 

Normal (12) | 41° and 55° N > 6 37 117 72 | | ; 
Normal (18) | 46° N = 3 61 115 MIEr42 4] di | oO og 
Loaded 47° $ 23 19 44 11 173 | 400 100 


(2) J. H. Noon and M. F. Kapton: Phys. Rev., 97, 769 (1955). 
(18) P. H. FowLER, R. R. HirLieRr and C. J. WADDINGTON: Phil. Mag., 2, 293 


(1957). 
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to 100 in each case. These results indicate that about 50% of the events in 
these loaded emulsions are due to the oxygen and hydrogen nuclei introduced 
by loading. It is of interest to note that for the flight on 26 March 1956 the 
loaded emulsion area of 116 cm? gave 24 events whereas one third of this area 
also carefully area scanned in the normal emulsion gave only two events. If 
one assumes that the gelatine and water nuclei have the same cross-section 
for the fragmentation of heavy primaries, the data of Table III combined 
with emulsion compositions show that of the events observed in the loaded 
emulsions 75% occurred in light nuclei; the corresponding value for the normal 
emulsion being 45%. The percentages calculated from the emulsion compo- 
sitions are 70% and 40% for the loaded and normal emulsions using the cross- 
section. formula given by BRADT and PETERS (14): 


(3) o=nm(R;+R,—2AR}, 


R; and R, being the nuclear radii of the incident and the target nuclei given 
by 7,A*, ro = 1.45-10-1 em, AR = 0.85-10-1* cm. The same process may be 
applied to events with N,=0 and 1. It is found that 33% of the events in 
the loaded and 16% in the normal emulsions are due to interactions with 
hydrogen nuclei as target whereas the theoretical formula of Bradt and Peters 
gives 25% and 14% respectively. In the loaded emulsion 90% and in the 
normal emulsion 75% of the events with N,<6 are interactions with light 
nuclei. The hydrogen events account for about 85% of the normal emulsion 
events and 95% of the loaded emulsion events with N,= 0 and 1. Thus it is 
reasonable to attribute the peaks at N,= 0 and 1 and N,% 5 in the distri- 
butions to interactions in hydrogen and C, N, O nucleus targets respectively. 


5. — Fragmentation probabilities. 


Since 90% of the fragmentation events with N, < 6 observed in the loaded 
emulsions took place in light nuclei, the fragmentation probabilities derived 
from these must be characteristics of the collision of a heavy primary particle 
| with a light nucleus of the emulsion or water. Detailed knowledge of frag- 
mentation probabilities with light nuclei is required to estimate the nature of 
the primary flux from measurements of flux at balloon altitudes. The frag- 
mentation probabilities P,,, i.e. the number of J-type nuclei emitted in the 
fragmentation of an J-type primary particle, are given in Table IV. The 
nomenclature of BRADT and PETERS (14) has been followed: H-nuclei, 10 <<; 
M-nuclei 6<2<9; L-nuclei 3<2< 5. 


(14) H. L. BRADT and B. PETERS: Phys. Rev., 77, 54 (1950). 
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TABLE IV. — Fragmentation probabilities P,, from events with N,<6 in loaded emulsions. 


Gev/nucleon 


Number —%» 


0 5 1.0 5 2.0 
GeV /nucieon 


Gev/nu cleon 
A 
1.5 


2.0 


Fig. 3. — (a) Differential energy spect- 
rum of primary particles which pro- 
duced fragmentation; (b) lower end of 
spectrum of Fig. 3a; (c) energy spect- 
rum at top of atmosphere of particles 
which ended in the emulsion. 


J 

ee? H M L x p 
H 0.19 + 0.11 | 0.25+0.13 | 0.31 + 0.14 | 2.25 + 0.38 | 4.44 + 0.53 
M — 0.11 + 0.07 | 0.47 + 0.16 | 0.79 + 0.20 | 2.90 + 0.39 
L — — 0.40 + 0.20 | 1.00 + 0.32 | 0.90 + 0.30 


The probability for L-nucleus emis- 
sion in the fragmentation of H- and 
M-nuclei P;,yx is 0.40 + 0.11 and 
the «/p ratio is 0.45. 

The values may be compared with 
those obtained by NooN and KAPLON (1?) 
and by FOWLER et al. (1) for events. 
in normal emulsions with N, <7. The 
probabilities of H- and M-nucleus emis- 
sion are not significantly different from 
those listed by these workers, whereas 
for the L-nucleus and «-particle emis- 
sion our results agree with those of 
Noon and KAPLON and correspond to. 
a higher probability for L-nucleus emis- 
sion and a lower probability for «-part- 
icle emission than the Bristol results. 


6. — Energy spectrum and geomagnetic: 
cut-off. 


The differential energy spectrum, 
where for each particle the ionization 
loss has been corrected to correspond 
to energy at the top of the atmosphere, 
is Shown in Fig. 3a. The lower end of the 
spectrum is given in 0.1 GeV/nucleon 
intervals in Fig. 3b, the number inside 
each square giving the charge of the 
particle. Seven tracks of energetic heavy 
nuclei have been observed to end 
in the emulsions. Their charges were 
determined from the thin-down lengths 
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and energies estimated by the range-energy relation. The energy distribution 
after correction for ionization loss in the overlying air on the assumption that 
they occurred at the ceiling height and suffered 
no loss of charge before reaching the detector is 


given in Fig. 3c. al 

The plates were attached to hang vertically 12 
and zenith angles were measured on the as- i 
sumption that they did not swing during the gb 
level part of the flight. Visual observations sup- 8 
port this contention. In addition, the measured ST 
zenith angular distribution of the primary part- 3.6) 
icles which produce fragmentation in loaded dl 
emulsions is shown in Fig. 4, and it is seen 3 
that very few particles had calculated zenith 2 [ 


—L LP _L__L__ L 


angles greater than 60°. The smooth curve is Cea ani 
the calculated angular distribution of the in- SE 

tegrated track length based on the assumpt- Fig. 4. — Zenith angular distrib- 
ion of isotropy at the top of the atmosphere  Ution of the primary particles 
and normalized to the total number of inter- E EL eo ia 

3 E Thesmooth curve corresponds to 
actions. This result supports our argument the calculated angular distrib- 
that the path of the primary particles may ution, assuming isotropy at the 
be assumed straight in passing through the top of the atmosphere. 
bulk of the atmosphere above the plates. 

The spread in the lower end of the energy spectrum shown in Fig. 3 may 
be due to the following: (a) the geomagnetic cut-off energy varies with di- 
rection; (b) the correction for ionization loss is made on the assumption that 
the primary particles reached the plates without loss of charge, but some were 
the secondary particles of some heavier nuclei. The corrections for these would 
be inadequate, since the ionization loss in energy per nucleon depends on 2R, 
the product of charge and range. That this is an important factor giving rise 
to the observation of particles with energies apparently lower than the cut-off 
may be seen from Fig. 3b, where most of the lowest energy particles appear 
to have low charges. On the assumption that most of these were products 
of the fragmentation of the H-nuclei incident with near cut-off energies, this 
effect would account for about 5 of the low energy particles. 

Thus from our results we conclude that the geomagnetic cut-off at geomag- 
netic latitude 47° S and geographic longitude 145° E occurs at 0.6 —0.8 GeV /nu- 
cleon. This cut-off is in agreement with the energy distribution of heavy nuclei 
arrested in the stack shown in Fig. 3c. 

The smooth curve in Fig. 3a is given by 


dN/dE ~ (E + 0.94)? ; yw 2A 
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the area defined by the curve and H= 0.7 GeV/nucleon having been normalized 
to equal the total number of particles. The exponent is in agreement with 
the values given by other experimenters (7). 

In view of the recent results of FOWLER et al. (13) which suggest that the 
mean kinetic energy of evaporation «-particles in the C-system of the incident 
nucleus increases with the charge of the incident particle, we have divided 
the events studied into two groups according to the primary charge 3 < 2 < 7 
and 2 >8 and considered their distributions in energy after ionization cor- 
rection. The two distributions were not significantly different and gave the 
same cut-off energy. We have also considered separately events with WN, = 0 
and 1 which are mainly due to collisions with hydrogen nuclei and found this 
gave the same cut-off energy within the statistical error and this justifies the 
extension of the method of estimating energies to events N, > 1. 


7. — Flux of primary particles with 2 > 6. 


It is possible to make a rough estimation of the flux of primary particles 
with 2 > 6 above Victoria from the number of fragmentation events produced 
by these particles using the interaction cross-sections calculated with the for- 
mula of Bradt et Peters (Eq. (3)). The calculated mean free paths of the 
M- and H-nuclei (3— 7 and 3— 15) are respectively 26.5 and 18.0 g em? in 
air (12) and 33.2 and 22.0gem-= in the loaded emulsions exposed on 
26 March 1956. 12 M-nucleus and 6 H-nucleus events were observed in an 
area of 97 cm? of emulsions of this stack swollen to 2.62 times the original 
thickness. With the fragmentation probabilities tabulated in Table III, which 
serve aS an approximation to the case of interaction with air, the fluxes of 
primary M- and H-nuclei at the top of the atmosphere are evaluated as 
(5.7+1.6) and (3.7+1.5) particles m°? sr-! 8-1 respectively and these would be 
sufficient to create all seven L-nuclei events observed on the plates. The 
number of events occurring during ascent and descent has been estimated to 
be only 2% for these particles, which justifies the assumption made in the 
ionization loss correction that the events took place at the ceiling height. 


8. — Discussion. 


In a recent experiment of WADDINGTON (?) in which two stacks of plates 
exposed at the same geomagnetic latitude 55° N, one over England and the 
other over North America were studied, it was shown that the flux of primary 
a-particles over England was significantly lower than the flux over America. 
The cut-off energy over America from scattering measurements was estimated 
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as less than 0.15 GeV per nucleon and in England the method gave a value 
between 0.6 and 0.7 GeV per nucleon. The difference in flux could be attri- 
buted to the difference in cut-off energy. The cut-off energy calculated from 
the known geomagnetic latitude was 0.3 GeV per nucleon. If a shift in the 
geomagnetic latitude was entirely responsible for the change in the cut-off 
energy, it was pointed out that the observations would be consistent with a 
lowering of 5° and an increase of 3° or more in the accepted geomagnetic lati- 
tude of Europe and North America respectively. Later measurements (FOWLER 
and WADDINGTON (!)) carried out in North Italy gave a cut-off energy of 
1.6 GeV per nucleon whereas from the known geomagnetic latitude (46° N) 
the expected value would be 1.0 GeV/nucleon. 

SIMPSON et al. (*) studied the position of the effective geomagnetic equator 
as defined by sea level cosmic ray intensity minima for different longitudes. 
These workers concluded that the effective geomagnetic equator did not co- 
incide with the geomagnetic equator obtained by surface magnetic measure- 
ments. They suggested that a westward shift of the inclined magnetic dipole 
of the earth, whether centred or off-centred by about 40 — 45° without requiring 
an appreciable change in the angle of incidence would fit the sea level measu- 
rements made by them and other experimenters. This suggestion however did 
not fully explain the high altitude «-particle cut-off energy and flux measu- 
rements. 

The cut-off energy at a given location depends strongly on the angle of 
incidence so that a sharp cut-off value is not obtained. It will be assumed 
that the measured cut-off is given approximately by the cut-off for vertical 
incidence. The Stoermer cut-off for vertical incidence of particles with 2 > 2 
at 47°S is 0.9 GeV/nucleon. On the assumption that the centred dipole field 
of the earth is displaced 45° westward, the vertical cut-off at our location 
would be 1.8 GeV/nucleon. The extrapolation is more complicated if the dipole 
is displaced from the geographic centre. The effect of dipole eccentricity has 
been evaluated by Jory (17). At the geomagnetic longitude of Victoria 
(135° W) the change of vertical rigidity due to this effect is about — 14%. 
Thus the cut-off energy for Victoria based on the model of SIMPSON et al. and 
taking into account the dipole eccentricity is 1.5 GeV/nucleon. Alternatively, 
using the geomagnetic co-ordinates obtained by the surface surveys and 
allowing for dipole eccentricity a cut-off of 0.7 GeV/nucleon is expected. Our 
measurements thus give a cut-off value ranging from 0.6 to 0.8 GeV/nucleon, 


(5) R. E. DANIELSON, P. 8. Friern, J. E. NAUGLE and E. P. Ney: Phys. Rev., 
103, 1075 (1956). 

(16) P. H. FowLER and C. J. WADDINGTON: Phil. Mag., 1, 637 (1956). 

(MR JOR Phys. Rev. 1021671956) 
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which is in agreement with values derived from the accepted geomagnetic 
co-ordinates but disagrees significantly with the Simpson et al. model. 

No significant difference in cut-off value was detected over flights separated 
by several months, so it is possible that the best way of studying the effective 
magnetic field of the earth would be to measure cut-off energies at widely 
carefully selected points. To avoid any doubt of absolute energy determination 
it would be of value to determine the lines of constant cut-off energy and to 
compare these with the geomagnetic latitudes. The use of liquid loaded emul- 
sions has considerable advantage over normal emulsions for such a study as 
they give a much greater density of fragmentation events for determining 
cut-off energies. 
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RIASSUNTO (*) 


Per studiare il cut-off, lo spettro energetico e il flusso delle particelle primarie alla 
latitudine e longitudine geomagnetica 47° S e 135° W si sono usate emulsioni imbibite 
d’acqua. I risultati si accordano con l’energia di cut-off di Stôrmer per la latitudine. 
Si dimostra che la maggior parte delle interazioni nelle emulsioni imbibite d’acqua 
avvengono nei nuclei leggeri. 


(*) Traduzione a cura della Redazione. 
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K Meson Exposures in Water-Soaked Emulsions. 


G. AscoLi, R. D. Hm and T. S. Yoon 
Physics Department, University of Illinois (*) - Urbana, IU. 


(ricevuto il 12 Novembre 1957) 


Summary. — A technique is described whereby pellicle strips of G5 emul- 
sion are soaked in water and exposed as a stack of pellicles in intimate 
contact with one another. This technique allowed the normal process 
of tracing-through to be carried out and complete analysis of events 
was observed to be possible. The performance of the soaked emulsions, 
when they were exposed to a K -meson beam, is described and compared 
with that of unsoaked G5 emulsions exposed in the same beam. 


In the following note we have been concerned with the problem of how 
to increase the frequency of K -free proton interactions in nuclear emulsions. 
Many experimenters (1%) have already used soaked or impregnated emulsions 
for this purpose, but as far as we are aware these experiments have been 
performed with single plates or pellicles. In order to retain as complete 
information as possible on the K -free proton interactions we set out to 
employ a stack of emulsion strips. 

Two stacks of 30 pellicles of two different thicknesses, one of 600 um and 
another of 1000 um G5 emulsions, were prepared, exposed and processed. 
The dry pellicles, immediately after unpacking were weighed and measured 
for size. The pellicles were then placed in distilled water at 4°C. The water 
was slowly circulated over the emulsions and nitrogen gas was bubbled through 
the water at the inlet. It is not known whether this was essential but it was 


(*) Assisted by the joint program of the U.S.A.E.C. and the U.S.O.N.R. 

(1) G. GOLDHABER: Phys. Rev., 74, 1725 (1948); S. GoLDHABER, J. C. BIELK, E. M. 
FRANKL and G. GOLDHABER: A.H.C. Document N.Y.0.-6138, 1953 (unpublished) . 

(2) DULKOVA, ROMANOVA, SOKOLOVA, SUKHOV, TOLSTOV and SHAFRONOVA: Proc. 
Acad. Sci. USSR, 107, 43 (1956). 

(3) C. J. Barry: Nuclear Inst., 1, 138 (1957). 
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felt that the dissolved oxygen in the water should be expelled, especially as 
it might lead to fading of the image between the time of exposing and pro- 
cessing. Care was taken to turn the pellicles over several times especially at 
the start of the soaking. This tended to prevent sticking of the pellicles to 
the paper on which they rested and also might have allowed more uniform 
swelling of the emulsion. i i 

After a soaking period of 12 hours, the pellicles were loaded into a lucite 
box which was specifically designed so as to contain the pellicles without 
allowing space either around the edges or above the pellicle stack. As each 
pellicle was transferred to the box, its area was measured and the average 
thickness of a pellicle was obtained by measuring the total thickness of the 
stack loaded into the box. It may be pointed out that the pellicles are readily 
handled as long as their temperature is kept low. Following soaking, the sur- 
faces of the pellicles become very smooth and slippery and there is no ten- 
dency for the individual pellicles to stick to one another. In fact, a warning 
should be given that unless the pellicles are rather positively compressed within 
the box, they may slide over one another during transport and even between 
the time of exposing and processing. When all the pellicles had been trans- 
ferred to the box and all the interstices had been filled with water, the box 
was tightly sealed with a rubber top plate and a lid of lucite which was screwed 
down. 

During transport and exposure, the box of pellicle was maintained at a 
temperature of approximately 5 °C. During exposure the temperature was 
controlled by two copper plates which were clamped to the box and through 
which was circulated water from a reservoir of melting ice. After exposure 
the emulsions were X-rayed with a grid of lines spaced accurately 0.5 inches 
apart across the surfaces of the emulsion strips. We used X-rays of 150 kV 
filtered by 3 mm of copper and observed a penetration of approximately 
4.5 cm into the wet stack of emulsions. For the X-ray slit system we found 
50 um lead slits of 1 inch thickness followed by 500 um steel slits of 1 inch 
thickness to give satisfactory X-ray lines in the developed emulsion of the order 
100 um width. 

Processing of the emulsions was started as soon as possible after exposure. 
No pre-soaking of the emulsions was, of course, necessary and a normal de- 
veloping procedure was followed, except that necessarily the pellicles were 
unmounted on glass throughout. Although this has the advantage that the 
developing process, because of diffusion through both sides, is speeded up, 
and also any trouble of sticking-down is avoided, it introduces problems of 
uniform drying and of fixing down to glass backing in the later stages of the 
processing. Concerning the finished product of these actual two emulsion stacks 
it must be admitted that trouble was experienced with dichroic fog. However, 
in several previous tests using similar techniques no serious dichroic fog was 
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experienced. A considerable improvement in the transparency of the emulsions. 
was brought about by refixing in ammonium thiosulphate. 

On the average, the shrinkage factor of the dried emulsions measured from 
the wet state in which they were exposed was found to be 5.2 + 0.2. This. 
figure undoubtedly depends to some extent on the manner of drying the pel- 
licles; ours were soaked in increasing concentrations of alcohol and water and 
then finally dried on frames. The dried pellicles were mounted finally on glass. 


TABLE I. — Characteristics of soaked and normal emulsions. 

Average values 1000 um stack 600 um stack 
DhICKNESS AT CN 1.04 mm 0.64 mm 
CATANIA Yi M EN tre 8.5 em x 11.6 cm 8.3 cm x 10.8 cm | 
Densityeedny 0. wet. See er 3.85 g/cm? 3.91 g/cm? 
Lthieknesssswetr > ene 2.16 mm 1.35 mm 
ATEN OCR ey ae TT 10.4 cm x 14.1 em 10.0 em x 13.1 em 
Density wet (5) EE 1.92 g/cm8 1.94 g/cm? 
Volume swelling factor . . . 3.1 3.1 
Thickness swelling factor . . ell 2.1 


(*) The density of the wet emulsion was estimated from volumetric measurements only 
and by assuming that the volume increase of the emulsion due to water uptake was 1 cm* per g 
of water. The validity of this assumption is open to question as it has been shown experi- 
mentally by WALLER (*) and by BATTY (*) thaï the figure of water-uptake varies between 0.84 and 
1.00 cm* per g of water. Our preliminary experiments appeared to show for the batches we 
used that the value was, within 5 %, 1.0 cm? per g. It may also be possible that emulsions in 
pellicle form take up water somewhat differently from those on glass backings, such as have 
been used in other experiments. It is abundantly clear from the information of Table I that 
the swelling of emulsion pellicles is anything but isotropic. 


photographic plates on which had been printed grids of 0.5 mm spacings. The 
X-ray marks on the pellicles were located accurately for each individual pel- 
licle on the sub-stratum of the 0.5 mm grid. The pellicles were cemented into 
position on the glass by a thermo-setting plastic used in a manner similar to 
that already described by YAGODA (4). It has been found that tracing through 
of the tracks from one emulsion to the next is straightforward; rarely are 
the beam tracks displaced from one emulsion grid to the next by more than 
a field under X22 power. 

Data on the uptake of water and swelling of the emulsions are given in 
Table I, and approximate range-energy values, calculated from values for 


(4) H. YAGODA: Rev. Sci. Instr., 26, 263 (1954). 


Usi 


568 G. ASCOLI, R. D. HILL and T. S. YOON 
dry G5 emulsions, are given for protons in our wet emulsions in Table II. 
The minimum grain density observed in the soaked emulsions was (16 +2) grains 
per 100 um, as compared with a normal grain density in dry G5 emulsion of 
the order of 30 grains per 100 um. In some earlier tests on different batches 
of G5 emulsions we observed minimum grain densities in wet emulsions as 
high as (20-22) grains per 100 um. The calculated scattering constant for 
100 um cells in the wet emulsion was approximately (read 14°-MeV/c), and this 
value has been found to agree with the observed multiple scattering of known 
particles of known momenta. 

Some information on the frequencies of events produced by K7 mesons 
in wet emulsions is given in Table III. Also included in the table, for the 
purpose of comparing wet with dry emulsion performances, are some data 


TABLE II. — Comparison of proton ranges in wet and dry G5 emulsion. 
E (MeV) 1.0 5.0 10 20 100 
Range, dry G5 10.5 um 170 um 550 um 1.9 mm 31 mm 
Range, wet G5 (*) 20 um 255 um 850 um 2.85 mm 50 mm 


(*) These values are computed for the compositions of wet emulsions described by the 
characteristics of Table I and in accordance with the assumption discussed in the footnote of 
Table I. They have been checked satisfactorily against the u-meson range from observed 
m-u. decays and also against the observed range of K-mesons of incident momentum of 300 MeV/c. 


TaBLe III. — Comparison of events produced in soaked and normal G5 emulsions by 
K -mesons of initial momentum approximately 300 MeV/c. 
Stars Zero-prong Stars « Stoppings » Elastic 
Type of event at rest stars in flight in flight hydrogen 
() at rest (*) interactions 
Water-soaked G5 500 143 128 32 5 
(62.3%) (17.8%) (15.9%) (4%) (0.6%) 
Normal G5 2557 516 330 69 14 
(74 %) (14.7%) ( 9.5%) (2%) (0.4%) 
-Ratio 0.8 1.2 197 2.0 1.5 
(*) Stars at rest and stars in flight have one or more visible heavy prongs. Stars in flight 
also include decays in flight. 
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by the Bristol workers (*) who used an identical K~-meson exposure with 
normal G5 emulsions to the Barkas beam at the Bevatron, University of 
California. The percentages in parentheses are those of the total numbers 
of K -particle events observed. It is seen from the table that the ratios for 
those types of event which might be expected to show interactions between 
K-particles and hydrogen nuclei are consistently higher for wet emulsions than 
for dry emulsions. This increase is approximately what might have been 
expected from our emulsions which have 5.55-1022 atoms of hydrogen per cm? 
of emulsion as compared with 3.2-10?? atoms per em? of normal G5 emulsion; 
a ratio of 1.74. This figure should also be multiplied by the ratio of the average 
path length of K~-mesons studied in our emulsions to that studied in the dry 
G5 emulsions. The range of a 300 MeV/c K-meson in normal G5 emulsion 
is 3.8 em. The average path length of the K-mesons measured in our emulsion 
was only 3.92 em. Although the range of a 300 MeV/c K-meson in the wet emul- 
sions is 6.0 em, we had a significant loss of range due, firstly, to penetration 
of the K-particles through the wall of the box, secondly, to the spoilage of the 
edges of the emulsions where they were clamped in the drying frames, and thirdly, 
due to the fact that the scanning for the particles was commenced in from 
the edge by about 0.5 cm. An estimated overall figure of improvement of 1.8 
in favor of the wet emulsions over the dry emulsions in yielding K-free proton 
interactions is therefore reasonably consistent with the observations. 

In conclusion it seems clear to us that on account of a number of adverse 
factors, such as the large shrinkage correction required, the presence of larger 
distortions to the tracks, the less well-known range-energy values involved, and 
the lower grain densities of tracks exhibited, the use of soaked emulsions 
for increasing the number of K -hydrogen interactions does not justify the extra 
labor involved in the exposure and processing of the emulsions. It might be 
possible, however, that for special applications such as loading with deuterium 
or other elements the soaking techniques described here are advantageous. 


(*) We are indebted to Dr. D. J. PRowse and other workers at the University of 
Bristol for letting us have these figures from the K -stack collaboration exposure 
which was reported on at the Conference in Venice and Padua, September 1957. 


RTEMSISENARON() 


Si descrive una tecnica per la quale strisce di pellicole di emulsione G5 si fanno 
imbevere d’acqua e poi si espongono come un pacco di pellicole in intimo contatto 
reciproco. Tale tecnica consente di seguire le tracce e di fare l’analisi completa degli 
eventi osservati. Si descrive il comportamento delle emulsioni imbevute d’acqua e lo 
si confronta con quello di emulsioni G5 non bagnate esposte allo stesso fascio. 


(*) Traduzione a cura della Redazione. 
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The Electron-Electron Scattering 
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Recently several experiments have been carried out in order to verify the sug- 
gestion of LEE and YANG, namely the non-conservation of parity in weak inter- 
actions. The purpose of these experiments was to prove the asymmetry taking 
place as a consequence of the parity non-conservation as well as the polarization 
of particles produced in the interaction. From the available experimental data 
no unambiguous conclusions can be drown on the nature of weak interactions, 
mainly of 8-decay. Thus it is most desirable to carry out further experiments on 
several B-active nuclei. Recently a new experiment has been suggested to inves- 
tigate the polarization of electrons released in f-decay of oriented nuclei. Namely 
the emerging electrons have to be scattered on polarized electrons (scattering by 
magnetized iron) and the scattering cross-section have to be measured. From the 
difference of this result and the result obtained from Moller-scattering the polariz- 
ation of electrons emerging from B-decay can be determined (*) (+). 

The idea of this experiment has suggested to investigate theoretically how the 
cross-section of electron-electron scattering will be modified if the electrons are 
longitudinally polarized. The calculations have been carried to the first non- 
vanishing approximation of the S-matrix. The graphs of the problem in the used 
approximation are 


(*) This suggestion occured simultaneously but independently to H. FRAUENFELDER (University 
of Illinois) and to I. BERKES (Central Research Institute of Physics, Budapest). 

(+) After the completion of the present calculation we were informed about the experiment 
performed by H. FRAUENFELDER et al. [Phys. Rev., 107, 643 (1957)]. In this paper the authors 
mentioned, that this calculations were carried out by Dr. BINCER. 
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where p, q are the four-momenta of incident electrons, and p’, q’ the four-momenta 
of the emerging electrons. 


The transition matrix element of the reduced S,-matrix is (1) 


(1) (p', g'|S21p, 9) = Mô(p'+ q'—p— 9), 


(2) M = oie Fr) pp} DA vo) u(p')y,v(q)v(q')y,u(p) 
c (Pip) (gp)? 


The transition probability per unit time is 


(3) w=— > |(p,g'|S81p,9)l2. 


fh final states 
By summing over the spin and momenta of electrons in final states, we find 


iet 
(4) w= aap [raw mc?) 0(q'2 + m2c?) « 
; {(p' — p)-* Uy, (para + ime)yyuoy,(qvo + ime)y,0 + 
+ (a — p)-* Uy, (qua + ime)y,uoy,poro + ime)y,v — 
— (p' — p)*(q' — p)-? Uy, (py, + ime)y,00y,dqYe + ime)yyzu — 


— (p’ — p)-2(q' — p)?uy, (qava + ime)y POY, (Doro + ime)y,u} . 


It is assumed that the electrons in the initial state are polarized longitudinally, 
i.e. the spins are directed in the direction of the motion or in the opposite di- 
rection. Let us first consider the case when spin and momentum are parallel. The 
operator of the spin-projection (1+0,)/2 be introduced as follows: 


1 + op 1 + op 


u(t) = 1); 


2 2 


(5) 


Here u(t) is the spinor belonging to the spin state parallel to the momentum and 
u(}) belonging to the antiparallel state, o, is the scalar product of the spin-vector © 
and the vector-unit p/|p|. By using the spin-projection operator (1+ o,)/2 and 
the energy-projection operator 


(6) Ae ama à 


(1) J. M. Javon and F. RoHRLICH. The Theory of photons and electrons, (Cambridge Mass., 
1955), p. 253-257. 
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where O is an operator composed from y, matrices. Here the relations 


RE 1 
(8) (2) - Dadi A? =A, and PE, |=0 


2 


were used. Using equations (7) and (8) the transition probability (4) can be re- 
written in the form: 


set 1 
aaa ! §(p/2 202) 5 (012 DONC MES. 
Ca saa Le PI Et 


-{(p' — p)-4 Spur ((1 + op)(Pava + ime)y,(paya + ime)y,): 
‘Spur ((1 + 0,)(937g + ime)y,(qove + ime)y,) + 
+ (q' — p)-4 Spur ((1 + op)(Pav_ + ime)y, (Gav, + imc)y,) 
‘Spur ((1 + 0,)(qgvg + ime)y,(Povg + ime)y,,) — 
— 2(p' — q)-*(q' — p)-® Spur ((1 + 0,)(p,7, + ime)y,(pyy_ + ime)y,,: 


“(1 + o,)(dgvg + ime)y, (dore + iME)yy)} - 
From the transition probability we get for the cross-section: 


00 


(10 o . 
! Vit Ve 


v,+ v, is the relative velocity of the particles. In the following this formula will be 
specialized to the center-of-momentum system (CM-system), which is characterized 
by the relations: 


(11) P=-q@ P=-C: |pl=|p'|;: P=%= P= %- 


Owing to the 0-functions the integral in (9) can be carried out easily. Using 


the relation 


(12) Valy + Valy = 25 yy 


the spures occured in (9) can be calculated in the well known manner. Finally we 
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get for the differential cross-section: 


d 2 dd he 
(13) Sh do als 2 [2x Db 2) (x4 nl A 
d2 \dQ/Moner (#2 — 1)? | sin? d 
(in CM-system) 
where 
AME 
Ì dQ/moner  **(@—1î?| sint 6 sin? # | 


is the differential cross-section of Maller-scattering and we used the following ab- 
breviations: 


(15) ied a: È p= vje 
i ~ 1\4rme) ’ 3 Vea gas : 


v is the velocity of electron, 4 the angle between the vectors p, p' in CM-system. 
If the spin of the two electrons is antiparallel to the direction of motion, ope- 
rators (1 —o,)/2 and (1 — 0,)/2 have to be used instead of (1+0,)/2 and (1+0,)/2, 
respectively. The result of the calculation will be the same as in the former case. 
If the spin of one of the electrons is parallel to the direction of motion while 
that of the other is antiparallel, we get 


iw do (do oe [2x2 — 1)(4x2 — 3) 
doe dQ Moller 2? (4? — 1)2 | sin? 9 


a_n 
(x aC 


These two cases can be expressed in one form, namely: 


È do\ (do oc? | (2x? — 1)(4x2 — 3) n all 
oe dQ}, \do re (a — 1)? | sin? Ÿ 4 JE 


The notation (+) indicates the case when both spins are parallel to the momenta 
and notation (—) the case when only one of them is parallel. 

From the experimental point of view the expression (17) has to be also transformed 
into laboratory system. The relations between the corresponding quantities of 
laboratory- and center of momentum-system are as follows (1): 


2A 3) RIT Ve 
2 + (y — 1) sin? 6 


cos Ÿ — 


7 


8 cos 0(y + 1) 


= PS : == 
(18) ea 


Here 0 and d refer to the laboratory-system. We u el y for the abbreviation of 
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the following expression 


1 
(19) pe === 
SAT v2/C2 
v is the velocity of electron in laboratory-system. 


By using the relations (18), for the differential cross-section in laboratory-system, 
we get: 


n do _ (do n 6402 cos 0 | Jey» (y + Dieu 
do \do/moner (—1°[2+(y—1)sin?0]?| 1— a? 4 


(in L-system). 
* K * 


The authors wish to express their gratitude to Dr. G. Marx for suggesting this 
problem. 
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The first one hundred an fifty years; 
A history of John Wiley and Sons, 
Inc. 1807-1957. 1 Volume in-16° 
con 32 figure. J. Wiley & Sons, 
Ine. 1957, pp. xxv+242. 


Questa rassegna che la nota casa edi- 
trice J. Wiley and Sons ha pubblicato 
sulla sua centocinquantenaria attività of- 
fre lo spunto a talune brevi riflessioni su 
la storia della scienza e della tecnica 
dall’inizio dell’ottocento ad oggi. 

Sfogliando le pagine di questo volume 
dedicate alle pubblicazioni riguardanti le 
varie scienze si può notare innanzi tutto 
come le opere di fisica compaiano più 
tardi (1839) che non quelle di altre 
scienze quali ad esempi la matemati- 
ca (1817). Questo probabilmente perchè 
allora. scarsi erano ancora gli interessi 
culturali e pratici connessi a questa nuova 
scienza che riuscì durante l’ultimo cin- 
quantennio a raggiungere e sorpassare le 
altre discipline. 

A tale osservazione se ne aggiunge 
un’altra relativa al notevolissimo incre- 
mento della cultura, e quindi della edi- 
toria americana, sia sotto l’influsso del- 
l'imponente processo di industrializzazio- 
ne del paese sia per lo stimolo della nu- 
merosa emigrazione intellettuale europea 
dell’ultimo trentennio. 

I volumi ricordati più interessanti 
per il fisico sono infatti per lo più di 
recente pubblicazione e ancora largamen- 
te usati, quali l’opera di LoeB su le 
scariche nei gas e i compendi di mecca- 
nica statistica del MAYER e di statistica 
per la fisica del Linpsay. Per la fisica 


10 
~ 
and 


nucleare bastera ricordare il volume di 
BLATT e WEISSKOPF e quello edito da 
EMILIO SEGRÈ. 

Resta infine da notare quel classico 
della teoria dell’informazione che è la 
cibernetica del WIENER. 

Spiace tuttavia che non compaiano in 
questa rassegna quadri statistici che docu- 
mentino in modo più accurato il lettore 
sul numero di opere pubblicate nelle di- 
verse discipline nei vari periodi. 


A. ALBERIGI QUARANTA. 


R. COURANT — Vorlesungen über Dif- 
ferential- und Integralrechnung; 
Zweiter Band (Fuktionen Mehre- 
rer Veränderlicher), Dritte Au- 
flage, pp. 468; Springer-Verlag, 
Berlin-Gôttingen-Heidelberg 1955. 


Questa terza edizione ripresenta il 
noto e ottimo trattato del COURANT con 
leggere modifiche sulla seconda edizione 
tedesca. In esso lA. tratta il calcolo 
differenziale ed integrale a varie varia- 
bili con notevole chiarezza di linguaggio 
ed esattezza dei concetti. Questi indubbi 
pregi sono dovuti anche alla mancanza 
di sottigliezze troppo formali: così il 
COURANT non trascina ovunque le n varia- 
bili ma tratta bensì con gran cura il 
caso di due variabili estendendo poi i 
concetti ed i risultati al caso più gene- 
rale con assoluto rigore. 

Nel primo capitolo lA. dà alcuni ele- 
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menti di geometria analitica e calcolo 
vettoriale. Nel secondo capitolo introduce 
i concetti di funzione di varie variabili, 
continuità e limite, per poi definire le 
derivate parziali. Nel resto del 2° capi- 
tolo e nel 3° sviluppa il calcolo differen- 
ziale, applicando concetti geometrici con 
particolare riguardo. Nei capitoli 4° e 5° 
tratta il calcolo integrale, curando spe- 
cialmente gli integrali curvilinei e super- 
ficiali e quindi i teoremi di Gauss, Stokes 
e Green; considera anche gli integrali 
impropri facendone poi applicazione agli 
integrali di Fresnel, a quelli di Fourier 
ed alla funzione gamma. Nel capitolo 7° 
introduce brevemente le equazioni diffe- 
renziali (ordinarie ed alle derivate par- 
ziali) soffermandosi maggiormente su al- 
cuni casi di particolare interesse. 

Questo testo è da raccomandare cal- 
damente, soprattutto agli studenti di 
Fisica ed Ingegneria che desiderino una 
solida preparazione matematica; infatti, 
oltre ai pregi già elencati esso ha quello 
di riportare numerose applicazioni di 
interesse fisico. Così, per esempio, lo 
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studio del calcolo degli errori, il con- 
fronto tra le rappresentazioni di Euler 
e Lagrange del moto dei fluidi, l’intro- 
duzione dei momenti ed applicazione a 
problemi fisici come quello del pentodo 
composto, le applicazioni delle equazioni 
del moto introducendo concetti di equi- 
librio e stabilità, lo studio del problema 
del potenziale, dell'equazione delle onde 
e delle equazioni di Maxwell nel vuoto ecc. 
Queste altre applicazioni inserite nel 
testo fanno sì che i concetti matematici 
vengano continuamente adoperati e quin - 
di chiariti, rendendo inoltre più interes- 
sante e piacevole lo studio. 

Bisogna però notare che è un peccato 
che non siano stati aggiunti i capitoli 7° 
e 8° dell’edizione inglese (calcolo delle 
variazioni e funzioni di variabile com- 
plessa) e gli interessanti problemi da ri- 
solvere che si trovano alla fine di ogni 
paragrafo di quella edizione. 

Numerosi ed utili i diagrammi e le 
figure, ottima la presentazione tipogra- 
fica. 
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